Physics of Strongly coupled Quark-Gluon Plasma 

Edward Shuryak 

Department of Physics and Astronomy 

University at Stony Brook,NY 11794 USA 

shuryak@tonic . physics . sunysb . edu 

September 15, 2008 

oo 
O 
O 

^ Abstract 

Oh 

J^ This review cover our current understanding of strongly coupled Quark-Gluon Plasma (sQGP), 

especially theoretical progress in (i) explaining the RHIC data by hydrodynamics, (ii) describing 
^_^ lattice data using electric-magnetic duality; (iii) understanding of gauge-string duality known as 

AdS/CFT and its application for "conformal" plasma. In view of interdisciplinary nature of the 
fH subject, we include brief introduction into several topics "for pedestrians". Some fundamental 

Q^ questions addressed are: Why is sQGP such a good liquid? What is the nature of (de) confinement 

o, and what do we know about "magnetic" objects creating it? Do they play any important role 

1^ in sQGP physics? Can we understand the AdS/CFT predictions, from the gauge theory side? 

Can they be tested experimentally? Can AdS/CFT duality help us understand rapid equilibra- 
tion/entropy production? Can we work out a complete dynamical "gravity dual" to heavy ion 
collisions? 



^ 



> 
m 
m 1 Introduction 

O Soon after discovery of QCD, people used asymptotic freedom to argue that very hot/dense matter must 
Q be weakly coupled [Ij and thus deconfined. My own entrance to this field started with the question: 
while vacuum fiuctuations of gauge field lead to anti- screening, the famous negative coefficient leading 
to asymptotic freedom, what their thermal fiuctuations would do? Explicit calculation |2], using the 
Coulomb gauge, have produced positive sign of the Debye mass, opposite to that of virtual gluons and 
the same as in QED! Thus I called this phase of matter "plasma", putting it even in the paper titlq^ 
One more important finding of that calculation was that static magnetic screening is absent, also like 
in QED: we will see that this conclusion would be seriously modified non-perturbatively. 

The next 25 years theory of QGP was based on pQCD. The program to do so, with a number 
of "signals" was proposed in my other paper j3] : those were based on high energy hadronic/nuclear 
collisions and included several perturbative process, from gluonic production of new quark fiavors to 
charmonium dissociation by gluonic "photoeffect" . A lot of efforts have been made to derive perturbative 
series for finite-T thermodynamics, and eventually all calculable terms have been calculated, see e.g. 
[1]. Convergence was bad, but we thought that some clever re-summation can still make it work. 



> 

X 



^A small anecdote is related to that: in those days plasma research in Russia was semi-classified and thus the paper 
returned to me with a note saying that it lacks proper permissions. I found a quick fix to this problem: a letter from 
renown plasma physicist saying this plasma has nothing to do with "real plasma physics". And perhaps it was true: I 
got an invitation to speak about it at some plasma physics meeting only 30 years later. 
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But whether QGP can or cannot be created experimentally was not at all clear. In fact all the 
way to the year 2000, when RHIC started, most theorists argued that nothing else but a "firework of 
mini-jets" can possibly be seen at RHIC. 

And yet, in my talk at QM99 I predicted based on hydro calculations, that elliptic flow at RHIC 
would be twice larger than at SPS. It did not take long after the start of RHIC operation to see that this 
is indeed what happens: in fact a rather perfect case for hydrodynamical explosion was made both from 
radial and elliptic flows. After lots of debates, this period culminated with the "discovery" workshop 
of 2004 and subsequent "white papers" from 4 experimental collaborations which documented it. 

Theory of QGP is still profoundly affected by this "paradigm shift" to the so called strong-coupling 
regime. We are still in so-to-say non-equilibrium transition, as huge amount of physics issues required 
to be learned. Some came from other fields, including physics of strongly coupled QED plasmas and 
trapped ultracold gases with large scattering length. String theory provided a remarkable tool - the 
AdS/CFT correspondence - which related heavy ions to the the fascinating physics of strong gravity and 
hlack holes. Another important trend is that transport properties of QGP and non- equilibrium dynamics 
came to the forefront: and for those the Euclidean approaches (lattice, instantons) we used before is 
much less suited than for thermodynamics. All of it made the last 5 years the time of unprecedented 
challenges. 

Because the issues we discuss incorporate several fields of physics, some introductory parts of this 
review are marked "for pedestrians". Indeed, heavy ion physics did not have much in common with 
string theory and black holes, or dilute quantum gases, so some basic definitions and main physics 
statements (made at an "intuitive" level) may be helpful to some readers. We start with two such 
introductory subsections, about classical strongly coupled plasmas and quantum ultracold gases, which 
we will not discuss in this review in depth. 

1.1 Strongly coupled plasmas for pedestrians 

By definition, plasmas are states of matter in which particles are "charged" and thus interact via long 
range (massless) gaugqj fields. This separate it from "neutral" gases, liquids or solids in which the 
interparticle interaction is short range. Sometimes plasmas were called "the 4-th state of matter", but 
this does not comply with standard terminology: in fact plasmas can themselves be gases, liquids or 
solids. 

Classical plasmas are of course those which does not involve quantum mechanics or h. Let me start 
with counting the parameters of the problem. There are 4 variables [j - the particle mass and density, 
the temperature and the Coulomb charge. Three of them can be used as units of mass, length and time: 
thus only one combination remains. The standard choice is the so called plasma parameter, which can 
be loosely defined as the ratio of interaction energy to kinetic energy, and is more technically defined as 

r = {Ze)y{awsT) (1) 

where Ze,aws,T are respectively the ion charge, the Wigner-Seitz radius awT = {S/AirnY^^ and the 
temperature: this form is convenient to use because it only involves the input parameters, such as the 
temperature and density, while average potential energy is not so easily available. The meaning of it 
is the same, and the values of all observables - like transport coefficients we will discuss below - are 
usually expressed as a function of F. 



^ Why only gauge and not scalar fields? Indeed, supersymmetric models have massless scalars which in many cases 
create the so called BPS situation, in which gauge repulsion is canceled by scalar attraction: and we will call them 
plasmas as well. However this is as far as it goes: a generic massless scalar, attractive in all channels and not restricted 
by supersymmetry, is just a recipe for instability and should not be considered at all. 

■^Recall that the problem is not only classical but it is also nonrelativistic: thus no Ti or c. 



Depending on magnitude of this parameter F classical plasmas have the following regimes: 
i. a weakly coupled or gas regime, for F < 1; 
ii. a liquid regime for F ?^ 1 — 10; 
iii. a glassy liquid regime for F ?« 10 — 100; 
iv. a solid regime for F > 300. 

Existence of permanent correlation between the particles is seen in the simplest way via density- 
density correlation functions 
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with A^ is the number of particles, Xi{t) is the position of the i-th-particle at time t. G{r, t) characterizes 
the likelihood to find 2 particles a distance r away from each other at time t. Here are some examples, 
from our own (non-Abelian) MD simulations |S|, which show that liquid regime demonstrate nearest- 
neighbor peaks, and crystals have peaks corresponding to longer range order. They also show "healing" 
of correlations with time in gases, but much less so in liquids and solids. 
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Figure 1: (Color online) Gd correlation function for F = 0.83,31.3, 131, respectively. Red circles corre- 
spond to t* = 0, and blue squares correspond to t* = 6. 



The case of small F is widely discussed in statistical mechanics courses: let me just remind the reader 
that it is in this case when one can use Boltzmann eqn, cascades and other simple tools appropriate 
for a gas. Unlike gases and solids, the interplay of local order and randomness at large distances makes 
liquids difficult to treat theoreticalljQ Thus, in spite of their crucial importance for a lot of chemistry 
in general and our life in particular, most physics and statistical mechanics courses tend to either omit 
them completely or tell as little as possible about them. It is possibly worth reminding heavy ion 
practitioners, that for liquids neither Boltzmann equation nor cascades can be used because particle are 
strongly correlated with several neighbors at all times. The very idea of "scattering" and cross section 
involves particles coming from and going to infinity: it is appropriate for dilute gases but not condensed 
matter where interparticle distances do not exceed the range of the forces at any time. 

''I heard an opinion, ascribed to a lecture of V.Weisskopf, that if theorists would invent the Universe from scratch, 
without any experiment, they would never think about liquids. 



Strongly coupled classical electromagnetic plasmas can be studied experimentally: they are not at 
all exotic objects. For example, table solt NaCl can be considered a crystalline plasma made of perma- 
nently charged ions Na^ and Cl~ . At T ~ lO^/f (still too small to ionize non- valence electrons) one 
gets a molten solt, which is liquid plasma with F ~ 60. A more famous object of recent experimentation 
is a charged dust: in space (e.g. at International Space Station) it has been put into a nice crystal and 
studied in depth. For example, one can get a particle piercing it and creating Mach cones, several if 
there are excitations other than the sound. One can find more introductory information and references 
on the subject in Mrowczynski and Thoma review [Ej. 

We will not discuss any theory of it, but just note that starting from about 1980's availability of 
computational resources get sufficient to use "Molecular Dynamics" (MD). This means that one can 
write equations of motion (EOM) for the interparticle forces and directly solve them for say 10^ particles. 
It is this simplest but powerful theoretical tool we will use to access properties of strongly coupled QGP 
by classical simulations. 

1.2 Strongly coupled ultracold gases 

Back in 1999 G.Bertsch formulated a "many-body challenge problem," asking: what are the ground state 
properties of a two-species fermion system in the limit that the scattering length (a^) of its interaction 
approaches infinity? (This limit is usually referred to as the "unitary limit" because the scattering cross 
section reaches its unitarity limit per s-wave.) Such problem was originally set up as a parameter-free 
model for a fictitious dilute neutron matter: recall that nn scattering length is indeed huge because of 
near-zero-energy isoscalar virtual state. 

The answer was provided experimentally by atomic experimentalists, who found a way to modify the 
interaction strength in ultracold trapped systems by the magnetic field, which can shift level positions 
till they cross and form the so called "Feshbach resonances". As a result, the interaction measure - 
the scattering length in units of interparticle spacing a * n^/"^ - can be changed in a wide interval, 
practically from — cx) (very strong attraction) to +oo (very strong repulsion). Systems of ultracold 
atoms became a very hot topic a decade earlier, when laser cooling techniques were developed so 
that temperatures got low enough (so that atomic thermal de Brogle wavelength get comparable to 
interparticle distances) to observe Bose- Einstein condensation (BEG) for bosonic atoms. Last years was 
the time of strongly coupled quantum gases made of fermions, possessing superfiuidity and huge pairing 
gaps. Never before one had manybody quantum systems with widely tunable interparticle interaction 
amenable to experimentation, and clearly it created a kind of revolution in quantum manybody physics. 

For some technical reasons, BEG in strongly systems (of interest to us in respect to monopole 
condensation/confinement) is not yet studied, while strongly coupled fermions have been investigated 
quite extensively. Let me briefiy review the questions which were of the main interest of the atomic 
community. One of them is whether two known weakly coupled phases - BGS superconductor for small 
and negative a and BEG of bound atom pairs at small positive a - join smoothly or there exist a 
discontinuity at the Feshbach resonance. The answer seems to be the former, namely two phases do 
join smoothly. Another issue is whether low-T strongly coupled fermions is in a gaped superfluid state: 
it has been answered positively, and phenomena as complex as Abrikosov's lattice of vortexes has been 
observed. 

Glearly, there are many fascinating phenomena in this field, but let us focus on just two issues most 
relevant for this review and sQGP: 

(i) whether transition from weakly to strongly coupling regime is refiected in an unusually small viscosity, 
seen in onset of hydrodynamical behavior in unexpectedly small systems; 

(ii) possible universal limits in the infinite interaction limit a * n^^^ -^ oo for pairing, in connection to 
color superconductivity of quark matter. 

Rather early it has been observed that the answer to the former question is clearly affirmative. When 
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Figure 2: (a) False color absorption images of a strongly interacting degenerate Fermi gas of ultracold 
^Li atoms as a function of time after release from a laser trap. From O'Hara et al.ff] (b) "Quantum 
viscosity" in strongly-interacting Fermi gas a = rj/hn (trap-averaged), (c) Same data for the shear 
viscosity as rj/hs in units of the entropy density s as a function of energy E. The lower green dotted 
line shows the string theory prediction l/(47r). The light blue bar shows the estimate for a quark-gluon 
plasma (QGP) while the blue solid bar shows the estimate for ^He, near the A- point. 



the trap is switched off, the atoms in a weakly coupled gas simply fly away with their thermal/quantum 
velocities, displaying isotropic angular distribution of velocities irrespective of the trap shape. However 
in strongly coupled regime, with the particle mean free path smaller than the system size, hydrodynam- 
ical flow develops. Deformed traps thus develop "elliptic flows" in the direction of maximal pressure 
gradient, in a way analogous to what is happening in heavy ion collisions: see FigMa). This is by no 
means trivial: the interparticle distances are about 1000 times the atomic size, and the total number of 
atoms is only ~ 10^ (only few times more than in central heavy ion collisions at RHIC): similar number 
of water molecules would not show any hydro! 

The next questions was whether one can use hydrodynamics quantitatively, find out its accuracy 
and quantify the viscosity. From "released traps" the experiments switched to quadrupole vibrational 
modes: as the trapped system has a cigar-shape with much weaker focusing along z axes compared 
to axial ones, there is softer z-vibrations and higher frequency "axial" mode. I will not go into vast 
literature and simply say that the value of the vibrational frequencies are indeed given by hydro, reaching 
near-percent accuracy at the Feshbach resonance point (where \a\ = oo). 

The first study of viscosity has been made by Gelman, myself and Zahed [8] : from available data on 
two different vibrational modes - z-mode and axially symmetric radial mode of a cigar-shaped atomic 
cloud - we tried to deduce viscosity and found that the values are roughly consistent with each other. 
Instead of going into details of that work, let us discuss what one would expect based on "universality" 



arguments. The main point is that if certain observable are finite in the a * n^^^ — > cxo hmit, the value 
can only depend on few parameters - e.g. the particle mass m, the density n and Plank constant h when 
the temperature is zero. By dimensional analysis, at T = 0, the energy density of infinitely strongly 
coupled gas can only be a number times the energy density of the ideal Fermi gas with the same density, 

e(n) = /3 (3) 

m 

with some universal constant (3. 

Assuming that there can be nonzero "quantum viscosity" at zero T for large amplitude oscillations, 
we [8] proposed to measure viscosity in units of a = rj/hn, the only combination with the right dimen- 
sion. Experimental data by Turlapov et al. [9] for viscosity in such units are shown in FigMb). It is 
plotted as a function of energy per particle in Fermi energy units, E/Ep, which basically characterizes 
the excitation temperature: E/Ep !^ .5 (corresponding to the left side of the figure) is close to the 
ground state T = 0. From the measured points it seems that a = rj/hn is actually vanishing in the 
ground state, about linearly in T. Another way to express viscosity - familiar from Black Hole physics- 
is to express it as rf/fis with s being the entropy density. Such ratio is shown in FigMc): and it looks 
like this one reaches nonzero value as T, s — > 0. Its magnitude makes strongly coupled fermionic atoms 
to be the "second best liquid" known, in between the sQGP (light blue band below) and the "bronze 
winner" (former champion) , liquid ife4, shown by dark blue. 

Few people (myself included) may find it fascinating, but most of atomic quantum gases commu- 
nity do not care and even unaware that some qualitatively new regime happens close to the Feshbach 
resonance: global parameters (like (3 in total energy which was subject of many theoretical and nu- 
merical works) are smooth there. Small viscosity of quantum gases does not yet have any microscopic 
explanations: certainly not by "unitary" cross section and kinetic theory. Linear behavior of rj and s 
in T reminds that of electrons in solids: but strong coupling actually destroys the Fermi surface (as 
experiments measuring momentum distribution show quite clearly) and makes a very strong supercon- 
ductor. The critical point in units is at E/Ep ~ .85: a look at the discussed figure reveals no changes 
in viscosity (nor in oscillation frequencies themselves) visible by an eye. 

Apart of being in general related to the issue of sQGP as "the most perfect liquid" , the strongly 
coupled atomic superconductor provides some valuable information on how large the pairing gaps can 
possibly become in cold quark matter with color superconductivity. Assuming that near deconfinement 
cold quark matter is also strongly coupled, with weakly bound diquarks playing the role of Feshbach 
resonance, in [TQ] I have used universality of the ratio T^/ep (the pairing critical temperature to the 
Fermi energy) and data on strongly coupled fermionic atoms to get upper limit on the corresponding 
transition to the color superconductivity. Fig|3] explains the argument. "Universality" tells that the 
critical temperature must be simply proportional to the Fermi energy 

T, = ar^Ep (4) 

with the universal constant ap^. In fact two values of T show some change: the phase transition to 
superfiuidity at ap^ = J^ = .35 and 02 = |^ ~ 0.7 — 0.8 which experimentalists (Kinast et al) interpret 
as a transition to a regime where not only there is no condensate of atomic pairs, but even the pairs 
themselves are melted. Using these two values as slopes of solid and dashed lines in Fig|3]^b), one can 
see that all chemical freezeouts in heavy ion collisions (points) are above possible domain of pairing, 
even at infinite coupling. At finite coupling the gaps and Tc is even smaller. 

Let me finish with one recent development, providing quite fascinating prospects if true. The idea 
is that infinitely strongly interacting atoms may have "gravity dual" description similar to AdS/CFT. 
Son fTT] and also Balasubramanian and McGreevy [12] have argued that since at a * n^/^ — > c)0 the 
interaction can be treated as a kind of boundary condition when two atoms meet, and this problem has 
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Figure 3: (a) Schematic phase diagram for QCD, in the plane baryon chemical potential - temperature. 
M (multifragmentation) point is the endpoint of nuclear gas-liquid transition. E is a similar endpoint 
separating the first order transition to the right from a crossover to the left of it. (Black) solid lines show 
phase boundaries, dashed lines are curves of marginal stability of indicated states. Two dash-dotted 
straight lines are related with bounds from atomic experiments we discuss in the text, they intersect with 
unbinding of diquark Cooper pairs (D) and most strongly coupled point (S), which is at the maximum of 
the transition line and is also a divider between BCS-like and EEC-like color superconductor. Other lines 
are zero binding for singlet gg, singlet cc, octet qg, gg and finally triplet qq states, (right) Compilation 
of experimental data on the chemical freezeout from different experiments: squares (circles) are for fits 
at mid-rapidity (all particles), respectively. Two solid lines are the phase transition lines with the quark 
effective mass Mi = and 100 MeV, two dashed lines show pair unbinding lines for the same masses. 

what he called "Schreodinger symmetry" , a subgroup of conformal symmetry. They found examples of 
the metric possessing this very symmetry and suggested it may be the desired "gravity dual" . (The 
reader should however be reminded that there is no proof of its existence or in fact any empirical 
confirmation of this idea so far.) 



2 Heavy Ion Collisions 

2.1 Heavy ion collisions and flows for pedestrians 

Let me start with brief history. At one hand, high energy physics was for all of its history been interested 
in "high energy asymptotic" s -^ oo. In practice this mostly was reduced to global features like total 
and elastic cross sections, as it was thought that multibody final state produced is too complicated to 
get some sufficiently simple theoretical treatment. Nuclear physicists started from nuclear collisions at 
nonrelativistic domain: their cross section was obvious and the main objective was understanding of 
"excited matter" . Experimental program at Relativistic Heavy Ion Collider (RHIC) at Brookhaven and 
Large Hadronic Collider (LHC) at CERN is so to say a brainchild of both communities. 

What do we mean saying that collisions produced "excited matter" ? Should one expect simplification 
if it is the case? The answer is that we are not interested in any excited system produced (e.g. in 
"elementary" pp collisions), but mostly in a macro scopically large fireball whose size (macro scale) R 
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greatly exceed the micro scale I of correlations inside it. 

R » I (5) 

In weakly coupled systems (gases) / is the mean free path length, or relaxation time. In strongly coupled 
setting such as AdS/CFT the temperature T would be the only dimensional parameter describing the 
microscopic physics: most dissipative phenomena have the scale 1/T as a relaxation scale. 

Have we reached this regime? Perhaps good illustration that the question is nontrivial are my 
own two unsuccessful attempts to apply macroscopic physics for high energy processes. In 1971 I 
proposed "spherical explosion" in e~^e~ annihilation into many pions [T3], only to be killed by discovery 
of asymptotic freedom in 1973 and jets in e~^e~ evens in 1976. In 1979 Zhirov and myself ^^ looked 
at fresh results from pp collisions at then-new ISR collider at CERNJj The general idea of the paper 
was to look for collective transverse flows. We argued that since secondaries have different masses, the 
kinematic effect of a collective motion (flow velocity) can be separated from their thermal spectra. More 
specifically, if the "matter flow" is only longitudinal, along the beam direction, the transverse momenta 
spectra of different secondaries would be just thermal 

dN , My 



rsj 



dvl '^^^"^^ ^'^ 

where the mass and transverse momentum combine into a "transverse mass" M^ = M^ +Pt- The ISR 
data indeed showed this "M^ scaling" : but there was no sign of the deviations from it due to transverse 
flow. (Lacking the effect, we argued that it is the "vacuum pressure" which stops it.) 

Heavy ion collisions are used to make the system larger. The macro scale - nuclear sizes - for Au 
(or Pb) used in Brookhaven (or CERN) have typical A ~ 200 and radius R ~ 6/m. At RHIC the 
temperature T changes from the initial value of about Tj ~ AOOMeV ~ l/{.5fm), thus RT ~ 10 >> 1 
can be viewed as a large parameter. In fact experiments with Cu beam, A = 64, also showed good 
macroscopic behavior, as are rather peripheral collisions. The boundary between macroscopic and 
microscopic systems should be somewhere, as pp collisions are not hydro-like: but the exact location 
of it is not yet known: small systems of several nucleons are subject of large fluctuations, and it is not 
easjjjto study them. 

Once produced, the fireball is expanding hydrodynamically up to the so called freezeout conditions. 
People familiar with Big Bang can recognize existence of multiple freezeout, for each reaction at its 
own time: and in the Little Bang it is similar. The so called chemical freezeout (at which particle 
composition gets frozen) is at T ~ 170 MeV at RHIC, while re-scattering continue to Tf ~ 90MeV 
(for central collisions). In units of QCD critical temperature, it is T/Tc variation from about 2 to 1/2, 
conveniently bracketing the QCD deconfinement/chiral restoration transition. 

Flow velocity is decomposed into longitudinal (along the beam) and transverse components. The 
latter is further split into radial flow (present even for axially symmetric central collisions) and elliptic 
flow which exists only for non-central collisions. The reason elliptic flow to be very important was 
pointed out by H.Sorge: it is developed earlier than the radial one, and thus most sensitive to the QGP 
era than other flows. 

The radial flow has Hubble-like profile, with transverse rapidity growing roughly linear till the edge 
of the fireball. The maximal radial flow velocity at RHIC is about .7 c. This radial flow is firmly 



^ L. van Hove and few others, including myself, put forward a proposal to accelerate heavy ions in ISR, to see what 
happens. CERN leaders basically rejected it; the experiments never went beyond the alpha-particles, and then this first 
hadronic collider was physically destroyed. As we know now, QGP could have been discovered and studied at ISR 20 
years prior to RHIC. 

^On top of that, there is experimental background, as very peripheral collisions of two beams is hard to tell from 
beam-residual gas interaction. 



established from a combined analysis of particle spectra, HBT correlations, a deuteron coalescence and 
other observables: so we will not show 

At non-zero impact parameter the original excited system has a deformed almond-like shape: thus 
its expansion in the transverse plane can be described by the (Fourier) series in azimuthal angle 

Vn{s,puMi,y,b,A) =< cos{n(f)) > (7) 

where average is over all particle^ At mid-rapidity y = only even harmonics are allowed, the second 
one characterizes the so called elliptic flow. Multiple arguments of the parameter V2 stand for the 
collision energy s, transverse momentum pt, particle mass/type Mj, rapidity y, centralitjjj b and the 
atomic number A characterizing the colliding system. 

If in high energy collisions of hadrons and/or nuclei a macroscopically large excited system is pro- 
duced, its expansion and decay can be described by relativistic hydrodynamics. Its history starts with 
the pioneer paper by Fermi of 1951 |T6] who proposed a statistical model to Lorentz-contracted initial 
state. Pomeranchuk ^^ then pointed out that initial Fermi stage cannot be the final stage of the 
collisions since strong interaction in the system persists: he proposed a freezeout temperature Tf ~ 171^^. 
L.D. Landau [15] then explain that one should use relativistic hydrodynamics in between those two 
stages, saving Fermi's prediction of the multiplicity by entropy conservation. 

(Before we go into details, a comment: often hydrodynamics is considered as some consequence of 
kinetic equations, but in fact applicability conditions of both approaches are far from being coincident. 
In particular, for the former approach the stronger the interaction in the system is, the better. Ki- 
netic approach, on the contrary, was never formulated but for weakly interacting systems: and as we 
repeatedly emphasize in this review, it is not so for QGP.) 

The conceptual basis of the ideal hydrodynamics is very simple: it is just a set of local conservation 
laws for the stress tensor {T^'^) and for the conserved currents (Jf ), 

df.T^" = (8) 

d^Jt = 

The assumption is in local form for T^'^ and Jf related to the bulk properties of the fluid and its 
4-velocity u'^ by the relations, 

T^"" = [t + p)u^u'' - pg^"" (9) 

Jf = niU^" 

Here e is the energy density, p is the pressure, Ui is the number density of the corresponding current, and 
u^ = 7(1, v) is the proper velocity of the fluid. In strong interactions, the conserved currents are isospin 
(Jj), strangeness (J5), and baryon number {J^)- For the hydrodynamic evolution, isospin symmetry 
is assumed and the net strangeness is set to zero; therefore only the baryon current Jb is considered 
below. 

Let me add a simple heuristic argument why the first term in the stress tensor has (e + p) and 
not any other combination. The point is e and p themselves are defined up to a constant ±B (which 
depending on context we call the bag constant, or cosmological term or dark energy). The combination 
(e + p) does not have it, and it is also proportional to entropy which is defined uniquely. How this 



^This is written as if the direction of impsict parameter of the event is known. I would not go into how it is determined 
in detail but just say that it usually comes from counters sitting at rapidity y w ±3 far from the region where most 
observations are done, at mid-rapidity y w 0. 

^In real experiments the measure of impact parameter is the so called the number of participants Np, which changes 
from zero at most peripheral collisions to 2 A at central collisions. This number is the total number of nucleons 2 A minus 
the so called spectators, which fly forward and are directly observed by forward-backward calorimeters. 
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argument complies with the last term in ([9])? Well it has B but without velocity, so in hydro eqns there 
is only pressure gradient and this B term disappears as well. 

In order to close up this set of equations, one needs also the equation of state (EoS) p(e). One should 
also be aware of two thermodynamical differentials 

de = Tds dp = sdT (10) 

and the definition of the sound velocity 

2 dp s dT 

and that e + p = Ts. Using these equations and the thermodynamical relations in the form 

Al = ^ (12) 

e+p s ^ ' 

one may show that these equations imply another nontrivial conservation law, namely, the conservation 
of the entropy 

d^^su^) = (13) 

Therefore in the ideal hydro all the entropy is produced only in the discontinuities - shock waves - 
which are not actually there is application we discuss. Thus the "initial entropy" is simply passed to 
the solution as an initial parameter, determined in the earlier (violent) stage of the collision: (this is 
similar to Big Bang cosmology, in which "entropy production" stage is also very different from stages 
of cosmological evolution we can observe by e.g. Nucleosynthesis.) 

Next order effects in micro-to-macro expansion is the domain of "viscous hydrodynamics" : we will 
discuss their applications to data description as well as their derivation from AdS/CFT settings. 

The simplest Bjorken 1+1 dim solution is a good example "for pedestrians" , reminding how to write 
hydro in arbitrary coordinates. 

T'"V = 0, J% = 0, (14) 

where the semicolon indicates a covariant derivative. For tensors of rank 1 and 2 it reads explicitly 

■i -i _i_ -pi -k rpik rpik , -pi rpmk , pk ijiim (^K\ 

J;p ~ J,p ' ^ pkJ } -'- ;P ~ -^ ,P "T -*- pm-' "T -*- pm-' ) \^'-') 

where the comma denotes a simple partial derivative and the Christoffel symbols T'fj are given by 
derivatives of the metric tensor g"-^{x)\ 

Tl^ = {l/2)g^^{g,,,j + (?,,,, - (7,,,^) . (16) 

As an example, let us do the following transformation from Cartesian to light cone coordinates: 

x^ = (t, X, y, z) — ^ x"" = (r, X, y, ri) 

t = r cosh 1] T = y/f^ — z^ 

t+z 
z = Tsmh.ri r] = (1/2) In . 

In the new coordinate system the velocity field (after inserting Vz = z/t) is given by 

tT = ^{l,%,Vy,Q) (17) 






with Vi = fjCosh?7, i = x,y, and 7 = 1/ Jl—vl 
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Figure 4: The left and right sides show the hydrodynamic solution at the SPS and RHIC The thin lines 
show contours of constant transverse fluid rapidity {vt = tanh(|/7^)) with values 0. 1,0.2,.. .,0.7 . The thick 
lines show contours of constant energy density, e^o denotes the energy density where T = 120 MeV. 
en and cq denote the energy density where the matter shifts from hadronic to mixed and mixed to 
a QGP, respectively. The shift to hadronic cascade is made at en- {yr) denotes the mean transverse 
rapidity weighted with the total entropy flowing through the energy density contours. Walking along 
these contours, the line is broken into segments by dashed and then solid lines. 20% of the total entropy 
passing through the entire arc passes through each segment. 

Now we turn to the metric of the new system. We have 



dt^ — dx^ — dy^ — dz^ 
dr"^ — dx^ — dy^ — r'^drf 



and therefore 



/ 1 Q \ 
0-10 
0-10 

V -r2 / 



(19) 



The only non-vanishing Christoffel symbols are 






T . 



The l+ld equations for boost-invariant solution can be written in the following way 

d d 

-—(s cosh-u) + 7^(s sinhu) = 



(20) 



(21) 
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- (T sinh y) + -Q-^{T cosh y) = (22) 

The so called Bjorken (TS] solutioir] is obtained if the velocity is given by the velocity u^ = 
{t,0,0,z)/T where t"^ = t^ — z^ is the proper time. In this 1-d-Hubble regime there is no longitu- 
dinal acceleration at all: all volume elements are expanded linearly with time and move along straight 
lines from the collision point. The spatial tanh~ {z/t) and the energy-momentum rapidity ytanh" v 
are just equal to each other. Exactly as in the Big Bang, for each "observer" ( the volume element ) 
the picture is just the same, with the pressure from the left compensated by that from the right. The 
history is also the same for all volume elements, if it is expressed in its own proper time r. 

Thus the entropy conservation becomes the following (ordinary) differential equation in proper time 

T 

dsir) s . . 

-f^ + - = 23 

dr r 

which has the obvious solution 

const , ^. 

s = (24) 

T 

Let us compare three simple cases: (i) hadronic matter, (ii) quark-gluon plasma and (iii) the mixed 
phase (existing if there are first order transitions in the system). In the first case we adopt the equation 
of state suggested in j22] c^ = dp/de = constij) ^ .2. If so, the decrease of the energy density with 
time is given by 

e{r) = em-f^'" (25) 



r 

In the QGP case the same law holds, but with c^ = 1/3 

In the mixed phase the pressure remains constant p = Pc , therefore 

e(r) = (e(0)+pe)(^^)-p, (26) 

r 

So far all dissipative phenomena were ignored. Including first dissipative terms into our equations one 

1 de_lds_ 1/ (4/3)r^ + C \ ^27) 

e + pdr s dr t\ {e + p)T J 

Note that ignoring ( one finds in the r.h.s. exactly the combination which also appears in the sound 
attenuation, so the correction to ideal case is (1 — Tg/r). Thus the length Tg directly tells us the 
magnitude of the dissipative corrections. At time r ~ F^ one has to abandon the hydrodynamics 
altogether, as the dissipative corrections cannot be ignored. Since the correction is negative, it reduces 
the rate of the entropy decrease with time. Another way to say that, is that the total positive sign 
shows that some amount of entropy is generated by the dissipative term. We will discuss "gravity dual" 
to this solution in the last chapter. 

2.2 Collective flows and hydrodynamics 

Any treatment of the explosion by hydro in general includes (i) the initiation; (ii) hydro evolution 
and (iii) freezeout. Geometrically, (i) and (iii) form two 3-hypersurfaces in 4 dimensions, together 
constituting a boundary of 4-volume region in which hydro eqns are solved. Only the stage (ii) deals 
with the Equation of State and transport properties of new form of matter, the QGP - while (i) 



^ We call it following established tradition, although the existence of such simple solution was first noticed by Landau 
and it was included in his classic paper as some intermediate step. The space-time picture connected with such scaling 
regime was discussed in rcfs [201111] before Bjorken, and some estimates for the energy above which the transition to the 
scaling regime were expected to happen were also discussed in my paper [3j as well. 
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and (iii) have to deal with stages at which hydrodynamics is not apphcable. The stage (i) is least 
understood theoretically: yet the uncertainties it produces are rather minor. Correct treatment of the 
freezeout stage (iii) is much more important for calculations of the final spectra which are compared to 
data: and that is where most intensive debates were. Fortunately the latest hadronic stage is not that 
mysterious, as some theorists think: it happens in a dilute gas-like medium made of mostly pions and 
their resonances. We understand their low energy scattering cross sections quite well, the corresponding 
cascade codes have been extensively tested using the low energy heavy ion collisions at AGS (BNL) and 
SPS(CERN). It just needs a bit of extra work. 

Solving hydro eqns is not simple - they are nonlinear and prone to instabilities - but solutions are 
easily controllable by the energy, momentum and entropy (for ideal hydro) conservation: thus I presume 
it was done correctly by all groups. At years 2000-2004 most groups used an approximation in which 
one out of 3-1-1 (coordinates-|-time) - the longitudinal spatial rapidity - was taken as irrelevant, thus 



switching to 2-1-1. Example of output from such hydro calculations [19], with properly chosen Eo^^° 
is shown in Fig|4| for two collision energies, s^^^ = 15,200G'ey per nucleon, marked SPS and RHIC 
and shown at left and right figure. The only thing I would like to mention is that the fraction of time 
spent in the QGP phase (below the eg curve) is not that large: but its existence crucially change flow 
pattern at RHIC. As one can see from the right plot, all lines of constant rapidity are nearly vertical 
and nearly equidistant: it means after QGP the flow gets a simple Hubble-like flow with Vt = Hrt and 
imve-independent constant H. 

Location of the freezeout surface is also a nontrivial task. In our flrst application of hydrodynamics 
for radial transverse flow at SPS [23] we developed rather tedious "differential freezeout", for each 
geometry and each secondary particle was applied. We calculated individual reaction rates for different 
secondaries and matched them with hydro expansion rates locally. In spite of relative success of the 
work, there were no followers to this approach at RHIC. Selecting hydrodynamics as a Ph.D. topic for 
my graduate student, Derek Teaney [19], we had in mind a different procedure for freezout, namely 
switching at the onset of hadronic phase {ch in the flgure above) to a hadronic cascade (RQMD) 
which automatically leads to earlier freezeout for smaller systems or particles with smaller hadronic 
The same approach was later used by Hirano et al, who generalized our rapidity- 



cross sections 
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independent hydro to the full 3-1-1 dimensional case, successfully reproducing also the V2 dependence on 
rapidity y [23]. They confirmed that V2 linearly decreases from its maximum at mid-rapidity linearly, for 
the same reason as it decreases toward smaller energies: there is less matter there, shorter QGP phase 
and also shorter hadronic phase due to earlier freezeout. Independently developed code by Nonaka 
and Bass [25], with a different cascade code UrQMD, later also confirmed these calculations and well 
reproduce data for all dependences. These three groups have basically covered all outstanding issues 
related to applications of the ideal hydrodynamics to RHIC data. 

One might think that, after a couple of groups checked the calculations themselves, the rest would be 
history if all heavy ion community would recognize/accept it. Unfortunately, it is not the case. Many 
hydro groups have not implemented hadronic freezeout, ending hydro at some (arbitrarily selected) 
isoterms T = Tf. There is absolutely no reason to think it is the right choice. Even thinking about 
freezout as a local concept, one finds that it is determined not by local density but by local expansion 
rate n = d^Un of matter. It is well known for more than a decade [^ that the isoterms do not resemble 
the lines of fixed expansion rate. Furthermore, while hydro solutions simply scale with the size/time 
of the systemPl the freezout conditions (involving the reaction rates) do not. It is thus not at all 
surprising, that many results based on unrealistic freezeout show quahtatively wrong dependencies. 

Let me start with few plots from PHENIX "white paper" [26] in Figjs] for protons and pions, to 



^°LH8 mentioned in the figure means the "latent heat" of magnitude .8GeV/fm^ at the QCD transition temperature. 
^^Note that e.g. meson and the nucleon have about the same mass of 1 GeV, but cross section of re-scattering on 
pions ranging from few mb for the former to 200 mb for the latter, at the A resonance peak dominating re-scattering. 
^^Hydro eqns have only first derivative over coordinates, which can thus be rescaled. 
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Figure 5: Top two panels: On the left, proton V2{pt)/^ vs. px for minimum-bias collisions at RHIC 
[2Zl EH] are compared with hydro calculations [121 [221 EDI EI] , and on the right is the same comparison 
for pions. Bottom two panels: {l/27i)(fN/ptdptdy. On the left, for protons, for 0-5% centrality bin 
collisions at RHIC [32] are compared with the same hydro calculations. On the the right, the same 
comparison for pions. 



illustrate some important physics points explaining different level of success of different authors. The 
lower part of FigjS] shows measured pt spectra of protons and pions, in comparison with different hydro 
calculations. The shape of those is very different mostly because heavy protons and light pions have 
different thermal motion at the time of freezout, in spite of the same collective radial flow. Note that 
nearly every group has a correct shape of the spectra (and thus correct radial flow velocity), but they 
don't aways have the normalization (for the nucleons) correctly: it is because of "chemical freezout" 
not implemented by some groups. There is no problem for hydro+cascade model [T9] (the red curve). 

Now we switch to elliptic flow, shown as the ratio of observed momentum anisotropy V2{pt) divided 
by calculated spatial elliptic anisotropy e, shown in the upper part of FigjSJ Although some calculations 
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are not too close to the data, the overall magnitude of the effect and its pt dependence is clearly 
reproduced. However, that is only true for the "good" dependence V2{pt)- (Actually there is another 
one, the dependence of V2{'m) on the particle mass, which everybody get right. See one nice example of 
that from Hirano in Fig. [6Fd).) 

Unfortunately other dependences of V2 are not so forgiving as V2{pt) and they show qualitative 
differences between models which do and do not include hadronic freezeout properly. Those include the 
dependences of the elliptic flow on (i) collision energy s^/^, (ii) centrality b or number of participants A''^, 
(iii) and rapidity y. Let me start with the collision energy: Heinz and Kolb in their large hydro review 
|30j give their excitation curve for elliptic flows shown in FigJ6](a). All variants of their prediction 
for the elliptic flow has rapid rise on the left side of the plot (at low collision energies) with about 
constant saturating values at higher energies (one variant even reaches a peak followed by a decrease 
). As one sees from see FigMb) from the same review, this is not the trend observed in the RHIC 
domain: the data show a slow rise without peaks or saturation. Kolb and Heinz thus concluded that 
hydro is not supported by the data, at all collision energies below RHIC. This lead to a myth about 
a "hydro limit" which was "never reached before RHIC" which was (and still is) repeated from one 
conference to another. Finally FigMc) from that review display rapidity dependence of V2{y) from a 
calculation by Hirano (before he switched to hydro+cascade): the conclusion was that hydro only works 
at mid-rapidity. 

All those results are for fixed-T freezeout, which is not based on anything and thus is simply wrong. 
Here what hydro+cascades approach finds for all of these observables. The energy excitation curves of 
V2 from [T9] and [23] are shown in Fig. [tI left and right. When correct freezout is implemented, the 
elliptic flow is rising steadily all the way from SPS to RHIC, as the data do. The reason Heinz et al 
(as well as curved marked 120 Mev and 100 MeV in left and right plots) strongly overshoot the data 
is simply because the freezout does not occur at the same T at different collision energy. In fact, it is 
independently measured (from radial flow for central collisions) that while the freezeout temperature 
at SPS is about 140 MeV, it is as low as 90 MeV at RHIC. The trend is well understood: the larger is 
the system, the hotter it is at the beginning, and the cooler it gets at the end of the explosion! 

Another "hydro problem" discussed in Kolb- Heinz review, the rapidity dependence of V2, also went 
away as soon as correct freezeout was used by Hirano et al. The results (circles in Figjs] ) are right on 
top of the data (triangles), without any change of any parameters. The reason is exactly the same as for 
the energy dependence: in fact one can check that ^2(2/) show good "limiting fragmentation properties, 
depending basically on y — Y, the distance to beam rapidity Y. One can see the difference in centrality 
dependence as well, in the left side of Fig|8] . 

Intermediate summary: these "problems" (and associated myths) were caused by wrong freezeout. 
Matter at this time is a dilute pion gas, which is not a good liquid, neither at SPS, not at RHIC and 



will not be at LHC, and cascade is the besl^ approach we have to describe it. As far as we can now 
test, two other evolution eras - sQGP and "mixed" or near-T^ one - can be surprisingly well described 
by the ideal hydrodynamics. 

The last statement should not be understood that the agreement with all details of the data is perfect. 
Theoretically, one should always ask about accuracy and applicability limits of this hydrodynamical 
description. As emphasized by Teaney |135j the answer should be obtained by calculation the role of 
viscous effects. Since we will have extensive derivation of "derivative expansion" in AdS/CFT language 
at the end of this paper, I will not go into details here, going directly into new developments. 

In order to get more accurate account of viscosity effect on flow, a new round of studies has been 
performed during the last year. Relativistic Navier-Stokes has some problem with causality, thus 'higher 
order" methods has been used. Apart from viscosity those methods have another parameter, the 

^^ This does not imply that we have complete confidence in many details of those cascades. To name one outstanding 
issue: the precise in-matter modification of hadronic resonances like p, A etc, dominating the cross sections, is being 
addressed but still far from been solved. 
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Figure 6: (a) Excitation function of the elliptic (solid) and radial (dashed) flow for Pb+Pb or Au+Au 
collisions at 6 = 7 fm from a hydrodynaniic calculation. [33] The collision energy is parametrized on 
the horizontal axis in terms of total particle multiplicity density dN/dy at this impact parameter, (b) 
A compilation of V2 data vs. collision energy from mid-central (12-34% of the total cross section) 
Pb+Pb and Au+Au collisions, (c) p^-integrated elliptic flow for minimum bias Au+Au collisions at 
s^/^ = 130 A GeV as a function of pseudorapidity, compared with data from PHOBOS and STAR. 
(d)Transverse momentum dependence of V2 for pions, kaons and protons. Filled plots are the results 
from the hybrid model. The impact parameter in the model simulation is 7.2 fm which corresponds to 
20-30% centrality. 

relaxation time r^, which is can either be used as a regulator - and its value put to zero at the end - or 
as a real representation of two-gradient terms. There are 4 groups who have reported solving 2+1 dim 
higher order hydrodynamics. P. and U. Romatschke [SB], Dusling and Teaney [2Z| , Heinz and Song^^ 
|38j . Molnar [39] have compared viscous hydro with some version of his parton cascade and found good 
agreement when the parameters are tuned appropriately: but cascade describe f2(pt) even at larger 
momenta, the solution is supposed to converge to that of the Navier-Stokes eqn, avoiding the causality 
problems. 

In Fig|9]we show Romatschke's results: literally taken they favor very small viscosity, even less than 
the famous lower bound. Now, is the accuracy level really allows us to extract rj/s'! The uncertainties in 
the initial state deformation [iQl UT] are at the 10% level, comparable to the viscosity effect itself. EoS 



^''This group originally found viscosity effects about twice larger than others: but it was found in their later work that 
this happened because of different account for some higher order term. In the t,. —>■ limit, in which Navier-Stokes limits 
is supposed to be recovered, all results are now consistent with others. 
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Figure 7: ( Two examples of the excitation function of V2 in Au+Au collisions: (a) from [19] and (b) from 
|24j . Points obtained with fixed-Tj (120 Mev in (a) and 100 MeV in (b)) freezeout strongly overpredict 
the data at low energies, but with cascade freezeout (red curves in both) those describe data quite well. 
(Two triangular points in (a) correspond to SPS and RHIC collision energies.) 
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can probably be constrained better, but I think uncertainties related to freezeout - not yet discussed at 
all - are also at 10 percent level, although they can also be reduced down to few percent level provided 
more efforts to understand hadronic resonances/interactions at the hadronic stage will be made. All of 
it leaves us with a statement that while literally fits require 77/5 ~ 0.1 or less, we can only conclude 
that it is definitely 77/5 < .2. q] Even so, sQGP is still the most perfect liquid known. 

In summary: hydrodynamics+hadronic cascades reproduces all RHIC data on radial and elliptic 
flows of various secondaries, as a function of centrality,rapidity or energy are reproduced till pt ~ 2GeV, 
which is 99% of particles. Contrary to predictions of some, CuCu data match AuAu well, so Cu is large 
enough to be treated hydrodynamically. New round of studies last year included viscosity and relaxation 
time parameters on top of ideal hydro: viscosity values is limited to very small value. 

2.3 Jets quenching and correlations 

Pairs of partons can collide at small impact parameter: in pp collisions this produces a pair of large pt 
hadronic jets, which are (nearly) back-to-back in transverse plane (because total transverse momentum 
due to "intrinsic' parton pt is small). We can use therefore those high-pt partons as a kind of x-rays, 
penetrating through the medium on its way outward and in principle providing its "tomography" . 

We will not go into this subject in depth (see e.g. PHENIX "white paper" J26j but just note that 
accurate calibration of structure functions have been made in pp and dAu collisions, as well as with hard 
photon measurements (which are not interacting with the QGP). Thus we know quite well how many 
jets are being produced, for any impact parameter. The number of hard hadrons observed at transverse 
momentum pt relative to those expected to be produced as calculated from the parton model is called 
RAA{pt)- If this quantity is 1, it means the jets are all accounted for and none is lost in the medium. 
This is what indeed is observed with direct photons, not interacting with the matter, see Fig, 10 a). It 



was quite unexpected that for mesons this ratio RAA{Pt) was found to be rapidly decreasing and then in 
a wide range of momenta pt > AGeV its value is only RAA{pt) ~ -2 for central AuAu collisions, which 
means that 80% of jets are absorbed. 

^^Unfortunately I am skeptical about magnitude of systematic errors of any lattice results for rj/s (such as 42J): while 
the Euclidean correlation functions themselves are quite accurate, the spectral density is obtained by rather arbitrary 
choice between many excellent possible fits. 
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Figure 8: (From [21]) (Left) Centrality dependence of V2- The solid (dashed) hne results from a full 
ideal fluid dynamic approach (a hybrid model). For reference, a result from a hadronic cascade model is 
also shown (dash-dotted hne). (Right) Pseudorapidity dependence of V2- The solid (dashed) hne is the 
result from a full ideal hydrodynamic approach with T^^ = 100 MeV {T^^ = 169 MeV). Filled circles 
are the result from the hybrid model. All data are from the PHOBOS Collaboration. 
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Figure 9: (From [3H]) Elliptic flow V2 dependence on (Left) centrality and (Right) pt, compared to 
viscous hydro with variable viscosity. 



Theory of quenching mechanisms included gluon radiation on uncorrelated centers (with Landau- 
Pomeranchuk-Migdal effect) [43], synchrotron-like radiation on coherent fields |Hlll5], as well as losses 
due to elastic scattering. Comparing two radiation mechanisms in general let me just remind that 
the synchrotron-like radiation gives the energy loss dE/dt ~ E"^ growing quadratically with energy, 
it is stronger at high energies than radiation from uncorrelated kicks which gives only the first power 
dE/dt ~ E: but then correlation length in "GLASMA" fields and their lifetime is limited. As for elastic 
scattering losses, it depends on what are the couplings and especially masses of quasiparticles (quarks, 
gluons or maybe monopoles near Tc) on which scattering occurs. The rate of energy loss itself is an 
order of magnitude larger than pQCD predictions. Multiple phenomenological fits for the RAA{pt) were 
made: but they depend on models. 

Furthermore, as noted in my paper |16], any such model has predictions for ellipticity V2 in the range 
below its "geometric limit" for infinite quenching, while the data showed f 2 exceeding such limits for 
all models used. In simpler words, those models could fit RAA{pt) but not a double plot RAA{Px,Py)- 
The reason of such large f 2 is not yet found, to my knowledge. 



10: 



10" 



I 

o 



PHENIX Au+Au (central collisions); 

■ Direct y 

A 11^^ PreJiminary 

• 1 

GLV parton energy loss (dN'/dy ^1100) 



4a 



•I' 

Hi 



^mr^: 



S^f^f^^^^^'^'^'Ii^^ - 



AT 

I. 1 







S 

< o 





+ 



.8 
.6 
4 
.2 
1 
.8 

.6r 

4 



L I ' ' ' I ' ' ' I 



6 8 10 12 14 16 18 20 

Pt (GeV/c) 



I ' ' ' I ' ' ' I ' ' ' I 
STAR Prelimary 



.2r • 




^/'^'^^ 



t 



2 4 6 8 10 12 

p^ (GeV/c) 



§0.2- 



^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 



i — p+p min. bias 
r * Au+Au central 



(bi: 




J I I ! I I I I I I I I \ I L 



^2 



A(t) (radians) 



^ ^ Armesto et al. (1) 

I van Hees et al. 

,.. I 3l(2nT) Moore & 

1 12/(27iT) TeanGy (111) 




Au+Au @ \|sj^ = 2O0 GeV 

. . I . . . . I . . . . I :■ . . , ll 



7 8 9 

p^ [GeV/c] 



Figure 10: (left) Raa from PHENIX for 7, vr*^, rj and STAR for protons and charged pions. (right top) 
Two particle correlation function from STAR; (right bottom) Raa from PHENIX for single leptons 
(closed points) and tt" (open points) compared to theoretical calculations 



Further crucial test of this theory came from experimental observation of "single lepton" quenching 
and V2: those leptons come from semileptonic decays of c,b quarks. At the same pt heavy quarks 
have smaller velocity, and if the main quenching mechanism be radiative, it should reduce quenching 
accordingly. The data however do not show any serious reduction, with the same RAA{pt) ~ -2 value for 
single leptons as for pions (coming mostly from gluon jets). This fact cast doubts at any perturbative 
mechanism of energy loss, since re-scattering of a gluon should be larger than that of a quark by the 
Casimir (color charge) ratio 9/4. 

Moore and Teaney [17| developed a general framework of dealing with heavy quark dynamics in 
QGP, by invoking Focker-Plank or Langevin eqn. They have provided a general argument that if quark 
mass relative to temperature M/T is large, relaxation of heavy quark is happening slowly and thus 
justify the Langevin's uncorrelated kicks assumption. 



dpi 
dt 



^i{t) - r]DPi 



(e.WO(0) = '^M(^-0 



(28) 



Here rjn is a momentum drag coefficient and ^i{t) delivers random momentum kicks which are uncorre- 
lated in time. 3k is the mean squared momentum transfer per unit time. The usual diffusion constant 
D in space is related to those parameters by 
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Figure 11: (from Moore and Teaney) (a) The charm quark nuclear modification factor Raa and (b) 
eUiptic flow for representative values of the diffusion coefficient given in the legend. In this model, the 
drag is proportional to the velocity, -£ oc v. 



In FigjTTjwe show the calculated dependence of quenching Raa and elliptic flow V2 for leptons, resulting 
from the Langevin process (calculated on top of hydro evolution). As one can see stronger coupling 
leads to smaller Raa and larger V2: comparison with data (value about .2 for Raa and yellow band 
for ^2) clearly favor the smallest diffusion constant, about DpnT ~ 1. Further work on heavy quark 
diffusion by Rapp and collaborators [48] have tried to specify the diffusion constant from data better, 
and also suggested its explanation using heavy-light resonances. 

The next RHIC discovery was associated with "jet correlations" , which means that in events triggered 
by one hard particle with large pt one look for a second "companion particle" correlatecj_]with it. While 
in pp and peripheral collisions one sees "back-to-back jets" , with two peaks and relative azimuthal angle 
A0 values near zero or vr, nuclear collisions typically show no (or strongly reduced) peak at vr. Further 
subtraction of flow in the correlation functions revealed new peaks shifted from the direction of the 
companion jet by a large angle A0 ~ 1.2rad ~ 80°, see Fig, 12 'b). 

Where the energy of the quenched jets go? Thinking about this question at the time we came with 
the answer: provided energy is deposited locally, hydrodynamics should provide a detailed prediction. 
Thus new hydrodynamical phenomenon ^^ suggested in ^91150]. - the so called conical flow - is induced 
by jets quenched in sQGP. The kinematics is explained in Figjl2|^a) which show a plane transverse 
to the beam. Two oppositely moving jets originate from the hard collision point B. Due to strong 
quenching, the survival of the trigger jet biases it to be produced close to the surface and to move 
outward. This forces its companion to move inward through matter and to be maximally quenched. 
The energy deposition starts at point B, thus a spherical sound wave appears (the dashed circle in 
Fig, 12 eft ). Further energy deposition is along the jet line, and is propagating with a speed of light. 



till the leading parton is found at point A at the moment of the snapshot. The expected Mach cone 



-'^^It means that thousands of particles not correlated with the trigger are statistically subtracted. 

^^In the field of heavy ion collisions Mach cone emission was actively discussed in 1970's by Greiner et al: but it turned 
out not to work because nuclear matter - unlike sQGP - is not a particularly good liquid. The nucleon m.f.p. in nuclear 
matter is about 1.5 fm is not much smaller than the nuclear size. 
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angle is given by 



< c{sound) > 

cos{eM) = (30) 

c 



Here angular bracket means not only the ensemble average but also the time average over the time from 
appearance of the wave to its observation^ 

Experimental correlation functions include the usual elliptic flow and the serious experimental issue 
was whether the peaks I just described are not the artifact of elliptic flow subtraction. By Quark Matter 



05 this was shown not to be the case, see Fig, 12 [ top right), which shows PHENIX data selected in bins 
with speciflc angle between trigger jet and the reaction plane: the shape and position of the maximum 
(shown by blue lines) are the same while elliptic flow has a very different phase at all these bins. The 



position of the cone is independent on angle relative to reaction plane Fig, 12 right top), centrality (not 



shown here) and pt Fig, 12 right down): so one may think it is an universal property of the medium. 



The angle values themselves are a bit different, with 1.2 rad preferred by Phenix and 1.36 rad by Star 



data: those correspond to amazingly small velocity of the sound wave ^^ < Cg >= 0.36 — 0.2, indicating 
perhaps that what we see was related to the near — Tc region. Another evidence for that is observation 
of conical structure at low collision (SPS) energies, reported at QM08 by CERES collaboration: at such 
energies near-Tc region dominates. 

At the last Quark Matter 08 large set of 3-particle correlation data ( a hard trigger plus two 
companions) have been presented both by STAR and PHENIX collaborations. Although those are too 
technical to be shown here, the overall conclusion is that Mach cone structure is more likely explanation 
of the data than other possibihties such as "deflected jets". 

A number of authors have by now reproduced the very existence of conical flow in hydro, see e.g. 
Baeuchle et al[5Tj: but really quantitative study of its excitation is still to be done. Casalderrey and 
myself [52] have shown, using conservation of adiabatic invariants, that flreball expansion should in fact 
greatly enhance the sonic boom: the reason is similar to enhancement of a sea wave (such as tsunami) 
as it goes onshore. They also showed that data exclude 1-st order phase transition, because in this case 
conical flow would stop and split into two, which is not observed. 

Antinori and myself [53] suggested a decisive test by b-quark jets. Those can be tagged experimentally 
even when semi-relativistic: the Mach cone should then shrink, till it goes to zero at the critical velocity 
V = Cs = l/-\/3. (Gluon radiation behaves oppositely, expanding with decreasing v.) 

The experiments with tagged 6-jets seem to be even more important in view of recent studies by 
Guylassy et al [193[ 1192] who found (using AdS/CFT results from [157] and speciflc model relating it 
to Cooper-Fry formula) that when the velocity of the jet v approach Cg the angle of the peak does not 
accurately follow the Mach angle but remains always larger. They have further found that the main 
part of the peak comes not from cones but from the non-hydrodynamical near-jet zone (they call the 
"neck"): what is the nature of those large angle emission remains unknown. This group have further 
studied weak-coupling (pQCD-based) version of the near zone [192] . flnding that in this case most of 
the flow remains at small angles, with very small but visible peaks at Mach angle, but no trace of 
large-angle emission predicted by AdS/CFT. 

Let me flnally mention the main open question, which is the absolute and relative amplitudes of 
excitation of two hydrodynamical modes, the sound (responsible for the Mach cone structure) and the 
^^diffuson" mode, which show matter co-moving forward behind the jet. As emphasized in our works, 
this question cannot be answered from hydro itself, as close to the jet it looses its applicability. As we 
will see later, this ratio was recently found from AdS/CFT: but we don't know yet of it does or does 
not agree with the data. 

^^Recall that the speed of sound changes significantly during the evohition, becoming small near T^ 

^^Note, it is the velocity not squared. The speed of sound squared can be seen as a dash curve in Fig 15 right). 
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Figure 12: (left) A schematic picture of flow created by a jet going through the fireball. The trigger 
jet is going to the right from the origination point B. The companion quenched jet is moving to the 
left, heating the matter (in shadowed area) and producing a shock cone with a flow normal to it, at the 
Mach angle cosOm = Cg/v, where v, Cg are jet and sound velocities. 

(right top)The background subtracted correlation functions from PHENIX experiments, a distribution 
in azimuthal angle A0 between the trigger jet and associated particle. Unlike in pp and dAu collisions 
where the decay of the companion jet create a peak at A0 = tt (STAR plot), central AuAu collisions 
show a minimum at that angle and a maximum corresponding to the Mach angle 
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2.4 Charmonium suppression 

Charinonium suppression is one of the classic probes: since charm quark pairs originate during early 
hard processes, they go through all stages of the evolution of the system. A small fraction of such pairs 
~ O(10~^) produce bound cc states. By comparing the yield of these states in heavy ion collisions 
to that in pp collisions (where matter is absent) one can observe their survival probability, giving us 
important information about the properties of the medium. 

Many mechanisms of J /ip suppression in matter were proposed over the years. The first was sug- 
gested by myself in the original "QGP paper" [3], it is a gluonic analog to "photo-effect" gJ /ip — > cc. 
Perturbative calculations of its rate (see e.g. Kharzeev et al [54j) leads to a large excitation rates. 
Indeed, since charmonia are surrounded by many gluons in QGP, and nearly each has energy sufficient 
for excitation, one may think J /%p would have hard time surviving. That was the first preliminary 
conclusion: nearly all charmonium states at RHIC should be rapidly destroyed. If so, the observed J /ip 
may only come from recombined charm quarks at chemical freezout, as advocated e.g. by Andronic et 
al f55]. 

However the argument given above is valid only if QGP is a weakly coupled gas, so that charm quarks 
would fly away from each other as soon as enough energy is available. As was recently shown by Young 
and myself (SH], in strongly coupled QGP the fate charmonium is very different. Multiple momentum 
exchanges with matter will lead to rapid equilibration in momentum space , while equilibration in 
position space is very slow and diffusive in nature. Persistent attractions between c and c makes the 
possibihty of returning back to the ground state for the J/ip quite substantial, leading to a substantially 
higher survival probability. For the sake of argument, imagine the matter so dense that any diffusion 
of c and c is completely stopped: then, after this situation changes by hadronization, one would still 
find them close to each other and thus J/ip - with its by far the largest density at the origin - will be 
obtained again. Thus, strongly coupled - sticky - plasma may actually preserve the J/ip. 

Matsui and Satz j57j have proposed another idea and asked a different question: up to which T 
does charmonium survive as a bound statel They argued that because of the deconfinement and the 
Debye screening, the effective cc attraction in QGP is simply too small to hold them together. Satz and 
others in 1980's have used the free energy potential, obtained from the lattice, as an effective potential 
in Schreodinger eqn. 

F(r, r) ^ -^^ exp(-Mz)(T)r) + F{T, oo) (31) 

6r 

They have shown that as the Debye screening radius M^^ decreases with T and becomes smaller than 
the r.m.s. radii of corresponding states X;"^'; J/i^^ ^" ^ ^') ^■■■, those states should subsequently melt. 
Furthermore, it was found that for J /ip the melting point is nearly exactly Tc, making it a famous 
"QGP signal". 

Dedicated lattice studies [5H1 159] extracted quarkonia spectral densities using the so called maximal 
entropy method (MEM) to analyze the temporal (Euclidean time) correlators. Contrary to the above- 
mentioned predictions, the peaks corresponding to ?7c, J /ip states remains basically unchanged with T 
in this region, indicating the dissolution temperature is as high as T^ ~ (2.5 — 3)Tc. Mocsy et al [00] 
have used the Schrodinger equation for the Green function in order to find an effective potential which 
would describe best not only the lowest s-wave states, but the whole spectral density. Recently [61] 
they have argued that a near-threshold enhancement is hard to distinguish from a true bound state: 
according to these authors, the above mentioned MEM dissolution temperature is perhaps too high. 
My view is that the coUisional width of states in sQGP is probably large and thus the discussion of 
MEM data is rather academic: in any case we dont observe J /ip in plasma but after it, and a survival 
is a real-time issue which cannot be answered by the lattice anyway. 

Let us now briefly review the experimental situation. For a long time it was dominated by the SPS 
experiments NA38/50/60, who have observed both "normal" nuclear absorption and an "anomalous" 
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Figure 13: (From [56]). (a) Distribution over quark pair separation at fixed T = 1.5Tc after 9 fm/c of 
Langevin evolution, with (red squares) and without (green triangles) the cc potential. Strong enhance- 
ment at small distances due to potential is revealed. (b)The points are PHENIX data for anomalous 
suppression of J/ijj in AuAu min bias collisions Ji'^^i^"-'-"^^ (^y = q). Two curves are Langevin model, 
with (solid, upper) and without (dashed, lower) feed-down. 



suppression, maximal in central PbPb collisions. Since at RHIC QGP has a longer lifetime and reaches 
a higher energy density, straightforward extrapolations of the naive J/ip melting scenarios predicted 
near-total suppression. And yet, the first RHIC data apparently indicate a survival probability very 
similar to that observed at the SPS. 



One possible explanation [^21 E3] is that the J/i/j suppression is (nearly exactly) canceled by a re- 
combination process from unrelated (or non-diagonal) cc pairs floating in the medium. However this 
scenario needs quite accurate fine-tuning of two mechanisms. It also would require rapidity and momen- 
tum distributions of the J/tp at RHIC be completely different from those in a single hard production. 

Another logical possibility advocated by Karsch, Kharzeev and Satz [M] is that J/tp actually does 
survive both at SPS and RHIC: all the (so called anomalous, or nonnuclear) suppression observed is 
simply due to suppression of feed-down from higher charmonium states, ip' and x- (Those are feeding 
down about 40% of J/ip in pp collisions.) These authors however have not explained why J/ip survival 
probability can be close to one. 

Young and myself [56] did exactly that, followed Langevin dynamics of charm quark pairs, prop- 
agating on top of (hydro) expanding fireball. The treatment basically is the same as that discussed 
in the preceding section, where heavy quark diffusion constant has been derived. One new important 
element though is the cc effective potential, which we found is slowing down dissolution of the pair 
quite substantially, 1 see Figjl3|^a) leading to "quasiequilibrium" situation in which ratios of different 
charmonium states are close to equilibrium ones at corresponding T, while the probability is continue 
to leak into unbound pairs which occupy slowly growing volume. The main finding of this work is 
that the lifetime of sQGP is not sufficient to reach the equilibrium distribution of the pairs in space, 
allowing for a significant fraction of J/ip ~ 1/2 to survive through ~ 5 fm/c of the sQGP era. This 
probability for charmonium dissociation in sQGP is much large than in perturbative estimates, or for 
Langevin diffusion which would not include strong mutual interaction. We have not yet answered many 
other questions: e.g. what happens during ~ 5 fm/c of the "mixed phase". (In view of it seem to be a 
magnetic plasma, as we will argue below, the mutual attraction of charmed quarks gets only stronger 
there,) 
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3 From lattice QCD to sQGP 

This section has been more difficult to write than others, because a connection between lattice results 
and the "strongly coupled" regime of QGP at T = (1 — 2)Tc remains indirect. Perhaps by itself it is 
rather unconvincing for a critical reader, as it was for many lattice practitioners, who are still quite 
reluctant to accept the "paradigm shift" of 2004. There are quite serious reasons for that. One (which 
we will discuss in the AdS/CFT section) is that the difference between weakly and strongly coupled 
regimes is deceivingly small in thermodynamical quantities. The second reason is that by performing 
Euclidean rotation of the formalism and correlators, one indeed gets rid of the unwanted phase factors, 
but a heavy price for that is extremely limited ability to understand real-time transport properties - 
diffusion constants, viscosity and so on - which turned out to be at the heart of this debate. As usual, 
we start with introduction for pedestrians which experts should jump over. 

3.1 The QCD phase diagram for pedestrians 

QCD phase diagram is quite multidimensional: apart from the temperature T one can introduce chem- 
ical potentials ^j for each quark flavors u, d, s. One however only consider 2 combinations of those, the 
baryonic ^i, = (l/Nj) J2f A^/ and isospin fij = fiu — A^d chemical potentials. Then one can vary param- 
eters of the theory itself, such as the number of colors Nc or quark masses m/: in many cases however 
we will only discuss certain limits, for example in our discussion below a shorthand notation "2 flavors" 
Nf = 2 would imply massless u, d quarks and infinitely heavy (or just absent) s, c, b, t quarks. 

Three main phenomena will be under discussion: (i) confinement, (ii) chiral symmetry breaking and 
(iii) color superconductivity. The minimalistic phase diagram may have only three main phases: (a) 
hadronic, at low T and fib, which is both confined and has broken chiral symmetry; (b) Quark-Gluon 
plasma (QGP) at high T and /i/, where all kind of condensates are absent; and (c) color superconductor 
(CS) at high fib and low T. There can of course be many more phases, as these features are not really 
exclusive- e.g. there can be coexisting chiral and CS condensates. 

Note that two last phenomena are due to different kinds of pairing, chiral breaking (ii) is due to 
quark- ant iquark pairing, its nonzero order parameter is the "quark condensate" < qq >; while color 
superconductivity (iii) is due to quark-quark pairing and its nonzero order parameters are a set of 
diquark condensates < g"g^ >, with color indices a,b and the flavor ones i,j arranged in various ways. 
(The spinor indices will be always suppressed because it is believed that whatever color-flavor structure 
of the condensates can be, the diquark spin remains zero, in most cases except very exotic ones.) 

By confinement we mean more specifically electric confinement, which means that electric color 
field is expelled into flux tubes, making quark-antiquark potentials linear. Magnetic component of the 
gauge field is not expelled from the low temperature T < Tc phase: we will in fact see that it actually 



dominates il^ Linear potential is synonymous to the "Wilson area law" of a large Wilson loop C 

< Wc >~ exp[-aArea{C)] (32) 

with the nonzero string tension. From pioneering lattice calculation by Creutz in 1980's we know that 
pure gauge theories indeed have this feature, although mathematical "proof of that remains famously 
elusive. Polyakov introduced in 1978 his famous loop and argued that it provides a disorder parameter, 
with nonzero value at T > Tc. We will discuss issues related to this below in some detail. 

The area law and Polyakov loop average are of course true order parameters only in pure gauge 
theory without quarks. Since quark pair can be produced, the lattice static potentials show linear 
behavior up to certain distances. Many people thus think that in QCD with quarks the confining phase 

^°We will not discuss AdS/QCD in this review: let me still remind its practitioners that popular models of confinement 
- hard or soft "walls" ending the holographic space and forcing strings to generate linear potential - are over-simplistic, 
because the same thing happens with electric and magnetic strings. 
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Figure 14: (a) The Polyakov loop, with am = 0.01 and 16^ x 4 lattice, (b) Susceptibility of the chiral 
condensate, with am = 0.01 and two different lattices: 8^ x 4, 16^ x 4. The data are compatible with a 
crossover at the chiral transition. In the inset, the chiral condensate within the same range of f3 values. 
A clear jump is visible at the deconfinement phase transition, (c) The p parameter, with am, = 0.01, 
Lt = 4, for two different spatial volumes. 



cannot be strictly defined and thus there should be no real phase transition separating it from the 
deconfined QGP. 

This is however still disputed. Another possible order parameter (albeit nonlocal) is based on 
the idea of the "dual superconductor" by tHooft and Mandelstamm |103j which suggested long ago an 
existence of a nonzero magnetic condensate of some bosonic objects. Di Giacomo and collaborators have 
implemented the corresponding observable - called "Pisa order parameter" - which is an operator adding 
one explicit monopole to the vacuum. If so, the deconfinement should be a true phase transition: but 
numerical lattice data for generic nonzero quark masses show only a rapid crossover. This contradiction 
is not yet resolved. 

To show one example in which all used order parameters are studied together, let me discuss the 
work by Pisa group on adjoint QCD [HS] shown in Fig, 14 In this theory a quark has the same color 
as gauge particles, and - unlike the usual quarks - it makes chiral restoration temperature distinctly 
different ^^ from deconfinement. (Larger color representation (charge) of quarks is believed to lead to 
stronger qq pairing and thus higher melting temperature for the chiral condensate.) The Polyakov loop 
shown in figure (a) behaves as disorder parameter at deconfinement indeed, the same point as indicated 
by Pisa parameter (c), while the chiral condensate and its susceptibility (b) indicate much smoother 
higher-T chiral restoration transition. Thus in adjoint QCD we see a presence of one more phase (on 
top of the minimal list of three given above), a deconfined but chirally broken phase which is usually 
called "a plasma of constituent quarks". It has all the requisites of a chirally broken phase, such as 
massless Goldstone bosons - pions, etc. 

In general it is believed that with growing fi all transition become sharper. Arguments about likely 
shift in real-world QCD from near-crossover at zero /i and T^ ~ 170 MeV to real second-order point (the 
so called QCD critical point) and then first order line has been put forward by Rajagopal,Stephanov 
and myself [66], as well as few proposals how one can experimentally search for it. RHIC specialized 
run with greatly reduced beam energy is planned for 2009: perhaps it will shed light on this issue. 

What happens if one changes another famous parameter, the number of colors by increasing Nc"! 
Because gluons are adjoint and their effects are ~ N^, they dominate all quark-induced ~ N^ effects. 
It has been recently argued by McLerran and Pisarski [67] that this implies that at large enough 
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^^Most lattice works agree on coinciding deconfinement and chiral symmetry in fundamental 
it is still debated and I refrained from taking sides and show one but not the other plot. 
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"quarkionic" phase, thinking about quark-filled Fermi sphere but with baryonic excitations at the 
surface. Glozman [68J further suggested even more exotic possibility: particle excitations in form of 
(chirally symmetric) baryons, without baryonic holes. The issue is not yet settled and such phases 



have not been seen on the lattice ^ My view is that both these exotic ideas are perhaps excluded in a 
confined phase, so it should be just baryonic. 

In the discussion above we have ignored short 1-st order transition line between the "vacuum-like" 
and "nuclear matter-like" regions, at /x ~ Mn/3,T < Tc ~ 8MeV. Since the usual A^";, = 3 nuclear 
matter is Fermi-liquid, it is not qualitatively different from the bulk of hadronic phase: thus one can go 
around this phase transition. This is not the case for larger Nc, as in this case baryons are becoming 
heavier Mb ~ Nc while the nuclear forces are believed to have a smooth limit 

V{r, Nc) ^ y(r, Nc = oo) (33) 

The obvious consequence of this is crystallization, as kinetic energy gets subleading to the potential 



one. Examples of specific calculations with the skyrmions are well knowr^'^ [70J. Is there a possibility 
that solidification happens even for physical Nc = 3 QCD? Old Migdal's pion condensation was of this 
nature, and for more recent study of crystalline quark matter see Rapp, myself and Zahed [71]. There 
can also be a crystalline color superconductor, known also as LOFF phase, see more in the review [72] . 
Instead of going any further into the zoo of possibilities, let me stop with a joke: perhaps QCD would 
not have less phases than water does, and this is quite a lot. 

Finally, for completeness, let me mention the opposite direction, increasing the number of light 
quarks Nj. At some point asymptotic freedom will disappear, but before that there should be a critical 
line (or more) at which Banks-Zaks infrared fixed point will make the theory conformal in IR. As the 
corresponding fixed point coupling gets less and less, chiral symmetry breaking and confinement must 
go: again either together or separately. All this territory is amenable to lattice analysis but remains 
largely unexplored. 

3.2 Main QGP properties from the lattice 

The thermodynamical observables - pressure and energy density p(T),e{T) - from the lattice is the 
simplest global observable, thus they were calculated more than a decade ago (and used in hydro 
calculations). Let me show two recent plots from Karsch [75] which depict them in a combinations 
which reveal somewhat more than standard plots of p(T),e(T). 

The first combination e — 3p = T^^ is related to the famous scale anomaly: the fact that its value 
relative to p, e themselves goes to zero means that QGP gets more and more conformal. However the 
way it goes down is not 1/T^ - as simple bag model predicts: this phenomenon is not yet explained. 

As shown in [73] and earlier papers, quark and gluon quasiparticles dominate at T > 1.5Tc. But 
what happens below that? Gluon/quark masses are too high M ~ (3 — 4)T to explain the peak of e — 3p. 
Chernodub et al [73] have shown that lattice magnetic objects - monopoles and vortices - reproduce 
the shapq^ of this curve well. 

The second combination ^^ is p/e - now conformity is seen as the place where this ratio reaches 1/3. 



^^Chirally symmetric baryons can certainly exist: in fact Liao and myself [69j have argued that those are needed to 
explain lattice susceptibilities in the usual Nc — i QGP near deconfinement at T '^ /i. 

^•^These works also found another phase at higher density, in which skyrmions "fragment" into objects of fractional 
topology: they interpreted those as a "chirally restored" baryonic phase. I however don't know if one can trust the model 
all the way to this phase. 

^"^Unfortunately not the absolute normalization: this issue deserves further studies as the lattice used in this work is 
rather small. 

^^Its special role in hydrodynamics was emphasized j23j and its minimum got a spatial name - "the softest point" . 
Indeed, the gradient of pressure provides a force and energy density a mass to be moved, so p/e it is proportional to 
hydrodynamical acceleration of matter. 
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Figure 15: The trace anomaly, (e — 3p)/T'* (left) calculated on lattices with temporal extent A^^ = 6,8 
and the ratio of pressure and energy density as well as the velocity of sound obtained in calculations 
with the p4fat3 action on A^,- = 6 (short dashes) and 8 (long dashes). 

One can clearly see that it is not yet reached at RHIC, but it is the case at LHC. This is one of the 
reasons why LHC experiments will be decisive in proving (or disproving) whether the AdS/CFT duality 
can (or cannot) be used in the conformal window of finite-T QCD. 

Let me now turn to the screening lengths. As I already mentioned, QGP got its name after it 
was found p] that thermal gluons - unlike virtual ones - lead to electric screening of the charge in 
weakly coupled regime (high T). The corresponding electric (or a Debye) mass is M {electric) ~ gT. 
Static magnetic screening does not appear via perturbative diagrams; but it has been soon conjectured 
by Polyakov [73] that magnetic screening should appear non-perturbatively, at the smaller "magnetic 
scale" M{magnetic) ~ g^T. 

To illustrate whether lattice results on the screening masses are or are not in agreement with that, 
we show their T-dependence calculated by Nakamura et al [76], see Fig|l6|^a). Note that at high T the 
electric mass is indeed significantly larger than the magnetic one, but it vanishes at T^. - here electric 
objects gets too heavy and "electric part" of QGP effectively disappears. However magnetic screening 
mass grows continuously toward T^: thus the two cross each other, around T = l.bTf.. 

These observation were in fact the starting point for Liao and myself in thinking about "magnetic 
scenario" for the near- Tc region. We had used the screening masses to get an idea about density of 
electric and magnetic objects, one should conclude that QGP switched from electric to magnetic plasma 
somewhere around of 



T, 



E=M 



1.5T, -SOOMeF 



(34) 



The static potential between quark and anti-quark is another traditional observable, by means of 
which quark confinement in Non-Abelian gauge theories was established. It was originally inferred from 
heavy meson spectrum and Regge trajectories, and then studied in great detail numerically, through 
lattice gauge theories, for review see e.g. [78]. It is usually represented as a sum of a Coulomb part 
V ~ 1/r, dominant at small distances, and a linear part V = ar dominant at large distances. The 
latter, related with existence of electric flux tubes, is a manifestation of quark confinement. The string 
tension in the vacuum (T = 0) has been consistently determined by different methods to be about 



(T^ac. ~ (426Me1/)2 ^ 0.92 GeV/fm 



(35) 
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Figure 16: (upper left) Temperature dependence of electric and magnetic screening masses according 
to Nakamura et al [76]. The dotted line is fitted by the assumption, rUg ~ g'^T. For the electric mass, 
the dashed and solid lines represent LOP and HTL re-summation results, respectively, (upper right) 
Temperature dependence of the effective string tensions of the free and potential energies ap, ay. (down) 
The energy and entropy (as TSoo(T)) of two static quarks separated by large distance, in 2-flavor QCD 
according to [77] . 



Studies of the static QQ potential have been extended to finite T. In particular, deconfinement 
temperature T^ is defined as a disappearance of the linear behavior as a signal of deconfinement at 
T > Tc in the corresponding free energy F{T,r). Bielefeld-BNL group has published lattice results for 
static QQ free energy, as well as internal energy and entropy 



V{T, r) = F - TdF/dT = F + TS 



(36) 



at T both below and above Tc, see [THl ISO] . 

Remarkable features of these results include: 
1. The linear (in r) part of the potentials. Their effective tensions are shown in Fig. 16 top right). 



While that for free energy vanishes at T^ (by definition), that for potential energy extends till at least 

about 1.3Tc, with a peak values about 5 times (!) the (T^ac.- Similar behavior is seen in entropy, while 

canceling in free energy. The widths of these peaks provide a natural definition of "near-T^" region as 

T/Tc = 0.8 - 1.2 

2. Although potentials at large distances r — > oo are finite V{T, cxd), near Tc their values reach very large 

magnitudes, see FigJl6|(down). The corresponding large entropy S(Tc, oo) Ri 20 means that really huge 

~ exp{20) number of states is involved ; 

The origin of this large energy and entropy associated with static QQ pairs near Tc, remains mysterious: 

many attempts (e.g. [8T|) failed to explain it. Below we will return to this phenomenon in connection 

with "magnetic plasma" scenario. 

Before looking for explanations, however, let us focus on physical difference between F and U, based 
on papers by Zahed, Liao and myself [821 183] , in which they are related to what happens for slow and 
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fast motion of the charges. To be specific, let's consider a pair of static charges held by external hands. 
Suppose they are close initially Lj„j. —>■ 0, and then are separated to some finite L. This can be done in 
two possible ways, adiabatically slowly or very fast. The difference between them in thermodynamical 
and quantum-mechanical contexts are known in many fields of physics. Perhaps the oldest is the so 
called Landau-Zener problem [81] of electron motion, following the motion of two nucleus in a diatomic 
molecule. While nuclei change their relative distance L with velocity f 12 = dL/dt, the electrons are in 
a specific quantum state with the energy depending in L. The issue is probability of the level crossings, 
which appear when there are two quantum levels Ei{L) ^ criL + Ci E2{L) ^ a2L + C2 crossing each 
other, at some separation Lq. When the two nucleus approach the crossing at Lq very slowly, then the 
electrons may jump from one energy state to another, always selecting the lowest energy state. However 
if the two nucleus move fast, there is large probability for the electron to remain in the original state. 
More quantitatively, Landau-Zener showed that this probability is given by 



Pfast = exp 



27l\Hi2 



2 n 



t;i2kTl -(72 



(37) 



where H12 is the non-diagonal transition matrix element of the Hamiltonian; but we will only need the 
limits of large and small f 12. The adiabatic limit obviously corresponds to free energy F{T, L) measured 
on the lattice. The "potential energy" V{T, L) means that no entropy is generated: this implies that 
there was no transition from the original pure state at T=0 into multiple states as level crossing occur: 
thus it corresponds to fast motion limit. The positivity of entropy means that V > F always. 

This discussion is very relevant for the problem of effective potential to be used in Schreodinger eqn 
for the bound states, e.g. in charmonium problem. Zahed and myself [H2] argued that in this case one 
should use the internal energy: provide much more stable bound states, delaying J/ip melting to higher 
T ~ 3Tc. Several authors (e.g. [85]) have used effective potentials in between those two limiting cases. 
With such potentials not only charmonium but also light quark mesons get bound, as also are baryons 
and "gluonic chains" [86] and also colored binary states. However, in a liquid with the parameters we 
expect from such interaction the number of nearest neighbors associated with one charge is expected 
to be ~ 4, and thus it is not clear what is the role of the binary states. Quantum manybody studies of 
these systems are not attempted yet, and we don't know if there is any sense to identify them, and if 
so how wide those states may be. 

3.3 Polyakov loop, "Higgsing" and deformations of QCD 

Physics of monopoles, to which we will turn shortly below, has been originally developed in the Georgi- 
Glashow model or Af=2 SUSY which has adjoint scalars. Those may have some nonzero expectation 
values - this phenomenon would be colloquially referred to as "Higgsing" below. If so, the color group is 
broken, generically to diagonal U{1)^ . QCD-like theories are much more difficult precisely because they 
do not have elementary scalars, making Higgsing much more subtle, with its role at finite T presumably 
played by the zeroth component of the gauge field Aq. 

The definition of the Polyakov loop [75] is a holonomy of the gauge field across the periodic direction 



/•/3 

U = Pexp{i / Aodr) (3J 

Jo 



where /3 = 1/T is the Matsubara time. If it has VEV one can think of < f/ > as a diagonal color 
matrix, with some eigenvalues vi...vn'- Polyakov's view on confinement is that the eigenvalues widely 
fluctuate and < Tr{U) >= 0: in the deconfined (plasma phase) these eigenvalues fiuctuate little near 
minima of the effective action, at Vk = e^'^*'^/^. At high T - that is in weak coupling - one can calculate 
this effective potential perturbatively. In zeroth order, Aq = const costs no energy, but its coupling 
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Figure 17: Schematic phase diagram for adjoint QCD with Shifman-Unsal deformations, on a plane 
deformation parameter a - temperature T. At zero deformation a = (left side of the plot) one has the 
usual adjoint QCD in which deconfinement happens below chiral symmetry restoration. However, as a 
grows and we deform into a theory without deconfinement, the order must be reversed, giving place to 
a new phase in the right upper corner, which is confined but chirally symmetric. 



to gluons in the heat bath leads to shifting the gluon statei 
effective action 
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and as a result one gets the following 



(39) 



This pushes TrlJ away from zero, to the minima mentioned. One set of lattice data on < P(T) > has 
been already shown in Fig, 14|^a): as one can see the breaking indeed happens and thus it is indeed a 
disorder parameter. 

Opinions on the role of the breaking of the Zfq symmetry vastly duffer: while some think it the very 
essence of confinement others think it has no dynamical meaning at all. To exemplify this polemics 
the reader may e.g. see Ref. [HZ| in which Smilga pointed out that it is highly suspicious that the 
pure gauge theory - which has no such symmetry and can be formulated as 50(3) gauge theory rather 
than SU{2), with Z2 explicitly eliminated. Smilga further argued using simpler examples that effective 
potential is gauge-dependent concept and should be treated with care, he emphasized that although at 
finite lattices/coupling one may apparently see Zj^ domains, only one minimum is physical and the so 
called Z{N) bubbles with nonzero surface tension are not really there. For comparison of these two 
lattice formulation see deForcrand [88]. The question is no longer debatable when there are fundamental 
fermions in the theory, as they see different Z^ phases and thus no longer respect the symmetry and 
one vacuum with real Polyakov line VEV is preferred. 

Before we discuss phenomenology of the Polyakov line and its effective potential in strong coupling 
(on the lattice) in real-world QCD, let us turn to interesting "deformations" of QCD recently dis- 
cussed by Unsal [89] and Shifman and Unsal [90] . The central idea is to deform some QCD-like theory 
continuously into something else to which the answer is either known or can be obtained perturbatively. 
Although they consider several different fermion representations, for simplicity we will only consider 
in this section the case of adjoint fermions, for which the interplay of gluons and fermions is simpler. 
Before going into details, here is the list of deformations: 
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Here we mean gluodynamics only: we turn to quarks a bit later. 
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(i) rotating fermionic boundary conditions along the time from anti-periodic to periodic 
(ii) putting the theory into a spatial box with variable size 
(iii) introducing additional potential for the Polyakov loop. 

The first deformation is obtained via the introduction of the phase into fermion boundary condition, 
allowing to interpolate smoothly between the periodic and anti-periodic boundary condition 

(— ) => exp{ia + m) (40) 

into fermionic contribution into a line, with a changing from (anti-periodic) to vr (periodic) continu- 
ously. The effective potential becomes then 

Vefs-Y.(Nfe^^-^-)--l)\-—^ (41) 

n "■ 

When a = vr {periodic boundary conditions) and the number of flavors Nf = 1 two terms in the 
effective action simply cancel: this happens because this theory is the Af=l SUSY gluodynamics with 



the supersymmetry remaining unbroker ^^ The famous argument based on Witten index applies in this 
case, telling that the number of vacua cannot be changed with any deformation. When Nf > 1 the 
fermionic terms dominate and the sign of the potential is reversed. It means in this case one has the 
theory in which U is not pushed to large values and there is no Zn breaking and thus no deconfinement 
even at weak coupling. 

Let us now think about deconfinement and chiral symmetry of the adjoint QCD for 2 < Nf < 5 (the 



upper limit from asymptotic freedom) in the a — T plane, shown in Fig, 17 When a = we have the 



usual (undeformed) adjoint QCD in which (see Fig, 14) T^ > Tconf- But when a = ir there should be 



no deconfinement phase at all, which means that Tconf (en) grows indefinitely before crossing the a = ir 
vertical line. There should however still be chiral symmetry breaking at any a: in fact the a = it 
theory with periodic fermions has well known dyons with 2Nf fermionic zero modes, which generate 
NJL-type interaction and chiral symmetry breaking provided T is small enough to get the coupling 



large enough ^^ As a result, there should be an intersection between the deconfinement line Tconf{c() 
and the chiral restoration line T^{a) at some acrossing, as shown in Fig, 17: after the deconfinement and 



chiral symmetry restoration lines have crossed one finds a qualitatively new phase, with confinement 
but with unbroken chiral symmetry . 

Now let us turn to the second Shifman- Unsal deformation of the QCD-like theories: if one formulates 
the theory on R^ x 5*1 space with (periodic for fermions) compact direction of the variable length L, one 
can gradually interpolate between 4-d and 3-d gauge theories. The major difference between those, as 
explained by Polyakov [93] many years ago, is that while 4-d instanton-antiinstanton interaction is short 
range 1/r^, the 3d instantons (that is, monopoles) interact by a long range magnetic Coulomb 1/r. The 
result is that the 3d theory is confined by monopoles-instantons, while the 4d theory is not confined by 
its instantons. Moreover, it happens even in the weak coupling regime, in which the instanton-monopole 
density is exponentially small. 

Specific mechanism for QCD on R^ x S^ was discussed by Unsal in a separate paper [91] for pe- 
riodic fermions, it is condensation of magnetic charge 2 "bions" - pairs of certain dyons - bound by 
fermion-induced forces in spite of mutual magnetic repulsion. The binding is analogous to instanton- 
antiinstanton molecule formation, the confinement at high T is like Polyakov mechanism in 3dim ^93]. 

The third Shifman- Unsal deformation is done by an addition to the QCD action of artificial potential 
for the U, e.g. 

Kdd~E«-l^^^"l (42) 



^^We remind the reader that normal thermal boundary conditions obviously break the symmetry between fermions and 
bosons. 

^^The physics of chiral transition is the same transition between "instanton liquid" and "instanton-antiinstanton 
molecules described by Ilgenfritz and myself [91j, see also references for later studies in the review [92j . 
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with some coefficients a„ chosen at will. The authors themselves argued that if the coefficients are 
chosen in order to oppose the VefjiU) generated by quantum fluctuations naturally, one should be able 
to delay deconfinement (increase Tconf), to the extent that it will occur in the weakly coupled domain 
and make it tractable in a (semiclassical) controlled approximation. 

Perhaps it is also interesting to go to another direction as well, decreasing Tconf, reaching for the 
regime in which electric theory is even stronger coupled than usual but its dual -magnetic theory of 
monopoles - will gets perturbative instead. Both deformations can easily be done on the lattice: a 
possibility to check continuity of the underlying physics of both deconfinement and chiral restoration 
all the way from strong to weak coupling will surely contribute a lot into our understand of both. 

3.4 Phenomenology of Veff{U) in QCD 

The early history of perturbative derivation of the effective potential for U (or Ao = const) can be 
found in the classic review [95]. Instead of repeating here well known results, let me refer to more recent 
attempts to combine known perturbative results with the lattice data include a paper by Pisarski P5] 
where one can find the details of the recommended effective Lagrangian. 

One more motivation to study the QCD deformation via adding extra potential for Polyakov loop on 
the lattice comes from heavy ion phenomenology. Dumitru and collaborators [97] have used this form of 
the effective Lagrangian to study real-time evolution of Aq. The main conclusion from their work is that 
that < Aq > belongs to the class of so called slow variables, and its evolution in heavy ion collisions has 
to be treated separately from the overall equilibration. They have performed numerical solution of the 
EOM for it, starting from "suddenly quenched" value corresponding to its vacuum form, moving toward 
its minimum at the deconfined phase at T = 2Tc. The main finding of this work is that the relaxation of 
this variable is very slow, taking about 40 fm/c or so. This time significantly exceeds the QGP lifetime 
at RHIC which is only about 5 fm/c or so, which suggests that in real collisions we should treat Aq as 
essentially random variable frozen at some value and color direction during hydro evolution. It means 
there is a chaotic out-of-equilibrium Higgsing, slowly rolling down, like in cosmological inflationary 
models: thus one would like to know as much as possible about phase transitions and EoS for all values 
oi<U>. 

Another active direction using effective potential Veff{U) is the so called PNJL model [98], which 
combines the Polyakov loop with well known Nambu-Iona-Lasinio model for chiral symmetry breaking, 
see also [99]. Quite impressive results for QCD thermodynamics were obtained along this path by 
the group of Weise jlOOj . Let me give their notations and the parameterization of the potential. A 
background color gauge field = ^4 = iAQ, where Aq = 6^q gA^ t"- with the SU{?>)c gauge fields A^ and 
the generators t" = A"/2. The matrix valued, constant field relates to the (traced) Polyakov loop as 
follows: 



^ = w7^ 



1 c 



V exp 




^Tre*<^/^ , (43) 



In the so-called Polyakov gauge one chooses a diagonal representation for the matrix (j), (j) = (f)^ As+^g Ag , 
which leaves only two independent variables, 03 and 08- The potential involves the logarithm of </($), 
the Jacobi determinant which results from integrating out six non-diagonal SU{'i) generators while 
keeping the two diagonal ones, 03,8, to represent $: 



y4 
with 



U{^,^\T) 1 

- = --a[ 



(T) $*$ + h{T) In [1 - 6 $*$ + 4 (<l>*^ + ^^) - 3 ($*$)^] (44) 



a{T) = ao + ai ( ^) + a2 (^) ' , 6(T) = h(^)\ (45) 
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Figure 18: Comparison of PNJL model predictions |lUU] with susceptibilities of different order, from 2 
to 8 derivatives over yU at /i = with the lattice data |101j . 



The logarithmic divergence of W($, $*,T) as $, $* ^ 1 automatically limits the Polyakov loop $ to be 
always smaller than 1, reaching this value asymptotically only as T — i> oo. The parameters Oj and 63 are 
determined to reproduce lattice data for the thermodynamics of pure gauge lattice QCD up to about 
twice the critical temperature^^ The values of these parameters are Oq = 3.51 , Oi = —2.47 
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15.22 , 63 = —1.75. The critical temperature Tq for deconfinement in the pure gauge sector is fixed at 
270 MeV in agreement with lattice results. 

After performing a bosonization of the PNJL action and introducing scalar and pseudoscalar auxil- 
iary fields, a and vf, the PNJL thermodynamic potential becomes: 



O" 



fi(T,/x„a,<l>,$*) = W(<l>,$*,T) + — -2iV^ 



/ ^ P (t In [1 + 3 $ e~'^^''~^''^^^+ 3 $* Q-^iEp-f^,)/T _^^-3{Ep-f,,)/T 
J r27rl^ ^ L 



+ Tln 



1+ 3 $* Q-iEp+f^g)/Tj^ 3 ^ ^-2{Ep+f,g)/Tj^ ^-3(Ep+f,g)/T 



+ 3AEpeiA' 



p')} 



(46) 



where the quark quasiparticle energy is Ep = y/p'^ + m'^ and the dynamical (constituent) quark mass 
is the same as in the standard NJL model: m = rriQ — a = itlq — G{ipip). The last term in the previous 
equation involves the difference AEp between the quasiparticle energy Ep and the energy of free fermions 
(quarks). 

In order to see how this model works, the reader can have a look at FigJTs] the agreement is very 
good. The key feature of the expression above is the suppression by the Polyakov loop, which appears 

^^At much higher temperatures, where transverse gluons begin to dominate, the PNJL model is not supposed to be 
applicable. 
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as $ in the single quark (or antiquark) term in square brackets, as $* for diquarks and finally there is 
no $ in the last term corresponding to 3 quarks in colorless combination. The effect of all of that is that 
quark and diquark terms are suppressed in the deconfined phase, while 3-quarks are not: this drastically 
improves the prediction of the original NJL model. The reader should note, that unsuppressed colorless 
3-quark term is basically nothing else but a baryon - thus the main lesson of this work is actually quite 
similar to that of Liao and myself [69], namely that the QCD thermodynamics at finite /i, both below 
and right above T^, is still dominated by haryons rather than by individual quarks. 

4 Electric-magnetic duality and finite-T gauge theories 

Electric/magnetic duality is repeatedly resurfacing during the history of physics, serving as a source of 
inspiration for a while and then going dormant again. I would put the first occurrence of that to famous 
Maxwell's guess about displacement current, guided by a "dual relation" containing time derivative of 
magnetic flux discovered by Faraday. From then on, vacuum Maxwell eqns are of course perfectly E-B 
dual. 

In the first years of quantum mechanics Dirac famously shown that the wave function can only 
be consistently defined if the electric and magnetic coupling constants are related by the celebrated 
quantization condition 

a{electric)a{magnetic) = n (47) 

where n is some integer. With advent of quantum field theory, renormalization and running couplings, 
this condition elevates into a requirement that these two couplings must run in the opposite directions: 

P{electric) + ^{magnetic) = (48) 

Thus when a{electric) = e'^/An is small (at high T), a{magnetic) = g'^/Air should be strong. We will 
discuss implications of that for QGP below. 

In 1974 't Hooft and Polyakov have discovered monopole solutions in gauge theories. The specific 
setting was Georgi-Glashow model, which is a gauge theory supplemented by adjoint scalar providing 
"Higgsing". Recall that when VEV is diagonalized, all Nc — 1 diagonal color matrices commute with it. 



which leave those U{1) "photons" massless^" (There is huge literature about monopoles in this setting 
and in supersymmetric theories, for non-experts the book by Shnir [102j can be recommended.) 

G.'t Hooft and Mandelstamm [103] famously suggested to explain confinement by a "dual supercon- 
ductor" made of Bose-condensed magnetically charged objects. Seiberg and Witten |1U4] have famously 
shown how it works in the N=2 super Yang Mills theory. They have shown spectacular example of "tri- 
ality" between electric, magnetic and "dyonic" descriptions near the corresponding singularities on the 
moduli space. Work continues on the supersymmetric front, populating our "topological zoo" by new 
species, for good pedagogical review see e.g. |1U5] . Let me just mention two recent additions: the non- 
Abelian magnetic strings and monopoles which live inside them as "beads" |1U6] or duality between 
gauge theory on domain walls and in the bulk |107] . All these developments are of course the pillars of 
our understanding of magnetic objects: which however remains quite incomplete for the very theory it 
is needed, the QCD. One obvious way to find and study dynamical role of magnetic objects in various 
observables is the detailed analysis of lattice configurations. Such approach was very successful for 
instantons (for review see e.g.pSJ)- Let me comment on those works related to monopoles, dyons and 
fermionic dyons, consequently. 

Lattice observation of the monopoles require fixing the gauge, as explained long ago e.g. in DeGrand 
|108] , and then calculating the total magnetic flux through all elementary 3-d boxes. Naively, since each 
link appears twice there should be always zero flux: but since the phase is counted modulo 27r there 



^°In SU(3) theory those are 2 - A3, Ag - massless and 6 massive gluons. 
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are nonzero fluxes in some boxes. It literally means that (in a particular gauge) we look where singular 
Dirac strings end. The monopole current can then be defined: monopole boxes form closed loops. 
Are those objects physical or just some UV noise in unphysical gauge choice? The answer to this 
long standing question in general is still not completely clear: however better posed specific questions 
recently produced so reasonable and consistent answers - apparently independent of the particular 
lattice parameters - that they should be physical (in my opinion). 

In particular, if one counts only the loops which are winding around the lattice at least once, then 
their density was shown to scale correctly with lattice parameters. As we show in the next section, the 
correlation function between those monopoles behave "physically", in a way consistent with what one 
expects from a running magnetic coupling. In vacuum, they run around confining electric flux tubes, as 
"magnetic coil" should do. Above T^ they show liquid-like correlation in quantitative agreement with 
MD for classical electric/magnetic plasmas. All of this prompts me optimistically assume that those 
objects are physical. (Yet we stiH don't quite understand the objects themselves, don't know their 
masses and structure, and we even cannot explain why in some other gauge one cannot see them at 
aU.) 

From the point of view of the deflnition, another kind of magnetic objects - the self-dual dyons - are 
much easier. Kraan-van Baal solution [lllj has shown that already at the classical level, in the presence 
of nonzero holonomy or "Higgsing" , flnite T instantons are decomposed into Nc self-dual dyons. 

Because these dyons are descendant of instantons, they may be seen "via the eyes of fermions" : 
and indeed such methods were developed by Gattringer, Ilgenfritz and others, with stunningly beautiful 
examples of such dyons in gauge conflgurations. Diakonov and collaborators have evaluated quantum 
one-loop correction to Kraan-van Baal solution [112j and discussed the ensemble of dyons [113j . Yet 
there are serious reasons to think that dyons are less important than monopoles: . (i) Their density is 
not peaking near Tc; (ii) They seem to be more massive than monopoles (we will return to this issue 
in the summary of the chapter); and - last but not least - (iii) Bose-condensed dyons would not create 
only electric conflnement, as needed, because they are selfdual and thus cannot tell E from B field. 

Finally, massless fermions have zero modes on topological objects, including the monopoles/dyons. 
Supersymmetric theories - e.g. the celebrated Seiberg-Witten M=2 theory - must have it as a re- 
quirement, since all monopoles/dyons should have fermionic superpartners. I would think that in QCD 
quarks should also be able to "ride on a monopole": but their masses, density and dynamical role is 
not yet investigated. 

4.1 A "magnetic scenario" for near-T^ region 

In the chapter related with applications of hydrodynamics it was called traditionally "the mixed phase" , 
as the model EoS used still had a phase order transition and Maxwell construction at intermediate (en- 
tropy) density. This was done along the traditional thinking, although we know for fact that none of 
the overheating-overcooling fiuctuation phenomena (known near any first order transition) are actually 
observed in heavy ion collisions. The so called event-by-event fiuctuations there turn out to be remark- 
ably small and simple. So, this is definitely not a QGP and hadronic phase mixed together in some 
kind of emulsion. 

The question is, what is the adequate picture for the near- Tc region? Liao and myself |114j proposed 
a new and very simple view on this question - "the magnetic scenario" - which basically suggests it to 
be a magnetic plasma of monopoles, in a liquid form. Another line of work based on lattice data lead 
Chernodub and Zakharov |115j to the same conclusion. 

Let me start with a qualitative discussion of the running couplings on the phase diagram, which is 
a crucial element of the argument. As QGP produced is cooling down and T decreases toward T^ the 
electric coupling grows and the magnetic one decreases, in accordance with asymptotic freedom and 



Dirac condition (47). The electric objects - quarks - gets heavier and more strongly coupled, while 
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the monopoles get lighter and less correlated: thus one should expect an equilibrium at some T, with 
dual language based on monopoles being simpler. A schematic phase diagram explaining these ideas 



is shown in Fig, 19 top left). The main assumption is that plasma becomes magnetic-dominated at T 
above the deconfinement line, interpreted as a BEC transition of magnetic media. 

In the beginning of this chapter we already defined the near-T^ region as T/T^ = 0.8 — 1.2, based 
on the width of various "specific heat" peaks in static dipole energy. Recent lattice data |109j provided 
dramatic conformation of this scenaric^^ Figjl9|(b) shows two sets of these data, and the correlation 



(and thus magnetic coupling) is indeed stronger at higher T. Furthermore, the correlation function for 
50-50 mix of electric/magnetic plasma obtained in our Molecular Dynamics (MD) simulation Figjl9|c) 
has the same shape and magnitude, provided one compare at the same value of the magnetic plasma 
parameter F = a{'magnetic) / {^^Y^^ /T: its extracted values are shown in Fig, 19|^d). It is very nice to 
find always F > 1, which means that magnetic component of sQGP is also liquid not gas, thus it does 
not spoil the "perfect liquid" at RHIC One may further think that viscosity has a minimum where both 
electric quasiparticles (quarks) and magnetically ones (monopoles) have similar difficulty propagating. 
We infer from lattice data that such electric-magnetic equilibrium is at T ^ l-5Tc, right in middle of 
the RHIC domain. 

We already reviewed the main properties of lattice potentials in the previous chapter. A good starting 
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Even more recent lattice work by Meyer |110j have found liquid-like peaks in the correlators of the energy as well. 
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point for development of such a model is to think along the line of "dual superconductor" by t'Hooft- 
Mandelstamm [lOSj . If the QCD vacuum at T=0 contains some magnetically charged condensate, it 
will expel electric flux between QQ into a flux tube, analogous to the Abrikosov-Nielsen-Olesen (ANO) 
vertex. 

An explanation of the linear rise of energy and entropy to very large values is due to persistence 
of a flux tubes even in the plasma deconfined phase. This point of view was proposed in |114] and was 
further developed in [116j , for infinitely long flux tubes. In the latest paper |117j we extended the model 
to finite-distance potentials, by developing an analytic "elliptic fiux-bag" model, with two static quarks 
at its focal points. Omitting all technical details, we will only present two main results for the tensions 
of "slow" and "fast" potentials. In the former case we obtained 

^ = 1.69 X ai^ X V]^^ (49) 

and the saturated value of transverse radius to be 

7^5 = 0.82 X aj,^ X Vll'^ (50) 

For "fast" potential, related to "normal" component of magnetic plasma, the corresponding results are 

^^ = 2.94 X a]^^ X n]^^ (51) 

and the saturated value of transverse radius to be 

TZf = 0.57 X a]j^ x n^^^ (52) 

Using these expressions and lattice potentials/tensions one can estimate the monopole density needed 



to produce the energy string tension at T^, (Te{Tc) ~ h.lha^ac- From (51) we obtain the monopole 
density should be nMiT^ ^ 4.65/m~^ assuming aE = 1, or UMiTc) ~ 6.58/m~^ assuming ue = 0.5. 
This again is a reasonable number as compared to direct lattice observations . The transverse size is 



then 7?.^ ~ 0.34/m for ue = 1 and TZf ~ 0.24/m for ue = 0.5. The results are shown in Fig, 20 for 
Bose condensed component (upper) and "normal" one (lower), the latter compared to lattice data on 
directly observed monopole paths (winding ones around time axis). We will return to these numbers at 
the summary of this chapter below. 

4.2 Molecular dynamics for magnetic/electric plasmas 

Gelman, Zahed and myself [5] proposed a classical model for the description of strongly interacting 
colored quasiparticles as a nonrelativistic Non-Abelian Coulomb gas. The sign and strength of the 
inter-particle interactions are fixed by the scalar product of their classical color vectors subject to 
Wong's equations. Details should be looked at the papers: let us just explain here its physical meaning. 
For SU(2) color group a color vector rotates around the direction of the total color field induced by all 
other particles at its position: same as magnetic moments would do in a magnetic field. For arbitrary 
group precession on a group is determined by the Poisson brackets of color vectors: {Q"', Q'^} = /"'"^ Q'^ 
which are classical analogue of the SU(Nc) color commutators. Thus for arbitrary gauge group one 
should use its structure constants /"'"^ describing "precession" of the color vector. 

For the non-Abelian group SU(2) the adjoint color vectors resign on a 3-d unit sphere: with one 
conserved quantity {{Q"')^) it makes S^ or 2 degrees of freedom. For SU(3) the group has 8 dimensions: 
with two conserved combinations {Q'^Y^d'^^^Q'^Q^Q'^ it is 6 d.o.f., and so on. Although color preces- 
sion equations do not look like the usual canonical EoM for pairs of coordinates and momenta, they 
actually can be rewritten as pairs of conjugated variables, as can be shown via the so called Darboux 
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Figure 20: (top left) free energy F (in unit of aE/lc) versus separation L/lc] (top right) monopole 
condensate energy density £^^/^ in unit of ^o^ac-, the two curves are for oie being 0.5(upper) and l(lower) 
respectively. 

(down left) The potential energy V (in unit of aEJlu) versus separation L/Im] (down right) ther- 
mal monopole density hm/T^, the two curves across boxes are for aE being 0.5(upper) and l(lower) 
respectively, and green curve across diamonds shows data for T > 1.3Tc from D'Elia and D'Alessandro. 



parameterization. Thus one can even define the phase space and use all pertinent theorems related to 
its classical evolution, if needed. 

The model can be studied using Molecular Dynamics (MD), which means solving numerically EoM 
for 10^ — 10^ particles. It displays strong correlations and various phases, as the Coulomb coupling is 
increased ranging from a gas, to a liquid, to a crystal with anti-ferromagnetic-like color ordering. There 

one can 



is no place for details here: so we simply jump to results on transport properties. In Fig, 21 



see the result for diffusion and viscosity vs coupling: note how different and nontrivial they are. When 
extrapolated to the sQGP suggest that the phase is liquid-like, with a diffusion constant D ^ 0.1/T 
and a bulk viscosity to entropy density ratio rj/s ~ 1/3. 

Transport properties for novel types of plasmas, including electric and magnetic charges, have been 
calculated by Liao and myself |114] : and t] is indeed minimal for most symmetric mixture 50-50%. 
Before we turn to these results, let me qualitatively explain why in this case the diffusion/viscosity is 
maximally reduced. Imagine one of the particles - e.g. a quark. The Lorentz force makes it rotate 
around a magnetic field line, which brings it toward one of the nearest monopoles. Bouncing from it, 
quark will go along the line to an antimonopole, and then bounce back again: like electrons/ions do in 
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Figure 21: The diffusion constant (a) and shear viscosity (b) as a function of the dimensionless couphng 
r. Three sets of data are for only electric particles (MOO), with the quoter (M25) and the half (M50) 
particles being monopoles. Note the decrease in viscosity as the admixture of monopoles grows. 



the so called "magnetic bottle" ^ Thus in 50-50 mixture all particles can be trapped between their 
dual neighbors, so that the medium can only expand/flow collectively. 
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Figure 22: Transport summary from |114] : Log\\- j {ji j s')\ v.s. Log\V j {2tiT D^ including results from our 
MD simulations, the AdS/CFT calculations, the weakly coupled CFT calculations, as compared with 
experimental values, see text. 



Our MD results are shown on viscosity-diffusion plane in Fig, 22 by three lines: they are compared to 
those from the AdS/CFT correspondence in weak and strong coupling as well as with empirical values 
from RHIC experiments ( gray oval). The dashed curve in the left lower corner is for M=4: SUSY YM 
theory in weak coupling: , where viscosity for this theory in weak coupling is from |118] and diffusion 
constant from[119j. The curve has a slope of one on this plot, as in weak coupling both quantities are 
proportional to the same mean free path. These weak coupling results are quite far from empirical 
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data from RHIC in the right upper corner. (Viscosity estimates follow from deviations of the elliptic 
flow at large pt from hydro predictions and diffusion constants are estimated from Raa and elliptic 
flow of charm .) The strong-coupling AdS/CFT results (viscosity according to |120] with 0(A~^/^) 
correction, diffusion constant from |121] ) are represented by the upper dashed line, going right through 
the empirical region. Our MD results - three solid lines on the right - are close to the experiment as 
well, especially the version with the equal mixture of electric and magnetic particles. 

4.3 Bose condensation: from liquid He to monopole plasma 

Classical approximation discussed above can explain many properties of liquids and solids, but it ob- 



viously ignores quantum effects '^^ It is well known that quantum effect at low temperatures may lead 
to qualitative changes in the system's behavior, such as superfluidity and superconductivity. A system 
of bosons may undergo Bose-Einstein condensationm (BEC). While for ideal gas BEC is a textbook 
material , it is a difficult problem for interacting systems. Liquid He^ remained for a long time the only 
example and the relation between its superfluidity and BEC phenomenon was hotly debated and even 
denied in many classic works. The dilute atomic gases were flnally cooled to BEC in 1990 's, but those 
are weakly coupled and we would not discuss them. 

Liquid He^ problem became amenable to direct numerical attack in 1980's, when simulation of the 
Feynman path integral via Monte Carlo algorithms was flnally technically possible. Such flrst-principle 
approach was a success: see e.g. Ceperley's review |125] . After that, other systems were studied, flUing 
the gap between weakly coupled gases and liquid He"^. And yet when we studied this literature - trying 
to understand conditions for BEC of the monopoles -it still looked that the net was not cast wide 
enough and many general questions remained open. To name one particular example, such question 
of principle - whether solid He^ is a supersolid at T = 0, and if not why and whether a bit different 
atomic potentials may still produce it. 

One of the flrst applications of the path integral method by Feynman himself were aimed at expla- 
nation of the He^ A-point. His classical papers ^123| 1124] introduced the idea of "polygon clusters". 
Starting from some conflguration of particles at Euclidean time r = 0, Feynman identifles such clusters 
as a group of atoms which exchange places during the Matsubara time t = (3 = h/T. The polygons 
of course in principle may have any shape, but since a minimal additional action for "jumping" atoms 
be required, the most probable paths are those where particles interchange places with their nearest- 
neighbors and thus the most probable jumps have length of one interparticle distance. And yet, the 
polygons themselves need not be small: as Feynman shown BEC implies that there should be infinitely 
large cluster. Furthermore, the BEC transition temperature can be deflned as such T at which the series 
over inflnitely long polygons start to be divergent. And, since the number of polygons grow with its 
length in a certain geometrical way, one may argue that there should be universal action per "jumping" 
particle AS** corresponding to all BEC transitions. 

Feynman's main approximation was the estimate of AS** using kinetic energy only, ignoring the po- 
tential one. He motivated by "rapid particle rearrangements" and wrote the thermal partition function 
in the form 



m*T 



e-f = j^T. j (^T'' M-'^T.iR^ - Pm']pm....RN)^'Rr...^'RN (53) 



where exponent comes from "jumps" done in a straight lines and containing corresponding action related 
to kinetic energy. To correct for potential energy somehow, Feynman introduced an "effective mass of a 
He atom" m* . The sum is done over permutation P of the particle coordinates, the function p includes 

^■^Some of quantum effects can be put into effective potential, see e.g. recent dedicated study related to charm quarks 
in sQGP [ng. 

^^ In fact "Bose-Einstein" statistics is due entirely to Bose while condensation entirely to Einstein. 
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Figure 23: Total extra action for "jumping" He4 atoms as a function of the temperature. 



effects of the interparticle interactions. His main idea was that the function p can be inferred from 
general properties of the liquid (or solid), its quasi-ordered local structure with peaked distribution 
over interparticle distances at some nearest-neighbor value d (in the case of cubic lattice d is the lattice 
spacing) . 

The relative amplitude of a term with permutation of n atoms is thus proportional to the n-th power 
of the "jump amplitude" 

m*Td'^- 



Hf = exp{—AS) = exp 



2h' 



(54) 



which should at the transition point be exactly compensated the divergence of the combinatorial pref- 
actors, describing a number of corresponding to non-crossing polygons on the corresponding 3d lattices. 
Feynman estimated that he expected the critical value to be in the range yc = 1/4 - 1/3. Kikuchi et 
al. |127] have argued that the critical action should be Sc ~ 1.9, but later studies in the eighties such 
as by Elser |129j determined numerically that the critical action is smaller Sc ~ 1.44. 

Let us see what critical temperature for He would follow from those actions, in Feynman approx- 
imation. We consider the distance d of a "jump" being fixed to the position of the nearest neighbor 
maximum in the static correlation function g{r) for liquid He, which is d « 0.35 nm . Kikuchi's critical 
action with the unmodified He mass leads to T^ = 3.57i^ while Elser's to T^ = 2.72K, still well above 
the correct position of the A-point T^ = 2.17 K. 

Cristoforetti and myself jl30i have tried to generalize Feynman criterion for interacting gases by 
including the potential energy effect "jumps" and their actions. We have used two approaches, (i) a 
semiclassical one, in which the path is found by minimizing classical equation of motion and the action 
given by this "instanton"; and (ii) numerical one, simulating path integral for a "jumping" particle in 
the potential created by "non-jumping" particles in the background. 

The setting for the calculation can be described as follows: imagine Helium atoms filling the best- 
packing lattice and consider all atoms to be non-jumping except of atoms in one raw of the lattice which 
are all "jumping" coherently by one step in the same direction. This is a good approximation because 
one should consider a very large polygon with infinite number of atoms. We found that many different 
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systems, from He atoms to Coulomb systems ^^ have effective potentials for jump which can be quite 
accurately be represented by universal periodic potential with only one harmonics V ~ cos{nx/a). It is 
simple to find the instanton path for this potential and its action. It improves substantially Feynman 
estimate, but since the value of the critical action (~ 1.5) is not really large, semiclassical estimates 
cannot be very accurate here and have mostly pedagogical value. 

In fact it is not difficult to perform numerical evaluation of the ratio of two path integrals, for 
"jumping" to "periodic" paths can be performed numerically: in this way also quantum fluctuation 
around the classical trajectory are taken into account. These calculations are very cheap in terms of 
computational power because we have a system with only one particle. We have done so for different 



values of the temperature, the result is presented in Fig. [23| The value of the jump action at T^ = 2.17K 
is Sc = 1.48, which is very close to the value Sc ~ 1.44 computed in [129] using the combinatorial 
method^ 

Solid He is a bit denser than liquid, and thus its action never reaches Sc ~ 1.44: thus - from the 
viewpoint of this criterion - it seems never become a "supersolid" . (The experimental situation is 
complicated with conflicting experiments about rotating solid He at very low T.) 

Let us now return to QCD near-T^ region. Before we return to our monopole, let me mentioned 
one school of thought (from |132] to |131j ) who think differently about confinement, putting 2d objects 
- vortices - as a primary objects randomizing the vacuum fields and leading to the area law for the 
Wilson loop. If so, the vacuum should have a condensate of vortices, which are after all also bosonic 
objects. However to my knowledge neither the order parameter - resembling Pisa one for monopoles - 
have been constructed, nor there is any understanding how one can address the issue of BEC for such 
objects, whose worldvolumes are 2-dimensional. 

The monopoles, on the other hand, are particles, and the Feynman criterion { S = S* aX T = Tc) 
discussed above should be applicable. Let us use it a la Feynman, using only kinetic energy first. The 
motion of the monopole which tries to occupy the position of another identical monopole at distance a 
away during the Matsubara time (3 = h/T = l/{.27 GeV) = .73 fm is relativistic, with the (Euclidean) 
velocity on the straight path v = a/ (3. Since we speak about tunneling, having imaginary action and 
velocity > 1 is in fact quite appropriate; no negative roots should appear. The Euclidean relativistic 
action is just the length of the path 

S, = rSu=nalds^rn0.JYT^- (55) 

should be S* ^ 1.44 at T = Tc. Since we know the monopole density (and thus a) and T^ we can use this 
criterion to estimate the effective monopole mass: this sets it in the range m = {0.2 — 0.3)GeV. Putting 
Coulombic repulsion and monopoles into appropriate cubic lattice, one can get sinusoidal potential for 
jumps and evaluate semiclassically or numerically the action for a "jump" more accurately. 

Since in principle lattice data on the whole path integrals of monopoles is known, in the future one 
should be able to learn the effective monopole masses and interactions in the near-T^ region and verify 
the details of their Bose-Einstein condensation. 

4.4 Summary of the magnetic scenario in the near-Tc QGP 

The main conclusion of this section is that lattice monopole-like singularities — located via endpoints 
of the Dirac strings in certain gauges such as maximal abelian gauge - are found to behave as physical 
particles. Their density is independent of the lattice details. The extracted magnetic coupling grows 
with T, exactly as expected. Their spatial correlation functions agree with the idea that their ensemble 

^^Let me remind the reader that below we hope to use this criterion of confinement for monopole plasma. 
^^The figure clearly shown a minimum at T = 1.94 K and then grows again but presumably the Feynman argument 
cannot be used at T < Tc. 
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is just a strongly coupled liquid driven essentially by the magnetic Coulomb interaction: they are quan- 
titatively reproduced by simple classical MD. In addition, the transport properties of such a magnetic 
plasma is in agreement with heavy ion phenomenology. 

But we don't yet know many important parameters, such as e.g. the masses of all quasiparticles 
involved. Current understanding of those masses is vague, the only direct lattice measurements have 
been done at T = l.ST^ and ST^ for quarks and gluons in |133j . both masses are about equal and 
M Ri .6 — .8GeV at both T. Let us assume that they are the same at T^ as well. 

We now proceed to self-dual dyons, a components of finite-T instantons. The mass is derived from 
the instanton action 

Mdyons = n%^) ~ 4r, ~ 1 GeV (56) 

where we used a typical instanton action to be around 12. Their density (for all Nc types together) can 
also be estimated from lattice data on instantons 

ndyons ~ Tfrj^ ~ 2 fm^ (57) 

So, although the dyon density is not negligible, it is not peaking near T^ and is few times smaller than 
that of monopoles. 

The monopole density numbers from Figures above does peak at Tc, indicating that they do play 
significant role in confinement. It is numerically very high - in absolute numbers it is about 

UmonoiT = T,) ~ 10fm~^ (58) 

This is larger than Bose gas of non-interacting bosons even with zero mass, but for strongly correlated 
liquid this is probably possible. Hopefully, the estimate of the monopole mass from the Feynman's 
Bose condensation criterion rrimono = (0.2 — 0.3)GeV will also work out. At the moment, all we know 
about magnetic plasma is its density and magnetic coupling. Those have been combined into Gamma 



parameter in the Fig, 19 (down right) with a conclusion that it is always above 1, so we do have a 
magnetic liquid. 

As explained above, we had classical MD simulations which lead to reasonable diffusion and viscosity 
of such plasmas: and a 50-50 mixture of electric and magnetic particles (occurring somewhere at 
T ^ 1.5Tc) is probably the optimal point for the lowest viscosity. How classical effects will modify 
these numbers we do not know yet. 

5 AdS/CFT duality 

AdS/CFT correspondence |134j is a duality between specific gauge theory known as Af=4: super- Yang- 
Mills theory (SYM) in 4 dimensions and the (10-dimensional) superstring theory in a specific setting. 
Before we describe the results obtained, here is a brief introduction, explaining its logic and what all 
these words mean aiming at non-experts. Those who know it may skip two next subsections. 

5.1 Black holes and AdS/CFT for pedestrians 

For applications we have in mind the most important feature of string theories is the existence of 
massless modes of closed strings, which include spin-2 "gravitons" as well as lower spin members of the 
multiplet, including vectors and scalars. String theory thus has aspiration to be "theory of everything" 
combining gauge and gravity interactions known in nature. This aspect of it is not discussed in this 
paper at all: all of these "bulk" fields in 10 dimensions will only play a role to mimic 4-dimensional 
QCD-like gauge theories. 
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One of several "string theory revolutions" of 1990's was caused by a discovery of solitonic objects of 
(nearly) any dimensions p (ranging p = —1 — 9) called Dp branes. One should think of them as some p- 
dimensional membranes, solitons made of strings and embedded in 9+1 dimensions of space+time. For 
our purposes the most important example is the D3 brane which has infinite extension in 3 coordinates 
xi,X2,X3, while it has no extension (and thus is just a point object) in the remaining 6 coordinates 

X4...Xg. 

Since branes are material objects with certain masses and charges, they create gravity/Coulomb/scalar 
fields around them in the bulk, like planets and stars do. As they have no extension in some coordinates 
(e.g. just mentioned X4...X9 coordinates for D^ brane), their fields are that of a point objects: in general 
relativity this naturally implies that they are black hole-like in such coordinates. As they are extended 
in other coordinates, the proper name of such objects is black branes. To find out their properties in 
classical (no loops or quantum fluctuations) approximation, one has to solve coupled Einstein-Maxwell- 
Scalar eqns, looking for static spherically symmetric solutions. This is no more difficult than to find 
the Schwartzschild solution for the usual black hole, which we remind has a metric tensor in spherical 

coordinates 

dv 

ds^ = g^^dx^dx" = -(1 - rhlr)dt^ + — + r^t/fi^ (59) 

(1 -rh/r) 

and the horizon radius (in fulll units) is rh = 2GjvM/c^, containing the mass M and Newton constant 
Gat. The horizon - zero of Qqo - is the "event horizon": a distant observer cannot see beyond it. (There 
900, grr change sign.) Similar expressions, with appropriate powers of the distance can be derived for 
p-branes. 

Note that at large distance r <^ Vh deviations from fiat metrics is a small correction - there 
the nonrelativistic potential description of Newtonian gravity is possible. Similarly far from a brane 
the fields are just given by Newton+Coulomb+scalar formulae in corresponding dimensions, and it 
is not hard to figure out what are their mutual interactions. Surprising (to anyone not involved in 
monopoles/supersymmetry research) the answer is that all these forces mutually cancels out, so there 
is no interactions between branes -they are so called BPS protected objects. As a result, a "brane 
engineer" may consider a number of parallel branes to be put at some random points: and they will 
stay there. Large (iV^ — > cxd) number of branes put into the same point combine their mass and thus 
create strong gravity field, justifying the use of classical Einstein/Maxwell eqns. Open string states, 
which keep their ends on some branes i and j lead to effective gauge theory on the brane with the U{Nc) 
group. 

When the black hole has a nonzero vector charge, the Gauss' law insists on constant flux through 
sphere of any radius, thus there are nonzero flelds at large r and Schwartzschild solution gets modifled 
into a "charged black hole". For asymptotically flat spaces there is the famous statement, much em- 
phasized by Wheeler: there are "no scalar hairs" of black holes, as there is no Gauss theorem to protect 
them. This happens to be not true in general, and for spaces which are AdS-like thus scalars should 
not be left oveiri 

The solution for D3 brane happens to be the so called extremal (6-dimensional) charged black hole, 
which has the lowest possible mass for a given charge. When the mass decreases to the extreme value, 
the horizon shrinks to nothing, and thus it is especially simple. Spherical symmetry in 6d allows one 
to separate the 5 angles (making the 5-dim sphere 5*5) from the radius r (in 6-dimensions orthogonal 
to "our world" space Xi, x — 2, x — 3) and write the resulting 10-d metrics as follows 



^^, ^ -dt + dx, + dx ^+dx, ^ ^ lyr^idr' + r'dni) (60) 

'1 + Lyr^ 



^^In fact there is a whole recent direction based on solution with a "scalar atmosphere" around black branes, providing 
gravity dual to superconductors on the boundary. 
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Again, at large r all corrections are small and we have the asymptotically flat space there. However 
at this point Maldecena told us to use instead further simplified version of this metric, in the "near- 
horizon region" at r << L, when 1 in the roots can be ignored. If so, in the last term two r^ cancels out 
and the 5-dimensional sphere element gets constant coefficient and thus gets decoupled from 5 other 
coordinates. Thus quantum numbers or motion in S^ becomes kind of internal quantum numbers like 
flavor in QCD, and it will be mostly ignored from now on. What is left is very simple 5-dimensional 
metric known as Anti-de-Sitter metric. Using a new coordinate z = L'^/r we get it into the "standard 
AdS^ form" used below: 

, , —dt'^ + dx'i + dxl + dxl + dz^ , , 

ds^ = ^ ^^ 5 61 

Note that z counts the distance from "the AdS boundary" 2; = 0. This metric has no scale and is not 
asymptotically fiat even at the boundary. Performing dilatation on these 5 coordinates we find that the 
metric remains invariant: in fact one can do any conformal transformation. It is this metric which is 
AdS in the AdS/CFT correspondence, and string theory in this background is "holographically gravity 
dual" to some conformal gauge theory at the boundary. 

So far nothing unusual happened: all formulae came straight from string and general relativity 
textbooks. A truly remarkable theoretical discovery is the so called "holography": the exact duality 
(one-to-one correspondence) between the 5-dim "bulk" effective theory in AdS^ to 4-dim "boundary" 
(r — i> 00) gauge theory. There is a dictionary, relating any observable in the gauge theory to another 
one in string theory: the duality implies that all answers are the same in both formulations. We will 
see below how it works "by examples" . 

The last step, which makes it useful, is the Maldacena relations between the gauge coupling, the 
AdS radius L and the string tension a' (which comes from the total mass of the brane set): 

L' = g'N^ia'f = X{a'f (62) 

It tells us that large gauge coupling A >> 1 corresponds to large AdS radius (in string units) and one 
can use classical (rather than quantum) gravity. At the same time the string and gravity couplings 
Qs ~ g^ may remain small: so one may do perturbative calculations in the bulk! 

At this point many readers are probably very confused by new 5-th dimension of space. One possible 
approach is to think of it as just a mathematical trick, somewhat analogous to more familiar introduction 
of the complex variables. 

(Suppose an Experimentalist measured some complicated cross section which is approximately a 
sum of Breit-Wigner resonances. His friend Phenomenologist may be able to write the answer as an 
analytic function with certain pole singularities in the complex energy plane, which will help for fitting 
and for evaluating integrals. Even better, their other friend Theorist cleverly developed a "bulk theory" , 
deriving the pole positions from some interaction laws on the complex plane. ) 

However, there is a perfectly physical meaning of the 5-th coordinate. One hint is provided by the 
fact that distance along it /^ dl = J^ dz/z = log{b/a) is the logarithm of the ratio. Thus its meaning 
is the "scale", the argument of the renormalization group. If one takes a bulk object and move it 
into larger z, its hologram at the boundary (z=0) grows in size: this direction thus corresponds to the 
infrared direction. The running coupling constant would thus be a z-dependent field called "dilaton" . 
Indeed, there are theories with gravity dual, in which this field (related to the coupling) does "run" in 
2;:unfortunately, known examples do not (yet?) include QCD! In spite of that, there are efforts to built 
its gravity dual "bottom-up" , introducing weak coupling at ultraviolet (small z) |139] and confinement 
in infrared (large z) |140[ I141J by certain modification of the rules. These approaches - known as 
AdS/QCD- we would not discuss in this review, except briefiy in the subsection on bulk viscosity. 

Let me briefiy remind few more facts from the Black Hole toolbox. Studies of how quantum field 
theory can sit in a background of a classical black hole metrics have resulted in two major discoveries: 
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the Hawking radiation and the Bekenshtein entropy, related to horizon radius and area, respectively. 
Hawking radiation makes black holes in asymptotically flat space unstable: it and heats the Universe 
till the black hole disappears. Putting black hole into a finite box also does not help: it is generically 
thermodynamically unstable and gets smaller and hotter till it finally burns out. Only in appropriate 
curved spaces black branes can be in thermal equilibrium with their "Universe", filled with radiation 
at some finite temperature T. As shown by Witten, all one has to do to get its metric is to consider 
non-extreme (excited extreme) black brane solution, which has a horizon. 

Let me now make a logical jump and instead of just giving the metric for finite-T solution let me first 
provide an "intuitive picture" for non-experts, explaining the finite-temperature Witten's settings in 
which mosl|^ pertinent calculations are done, shown in Fig,25 The upper rectangle is the 3-dimensional 



space boundary z=0 (only 2 dim shown), which is fiat (Minkowskian) and corresponds to "our world" 
where the gauge theory lives. Lower black rectangles (reduced in area because of curvature induced 
by 1/z'^ in metric) is the corresponding (same in x — l,a; — 2, x — 3) patch of the horizon (at z = z^) 
of a black hole, whose center is located at z = cxd. Studies of finite-T conformal plasma by AdS/CFT 
famously started exactly by evaluation of the Bekenstein entropy |138] . S = A/ 4 via calculating the 
horizon area A. 

Now comes the promised "intuitive picture" : this setting can be seen as a swimming pool, with the 
gauge theory (and us, to be referred below as "distant observers") living on its surface, at zero depth 
z = 0, enjoying the desired temperature T. In order to achieve that, the pool's bottom looks infinitely 
hot for observers which are sitting at some fixed z close to its coordinate z^. thus diving to such depth is 
not recommended. Strong gravity takes care and stabilizes this setting thermodynamically: recall that 
time units, as well as those of energy and temperature are subject to "warping" with qqq component of 
the metric, which vanishes at Zh- 

When astronomers found evidences for black holes and accretion into those, the physics of black hole 
became a regular part of physics since lots of problem have to be solved. Here important step forward 
was the so called membrane paradigm developed by many people and best formulated by Thorne and 
collaborators [136], known also under the name of "stretched horizons". Its main idea is to imagine that 
there is a physical membrane at some small distance e away from the horizon, and that it has properties 
exactly such that all the eqns (Maxwell's, Einstein's etc) would have the same solutions outside it 
as without a membrane but with a continuation through horizon. For example, a charged black hole 
would have a membrane with a nonzero charge density, to terminate the electric field lines. The fact that 
Poynting vector at the membrane must be pointed inward means some (time-odd!) relation between E 
and B: this is achieved by giving the membrane finite conductivity, which in turn leads to ohmic losses, 
heat and entropy generation in it. Furthermore, as shown by Damour back in 1980's, displacements of 
the membrane and relations for gravitational analog of the Poynting vector leads to nonzero viscosity 
of the membrane, and its effective low frequency theory take the form of Navier-Stokes hydrodynamics. 
For a bit more modern derivations of the effective action and field theory point of view look at Parikh 
nd Wilczek |137] . As we will see below in this chapter, all of those ideas have resurfaced now in AdS 
context, generating new energy of young string theorists who now pushed the "hydrodynamics of the 
horizon memebrane" well beyond the Navier-Stokes to a regular construction of systematic derivative 
expansions to any needed order. 

5.2 CFT for pedestrians 

The J\f=4: SYM theory is a cousin of QCD: it also has gauge fields with SU{Nc) color symmetry, but 
instead of quarks it has four "fiavors" of fermions called gluinoes, as their adjoint colors are the same 

^^The exception is heavy quark diffusion constant calculated by Casalderrey and Teaney[f2f^ which needs more com- 
plicated settings, with a Kruskal metric connecting a World to an Anti-world through the black hole. 
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as for gluons. There are also 6 adjoint scalars: with 2 polarizations of gluons it makes 8 bosonic 
modes, same as 2*4 fermions. This makes supersymmetry possible and leads to cancellations of power 
divergences. This theory is the most symmetric theory possible in 4 dimensions: it has conformal 
symmetry and its coupling does not run\ 

How do we know this? Of course, one had calculated first few perturbative coefficients of the beta 
function, and indeed see that negative gauge contribution is nicely canceled by fermions and scalars, 
order by order. But there are infinitely many coefficients, and one has to check them all! An elegant way 
to prove the case is based on another outstanding feature of the A/'=4 SYM: this theory is self — dual 
under electric-magnetic duality. As we discussed above, the Dirac condition requires the product of 
electric and magnetic couplings to be constant: and so in QCD and other gauge theories they indeed 
run in the opposite directions, electric becoming weak in ultraviolet and magnetic weak in infrared. But 
the multiplet of (lowest) magnetic objects of the A/'=4 SYM theory include 6 scalars (the monopoles), 
plus 4 fermions (monopoles plus one gluino zero mode occupied), plus 2 spin-1 (monopoles with 2 
gluinoes): this turns out to be exactly the same set of states as the original electric degrees of freedom 
(gluons-gluinoes-Higgses) . That means that an effective magnetic theory has the same Lagrangian as 
the original electric formulation: thus it must have the same beta function. Since two couplings cannot 
run in the opposite direction following the same beta function, they cannot run at all!) 

Recall that in QCD-like theories the scale Aqcd came about because of running coupling. If in the 
A/'=4 SYM theory the coupling constant does not run, it means that there is no analog of Aqcd in 
this theory, and since all the fields are massless and all the coupling dimensionless, thus the A/'=4 SYM 
theory has no scale at all. It is thus conformal field theory - the CFT in the AdS/CFT correspondence. 
One consequence is that the finite-T version of this theory (we will be mostly interested in) is the same 
whether T is large or small, since there is no other scale to compare with! This is similar to QCD 
plasma in the so called "quasi-conformal regime" : at high enough T > 2Tc all dimensionless ratios like 
(energy density) /T"^ are practically T- independent. 

Weakly coupled A/'=4 SYM theory can be studied perturbatively, hke any other gauge theory. What 
makes it unique is that AdS/CFT correspondence allows also to study it in the strong coupling limit, 
defined by a large value of the so called 't Hooft coupling, a combination of gauge coupling and number 
of colors which go together 

A = g^Nc » 1 (63) 

It is this combination which appears in physical effects, e.g. the Coulomb law. Thus, before we 
embark on studies of strong coupling regime A >>!, the reader have all reason to say: wait a moment, is 
it not followed from Klein-Gordon (or Dirac) eqns that for coupling larger than something two charges 
will fall at each other, as the square of the Coulomb potential ~ A^/r^ will dominate the centrifugal 
term /(/ + l)/r^? Zahed and myself [l42j worried about this, and their semiclassical approach to 
density of charge was even used later by Klebanov, Maldacena and Thorn |143] . When the spectrum 
of heavy quarkonia was eventually found from AdS/CFT, all states including the lowest s-wave ones 
were accounted for - they have small but positive masses, and string theorists were not surprised. But 
frankly I still don't understand why no falling actually happens. This is one of many puzzles which 
shows that what can be derived from gravity side may be very hard to understand in the gauge theory. 

5.3 The first example of AdS/CFT at work: new Coulomb law 

Let me start with our first example of the "AdS/CFT at work", related with the strong-coupling version 
of the Coulomb law calculated in |135] . The setting- to be called "the Maldecena dipole" - is shown in 



Fig, 24 'a), includes two static charges (heavy fundamental quarks) separated by the distance R. 



At weak coupling - the usual QED - we think of one charge creating the electric potential in which 
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the other is placed, leading to the usual Coulomb law which in our notation is 



ViD = -£ (64) 



is a sum of two Coulomb fields. 



p ( ^_ ( 9 .( Vm- {L/2)em Vm + {L/2)em \ , . 

^^yy) - ^4^^ 1^1^^ _ (L/2)e^|3 \y^ + [Ll^^eJ,^] ^'"^ 

Because of a cancellation between two terms, at large distances from the dipole the field decays as 
E ~ L/y^ . The corresponding energy density (and other components of the stress tensor) are thus of 
the order Tqo ~ g^l? jy^ , with certain "dipole" angular distribution. 

In the A/'=4 theory at weak coupling the only difference is that one can exchange massless scalars 
on top of gluons. It is always attractive, and for two heavy quarks leads to cancellation of the force, 
with doubling for quark-antiquark (we discuss now). The QED coupling g^ changes to 't Hooft coupling 
A = g^Nc proportional to the number of colors Nc- 

Now we turn to the AdS/CFT for the 7\/'=4 theory at strong coupling |144t 1145] . The electric flux 
in the bulk forms a singular object - the string (shown by the solid curve in Fig|24Fa)) - which pends 
from the boundary z = due to gravity force into the 5-th dimension, like in the famous catenary 
(chain) problempl The calculation thus follows from Nambu-Goto action for the string, whose general 
form is 

S = -^ / dadrJdetGMNdaX^dpX^ (66) 

where 2 coordinates a, t parameterize the string world line X^^(r, a), where M, N are space-time indices 
in the whole space (lOdim reduced to 5d in AdS/CFT). Gmn is the space metric and det stands for 2*2 
matrix with all a,p. In the AdS^ metric we need the components —Goo = Gu = ^55 = l/z"^, and we 
can think of a,T as our coordinates x,t: the string is then described by only one function z{t,x) and 
its action is reduced to 

S r^ f dtdx\^l + {dz/dxf - {dz/dty (67) 

We will use this action for "falling strings" below, and now proceed to further simpliflcations for static 
string, for which there is no time derivative and the function is z{x). Maldacena uses u{x) = l/z{x) 
and thus the Lagrangian becomes L = J{u^xY + "^^ with comma meaning the x- derivative. One more 
simpliflcation comes from the fact that x does not appear in it: thus an "energy" is conserved 

dL 

H = pq — L = — u X — L = E = const (68) 

which reduces the EOM from second order eqn to just (m,^)^ = u^{u^ /E"^ — 1) which can flnally be 
directly integrated to 

, , r du' , , 

xiu) = / , 69 

The minimum position of the string Um is related to E by the relation following from this formula at 
X = L/2. Plugging the solution back into action and removing divergence (which is independent of L) 
gives flnally the total string energy, which is the celebrated new Coulomb law at strong coupling 

-'-^^ - -nwx (™) 

^^ Another - more Einsteinian -way to explain it is to note that this is simply the shortest string possible: it is not 
straight because the space is curved. It is the same reason why the shortest path from New York to London does not 
look straight on the map. 
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The power of distance 1/L is in fact the only one possible by dimension, as the theory is conformal 
and has no scales of its own. What is remarkable is the (now famous) vA appearing instead of A in 
the weak coupling. (The numerical coefficient in the first bracket is 0.228, to be compared to the result 
from a diagrammatic re-summation below.) 

What is the reason for this modification? For pedagogical reasons let me start with two "naive but 
reasonable guesses" , both to be shown to be wrong later in this section: 

(i) One idea is that strongly coupled vacuum acts like some kind of a space-independent dielectric 
constant, e ~ 1/vA which is reducing the effect of the Coulomb field, similarly at all points, 
(ii) Perhaps such dielectric constant has nonlinear effects, and thus is not the same at different points: 
but the fields created by static dipole are still just the electric field E. 

A good feature of the AdS/CFT is that one can get many further details about the problem and test 
ideas like just expressed: but we cannot get all features. In order to understand the difference between 
dipole fields in a weakly and strongly regimes it be nice to calculate the electric and the scalar field 
distributions in both limits. Unfortunately we cannot do that: those fields do not belong to limited set 
of quantities one can calculate via AdS/CFT correspondence. In general AdS/CFT rules allows one to 
find holograms on the boundary of whatever happens in the bulk. This can be done by solving bulk 
wave eqns for massless bulk fields - which are scalars, vectors or gravitons. For example I will give 
our results for boundary stress tensor obtained from the (linearized) Einstein equation for the gravity 
perturbations, shown by the dashed line in FigJ24|(a). 

For example, in (relatively recent) paper Lin and myself |146] have found the stress tensor oi matter 
< T^u{y) > at any point y on the boundary, induced by the Maldacena dipole. The solution is too 
technical to be presented here and even the resulting stress tensor expression is too long^° let me just 
show only the leading terms far from the dipole y » L. 

Too = v/Xi= (^n^) ^"'* '''' 

where Ci, C2 are numerical constants whose values can be looked up in the paper and f{9) is the angular 
distribution shown by solid line in Fig|24Fb), quite different that in weak coupling (the dashed line). 
Note that in weak coupling the energy density from the dipole is just electric field ( [65| squared, leading 
to a different power Tqo ~ XL'^/y^. 

As we get glimpse of some first results from AdS/CFT we see that they are quite different from 
weak coupling and we would like to understand them. In fact we would like to do so both from the 
bulk (gravity) side as well as from the gauge theory side. It turns out the first is relatively easy. For 
example, both in the total energy and energy density we get a/A because this factor is in front of the 
Nambu-Goto Lagrangian (in proper units). The reason the field decays as y'^ is extremely natural: in 
the AdS^ space the function which inverts the Laplacean (analogously to Coulomb 1/r in fiat 3d) has 
that very power of distance 

Ps = ^ '' . (T2) 

with z being the 5-th coordinate of the source and r the 3-distance between it and the observation 
point 
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In order to understand the same results from the gauge side we will need a bit of pedagogical 
introduction: the resolution will be given by the idea of short color correlation time by Zahed and 



"^"^By he way, the stress tensor should always be traceless in conformal theory and have zero divergence: those conditions 
are used to verify explicitly that no mistake in the calculation was made. 

^"'^Please recall that in gravity there are no cancellations between different contributions: any energy source perturbs 
gravity with the same sign. 
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Figure 24: Setting of the Maldacena dipole: two charges at the boundary (black dots) are connected 
by the string (shown by sohd curve) pending under gravity toward the AdS center z —>■ oo. Classical 
graviton propagator (the dashed line) should be used to calculate the "hologram" - the stress tensor 
at the observation point y. The string is the gravity source; the point A has to be integrated over, (b) 
Angular distribution of the far field energy versus the polar angle {cos{6) = yi/\y\). Solid black line is 
the AdS.CFT result, compared to the perturbative dipole energy {3cos^9 + l)/4 (the dashed blue line), 
both normalized at zero angle. 



myself |142j . In QCD, with its running coupling, higher-order effects modify the zeroth order Coulomb 
field/potential. Since we consider static problem, one can rotate time into Euclidean time r = it, 
which not only makes possible lattice simulations but also simplifies perturbative diagrams. Let us do 
Feynman diagrams directly in the coordinate representation. The lowest order energy, given by the 
diagram in which one massless quantum (scalar or gauge component Ao) is emitted by one charge at 
time Ti and absorbed by another at time T2 is just 



ViLMTime) ~ -/ / 
Jo 



Time dTidT2 



g'^{Time) 



L2 + (n 



T-2 



(73) 



where 'Time' is total integration time and the denominator - the square of the 4-dim distance be- 
tween gluon emission and absorption represents the Feynman propagator in Euclidean space-time. The 
propagation time for a virtual quantum (ri — T2) ~ L, thus the Coulomb 1/L in the potential. 

Higher order diagrams include self-coupling of gluons/scalars and multiple interactions with the 
charges. A famous simplification proposed by 't Hooft is the large number of colors limit in which 
only planar diagrams should be considered. People suggested that as g grows those diagrams are 
becoming "fishnets" with smaller and smaller holes, converging to a "membrane" or string worldline: 
but although this idea was fueling decades of studies trying to cast gauge theory into stringy form it 
have not strictly speaking succeeded. It may still be true: just nobody was smart enough to sum up all 
planar diagrams ^^ 

If one does not want to give up on re-summation idea, one may consider a subset of those - the 
ladders - which can be summed up. Semenoff and Zarembo |148j have done that: let us look what 
have they found. The first point is that in order that each rung of the ladder contributes a factor 
Nc, emission time ordering should be strictly enforced, on each charge; let us call these time moments 
Si > S2 > S3... and ti > ^2 > ^3---- Ladder diagrams must connect Si to ti, etc, otherwise it is nonplanar 
and subleading diagram. Thus the main difference from the Abelian theory comes from the dynamics 
of the color vector. The (re-summed) Bethe-Salpeter kernel r(s,t), describing the evolution from time 



^Well, AdS/CFT is kind of a solution, actually, but it is doing it indirectly. 
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zero to times s, t at two lines, satisfies the following integral equation 

^(^■^' = i + 4^r''»r''*(737^r(M) (74) 

If this eqn is solved, one gets re-summation of all the ladder diagram. The kernel obviously satisfies the 
boundary condition T{S, 0) = r(0, T) = 1. If the equation is solved, the ladder-generated potential is 

VUL) = -l\m ;^r(T,r), (75) 

1 — >+oo 2 

In weak coupling F ?« 1 and the integral on the rhs is easily taken, resulting in the usual Coulomb law. 
For solving it at any coupling, it is convenient to switch to the differential equation 

"^" - '"'' nS.T) . (76) 



dSdT {S-Tf + L^ 
and change variables to x = {S — T)l L and y = {S + T)/L through 

Tix,y) = J2CrniU^)e--y/' (77) 

m 

with the corresponding boundary condition F(x, |x|) = 1. The dependence of the kernel F on the relative 
times X follows from the differential equation 

^ "^^ ^/'"'^7™(x) = -47-(x) (78) 



dx"^ x^ + 1 / 4 



For large A the dominant part of the potential in (78) is from small relative times x resulting into a 
harmonic equation |148] 

At large times T, the kernel is dominated by the lowest harmonic mode. For large times 5 ~ T 
that is small x and large y 

T{x, y) ^ Co e-^"'/^" e^^'^^ . (79) 



From ( 75 ) it follows that in the strong coupling limit the ladder generated potential is 



VUL) = -^ (80) 



which has the same param,etric form as the one derived from the AdS/CFT correspondence (70) except 
for the overall coefficient. Note that the difference is not so large, since I/tt = 0.318 is larger than 
the exact value 0.228 by about 1/3. Why did it happened that the potential is reduced relative to the 
Coulomb law by l/vA? It is because the relative time between gluon emissions is no longer ~ L, as 
in the Abelian case, but reduced to parametrically small time of relative color coherence r^ ~ l/LA^^^. 
Thus we learned an important lesson: in the strong coupling regime even the static charges communicate 
with each other via high frequency gluons and scalars, propagating (in Euclidean formulation!) with a 
super-luminal velocity v ~ A^^^ ^ 1. 

This idea -although it seemed to be too bizarre to be true - as we will see below to explain some 
of the AdS/CFT results. Klebanov, Maldacena and Thorn |143j have pointed out that the reason the 
stress tensor around the dipole is different from perturbative one by a factor ~ (L/r)/vA is actually 
explained by limited relative emission time by color coherence time. I am sure possible usage of this 
idea does not ends here. 

Before we leave the subject of Maldacena dipole, one more interesting question is what happen if one 
of the charges makes a small accelerated motion near its original position. Will there be a radiation? In 
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a somewhat different setting than used above, the answer was provided by Mikhailov [149j . He studied 
perturbations of the string and found that the radiated energy is described by familiar classical Lienard 
formula 



-2 rU -^nS 

°° X — \x X x\ 



AE = A ^^T-^^^ (81) 

•^-°° (1-f ) 

in which QED weak coupling constant A = |e^ is substituted by CFT strong coupling A = ^. This 



result is similar in its meaning to QCD synchrotron radiation we already mentioned [H] and would 
provide complete quenching of all jets beyond certain energy if the coherent fields in "glasma" would 
produce such acceleration, see discussion in 



5.4 Conformal plasma in equilibrium and the idea of relaxation 

In brief, the finite-T setting is different from zero-T AdS space discussed so far by existence of a horizon. 
Relaxation is due to the fact that all objects in the bulk sink toward it and their energy and information 
gets lost. 

After we have considered static heavy quarks and their strings, it is natural to proceed to the next 
stationary situation with a heavy quark, now in the finite-T setting. FigJ25ta) shows a setting of heavy 
quark quenching |15Ul I15H I152[ I153[ 1154] : a quark is being dragged (at some hight z^ related to the 
quark mass) by an "invisible hand" (to the left): its electric flux goes into the 5-th dimension, into the 
so called "trailing string". Its weight forces it to fall to the bottom (horizon). (Think of a heavy quark 
as a ship diligently laying underwater cable to the pool's bottom.) The cost of that is the drag 

v/2 



dP/dt = -nT'Jg^N, 



VT^ 
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connected to the diffusion constant via Einstein relation, a nontrivial successful check on two very 
different calculations. 

Another form of relaxation is studied via propagating "bulk waves" (b): massless ones may have spin 
S=0 (dilaton/axion),l (vector) or 2 (gravitons). Absorptive boundary condition at the horizon (black 
bottom)leads to spectra of "quasinormarj modes" with the imaginary part Im{uJn) ~ vrTn, setting the 
dissipation timescale of various fluctuations. 

We will give the equations themselves in the next subsection: let me just mentioned near- zero modes, 
corresponding to two propagating "surface waves", the longitudinal sound and transverse ^^diffuson". 
Absorption at the bottom (horizon) of both famously gives the viscosity rj/s = l/47r |120] . The waves 



may have real (timelike) 4-momentum or virtual (spacelike) one ^^ Rather complete spectra of quasi- 
normal modes and spectral densities for S=0,l,2 correlators are available, unfortunately mostly ex- 
tracted numerically. The case (c) - a "falling stone" - perhaps represent colorless (no strings attached) 
"mesons" , released to plasma and relaxing. 



The dashed lines in all Fig, 25 corresponds to the next-order diagram, describing back reaction of 
the falling bulk objects onto the boundary, the observation point denoted by a small open circle. This 
fields may also have spins 0,1 or 2, providing 3 pictures - the (4-d) "holograms" - of the boundary. 
Contrary to our intuition (developed from our limited fiat-world experience), the hologram is not a 
reduced reflection of more complete 5-d dynamics in the bulk, but in fact represents it fully. This 
phenomenon - the AdS/CFT duality - is a miracle occurring due to near-black-hole setting. 

Quasinormal modes are those which do not conserve the norm of the wave: it is hke decaying radioactive states in 
nuclear physics which are distinct from scattering ones, with real energies. 
^'^ This case, named DIS in AdS, is discussed here by E.Iancu. 
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Figure 25: Schematic view of the relaxation settings, a string (a), a wave (b) or a particle (c) fall into the 5-th 
dimension toward the black hole. 
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Figure 26: From |156] : hologram of the trailing string, the normalized energy density for one quark 
(supersonic jet) with v = 3/4 at nonzero T. 



These holographic images are what the surface observer will see. Image of the trailing string was 
calculated in |155tll56] : the recent example at nonzero T is shown in Fig,26|^b,c): it accurately displays 



hydro conical flow. For a hologram of the stone FigJ25|(c) see recent paper |157] : but to our knowledge 
the holographic "back reaction" of the falling waves remains to be done. 



How these predictions are related to experiment? Apart of those shown in Fig, 22, important test is 
whether the drag force indeed depends only on the velocity (rather than momentum): can be done via 
single electrons from c and b decays. Another challenge is to test if the effective viscosity 7]{k) is indeed 
decreasing with increasing gradients (or momentum k), as AdS/CFT nontrivially indicate. While in 
|158j only the issue of initial entropy production was discussed, experimentally better control of this 
effect can be inferred from studies of the viscous corrections to elliptic flow at more peripheral collisions, 
when the gradient in the direction of impact parameter k r^ R^ can be made large. 

Last but not least, the AdS/CFT correspondence predict a particular ratio of the excitation am- 
plitudes of the sound (conical flow) and diffuson mode. So far we seem to see experimentally only the 
conical flow, but not the second mode (which would result in a peak in the original direction of the 
jet). Much more detailed experimental and theoretical studies are required to see if indeed there is a 
potential discrepancy here or not. 



6 AdS/CFT and hydrodynamics 
6.1 Linearized hydrodynamics 



Son and collaborators |159] started their program of studies in linearized approximation: they wanted 
to see if small amplitude perturbations of the flnite T Witten's metric do in fact correspond to hydro- 
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dynamical modes: the sound and diffuson. The answer was positive, for example their famous value of 

viscosity turned out to be consistently the same in all tests they have made. This paper and its sequels 

are so clearly written that any reviewer can hardly hope to improve the presentation: so my advice is 

to look up these papers. 

Small perturbations are described by linearized Einstein eqns, and if time and 3-space derivatives 

are substituted by tu, k one has second order differential eqns in the "holographic coordinate" which for 

now we will denote u = z/zh- The three relevant eqns and their quasinormal spectra are to be found in 

Kovtun and Starinets |160] . the WKB approach to spectral densities for vector/stress tensor currents 

in Teaney |161] and also in another paper by Kovtun and Starinets [162]. The "master equation" has 

the form 

d^ d 

-^Za{u) + Pa{u)—Za{u) + qa{u)Za{u) =0, (82) 

where a = 1, 2, 3 labels the three symmetry channels, called shear, sound and scalar, respectively. Three 
sets of coefficients depend on the dimensionless frequency w = uj/2'n-T and momentum q = k/27rT and 

are 

. , (w2 - q2/)/ + 2u2w2 w2-q2/ 

^^("^ = ./(3w2 + q2(.2_3)) ' ^84) 

M/^(3w^ + q^(u^ — 3)) 

P3{u) = J- , q^W = 7^ , (86) 

where f = 1 — u^. The procedure is to look for solution satisfying the incoming wave condition at the 
horizon u = 1. Those can be written as a linear combination of two independent solutions with different 
behavior at m = 

Z,{u) = AaZi{u) + B,Zi\u) (87) 

and standard general argument express the retarded Green function in terms of the ratio of two coeffi- 
cients, now written in the original units 

2Ar2T^4-^a(^5^) (qq\ 



Ga{uj,q) = -7r'N^T 



Aa{uj,k) 



The quasinormal modes are poles of G or zeros ^ of A in the lower part of the frequency complex 
plane. 

To make the equations more familiar they may be easily put into the form of Schroedinger eqn. 
Let me just comment on some issues which may be confusing for non-experts. The issue I want to 
comment on is existence of the "horizon" in the thermal AdS metric at z = Zh, and from z > Zh 
no wave will return. The same situation also appears in quantum mechanics, e.g. for radioactive 
alpha decay problem. There are two ways to address the problem, which are orthogonal to each other: 
(i) to consider scattering formalism with real uj, k and multiple (but normalizable) scattering wave 
functions, and (ii) use of discrete "quasitationary states" which have complex u, k coming from "black 
wall" quantization condition. There is vast literature on how to use those: beware for example that 
if the sign of Im{uj) corresponds to decay of the system, Im{k) would produce a wave exponentially 



^^So that only better behaved solution Z[^{u) is present at u = O.We remind that a boundary condition at w = 1 was 
satisfied already, and thus it is standard quantization setting. 
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Figure 27: Sound dispersion (real and imaginary parts) obtained from the analysis of quasinormal 
modes in the AdS black hole background, from Ref. |160j . Gothic omega and q are frequency and 
momentum in units of 27rT. 



growing at large distances. So such states cannot be normalized, and need a lot of care to deal with. 
Quasinormal modes of black holes are also vast field, starting from original setting by Regge and Wheeler 
in 1950's for Schwartzschild solution, to AdS black holes in various dimensions: recent review by Siopsis 
is recommended jl64j . 

Let me give one of the examples of the Kovtun-Starinets paper, showing the real and imaginary 



part of th e low est quasinormal mode in the sound channel, see Figj27) we will return to its discussion 
in section 6.4, Note that at small momenta we have a standard sound, with velocity c^ = l/-\/3 and 
small dissipation described by bulk viscosity (thin solid lines) but it departs from this behavior for 
large q> 1. Although I do not show spectral densities of the appropriate correlators, they do contain a 
"sound peak" whose position is given by real part and width by the imaginary part of this quasinormal 
mode. Thus (at least the low frequency part of) bulk gravitons becomes the "gravity dual" to sound, 
in a very direct way. 



6.2 Bulk viscosity 

In the hydrodynamic discussions above we focused on shear viscosity 7] and now we will discuss recent 
progress related with the bulk viscosity C,. As their ratio C^/r] is known to be relatively small both at 
small and large T, one may naturally expect it to be maximal at or around Tc, when compressibility 
of matter is very small. The Kubo formula for C, includes correlator of pressures. Trace of the stress 
tensor or the dilatational charge can also be used (with care), and based on corresponding sum rules 
Kharzeev et al |165] predicted integral of the spectral density. These authors predicted a dramatic rise 
of the bulk viscosity near Tc, see also Karsch et al |166j who sugested that at Tc the peak value be as 
large as Q/ s ~ .3. 

From the phenomenological point of view,one may expect limits on the < C, > value (averaged 
over time in hydro evolution, in which the near-Tc region makes about 1/3 of the time at RHIC) from 
the same sources as limits on shear viscosity rj. One source of that is elliptic flow, which would be 
destroyed if < C > be too large: in fact the limits obtained by Romatschke |3S] and others is actually 
for a combination of r] and C,. While we expect r]/ s to be nearly constant and C,/ s peaked near T^ 
the latter gets a smaller weight in the time average: however since the QGP and "near-Tc" eras have 
similar duration at RHIC this factor is about (1/2) or so. The same argument applies to conical flow, 
and perhaps this condition is even stronger because (as we argued in our discussion of he conical flow 
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above) large angle implies that it is mostly formed exactly in the near-T^ region. Unfortunately, to my 
knowledge such phenomenological limits are not yet worked out in the quantitative form: my initial 
guess is that Karsch et al value just mentioned is about as large as possible. 

Violation of conformal symmetry in AdS/CFT setting would lead to (i) deviation of the thermo- 
dynamics, with nonzero (c^ — 1/3); and (ii) nonzero bulk viscosity (. Both effects were studied in 
|167[ 11681 1169] in a gauge theory deformed by a mass to J\f=2 super- Yang-Mills. Recently Gubser 
and Nellore |171] simplified the setting considerably, in the spirit of AdS/QCDrj scalar field with some 
tunable potential which basically mimics the beta-function of the desired theory. The relevant action is 



IottGs J 
and the corresponding equations of motion are 



'R-lid<pY-v{<p) 



^9) 



A0 = \/'(0) R^,, - -Rg^, = T^, (90) 



where the triangle stands for covariant Dalambetian and the stress tensor for the scalar is 

r^j,u = -d^^cpd^cp - -g^,u{d^T - ^g,.uV{4>) (91) 

The backgrounds of interest have the form 

ds^ = e^Mr) \-h{r)dt^ + dx''] + e2^('-) ^ $ = 0(r) (92) 

L -I h[r) 

The choice of radial variable r is arbitrary: re-parameterizing it leads only to a different choice of B. 

We will follow more recent paper by Gubser et al [172] , from which Fig|28] is taken. The idea of the 
method is explained at Fig. (a), which compares extraction of the shear (left) and bulk (right) viscosities: 
in both cases it is given by a probability of a reflected waves. The QCD-like potential is taken to be in 
th form 

V{(j)) = -—coshincj)) + 602 . (93) 

The parameter h can be adjusted so that the dimension A of the field theory operator dual to the bulk 
field (f) matches that of trF"^^ in QCD at a particular scale. 7, A can be tuned to get QCD-like behavior 
of the entropy shown in the top figure: the corresponding bulk viscosity are shown by blue lines in the 
two variants of the calculation shown in the bottom. The peak in bulk viscosity is indeed there, but 
it is about 5 times smaller than proposed by Karsch et al |166] . If so, the bulk viscosity is probably 
unimportant in comparison to the shear one at RHIC. 

6.3 Deriving gravity dual to (non-linear) hydrodynamics 

Janik and Peschanski |173] have proposed an approach which I call "top-down" : they proposed to take 
some well known hydro solution - e.g. the rapidity-independent Bjorken 1+1 dimensional expansion 
we discussed above - and look for its gravity dual by extrapolation. It can be said to be a combination 
of a Big Bang expansion in one "stretching" coordinate and a black hole in all others. 

They have used Fefferman-Graham coordinates, which is basically a gauge in which (755 = l/^^, g^i = 
0. In this gauge the so called "holographic renormalization" - the Taylor expansion of the metric near 
the boundary 

9AB = ^S + z'g^Il + z'g^^l (94) 

^^This large subject by itself is unfortunately outside the scope of this review: the reader is advised to look into |170j 
for a general introduction and other references. 
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Figure 28: (top)The normalized entropy density s/T^ as a function of T /Tc for two potentials of 
the form (pi) with parameters {7 ^ 0.606,6 ^ 2.06, A ^ 3.93} and {7 ^ 0.606,6 ^ 1.503, A Ri 3.61}. 



(bottom) The ratio of bulk viscosity to entropy density C,/ s for those two potentials, respectively. The 
result of the calculation is shown by the blue curve (the lower one, except at small T). The other (black) 
curve shows the upper bound suggested by Buchel. 



has the most direct physical meaning: the unperturbed term g'^^^ is just flat Minkowski metric, g^^f} = 



yfj,i/ 



< T^u > is related to the induced stress tensor on the boundary. If the 



on general grounds and g^^J 
latter is consider known, from some assumed solution to relativistic hydrodynamics, standard argument 
from the theory of differential equations ensures that for a second order differential equation a knowledge 
of 2 subsequent Taylor coefficients is equivalent to knowledge of all of them since the equation itself 
provides the recursive relation. So one can in principle extrapolate from z = into finite z, recovering 
the whole "gravity dual" solution in the bulk. 

Writing the metric as a late-time expansion of the form (schematically, with different functions for 
different pertinent components) 



9^^u 



tf'expE( 



n2/3 



y^Cniz'/r'/')] 



(95) 



with coefficients depending only on a single scaling variable z^/t'^^^, Janik and collaborators have put 
this series into Einstein equations and got multiple coupled eqns for these coefficients. The leading order 
solution indeed has a departing horizon located at Zh ~ r^/^, in agreement with the entropy conservation 
in ideal hydro. The next subleading terms 0{t~^^^) has been calculated by Sin and Nakamura |174j 
who identified them with the viscosity effects. The viscosity value was however only fixed by still 
further term calculated again by Janik, following the following principle: all physical observables - the 
invariants made of curvature tensors - should be nonsingular at the horizon. However still the next (the 



third) term calculated in |175] had bad logarithmic singularity which produced unphysical imaginary 
part in curvatures beyond the horizon. This singularity indicated that something is serious wrong with 
the asymptotic solution expressed as series by Janik et al. Although appearance of some singularities 
in the bulk are to be expected on general ground, one usually hope (as the next subsection would show 
directly) that all singularities are hidden from the distant observers by corresponding horizons, at which 
locally nothing happens to curvature invariants. (I argued even before that something is wrong already 
with Fefferman-Graham coordinates, as they cannot show metric beyond the horizon.) 

Very recently Heller et al |176] have argued that the problem was somehow induced by Fefferman- 
Graham coordinates, and in a different (Eddington-Filkenstein) coordinates, a non-singular asymptotic 
solution for Bjorken flow is naturally obtained. Explicitly the metric found is up to second order in 
gradient expansion. The reason why it works in this settings is better explained from the opposite 
approach to bulk "gravity dual" , which we will discuss now. general setting. 

"Down-up approach" to derivation of hydro is based on two papers by Bhattacharayya et al |178t 
I179J . see also |180j . who have worked out very intuitive and general derivation of hydrodynamics starting 
from (a parameterized) horizon singularity and solving Einstein eqns toward the boundary. I view it as 
the AdS/CFT incarnation of the old "membrane paradigm": the details are different (e.g. conformal 
symmetry now forces the bulk viscosity and many similar coefficients to vanish). 

Before we get technical, let me explain its quite intuitive idea for pedestrians. Recall our finite-T 
setting, which I compared to a pool heated from its floor. In equilibrium this floor was flat, located at 
the depth Zh which defined the temperature at the surface z = 0. Now we generalize this problem into 
a situation which appears in case of tsunami , when the ocean bottom is moving and also not flat. It 
is intuitively clear that it is easier to get waves on the surface starting from the known position of the 
bottom, rather than recalculate what happened on the bottom from the observations on the top. 

Simplification of the situation, allowing for a systematic study, is obtained due to the assumption 
of small derivatives (in 3 physical coordinates plus time) 

l^^^l^l^^— He«l (96) 

In the zeroth order in e one obviously have a simple "tubewise" approximation, in which T{t,Xi) on 
the boundary is given just by the bottom position at the same time and directly below the observation 
point. The authors worked out two next approximations, the orders e, e^ , corresponding to Navier- 
Stokes and the second-order hydrodynamics, respectively: apart of technical complexity of expressions, 
nothing seems to prevent calculation of further terms of this expansion. 

Mathematically speaking, the Einstein tensor in 5d has 15 components (thus equations), the metric 
after fixing the gauge has 10 components. As we restricted zeroth order solution to an ansatz, the 
first order (in derivatives) part of metric is unknown functions which is found from 10 eqns. The 
rest are "constrains" which should be valid without any freedom left: thus 4 functions in the original 
ansatz (the local temperature and flow velocity T{x^),u^{x^), remember u^u^ = 1) cannot be arbitrary. 
Remarkably, the constraints happen to be (zeroth order) hydrodynamical eqns! 

Let us return to "equations" for metric: as usual there are three versions of the z-eqns (for 4d 
spin-0,1,2 of the 5d gravity). Keeping derivatives of zero order ansatz as sources, all of them can be 
actually solved by straight integration for arbitrary sources. The starting point is the form of the metric 



the 'boosted black branes^^ 



ds'^ = -2 u^ dx^dr - r^ f{br) u^ u^ dx^'dx" + r^ P^,^ dx^'dx" , (97) 

with 

!(^)-'-7'- ^'-7T^'' ^'^TTh' <'*> 



^^The indices in the boundary are raised and lowered with the Minkowski metric u^ — rj^^, u^ . 
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where the temperature T = ^ and velocities jSi are all constants with jS"^ = (3j (3\ and 

P^''^ = u^'u" + 7]'"' (99) 

is the projector onto spatial directions. The zero of ^foo is so-to-say proto-horizon: the true event horizon 
is to be discussed below, after the metric is determined. 

If all 4 parameters - T and velocity - are promoted into slowly varying functions of space and time, 
the Einstein eqn are no longer automatically satisfied: yet it remains a good zeroth order approximation. 



with systematic corrections'^^ organized in power of e. At the end, a set of equations one get for the 
functions turned out to be the fluid dynamics eqns on the boundary. 

Let me omit the solution itself: just saying that second order eqns in holographic coordinate r are as 
usual classified into 3 classes, with spin 0,1,2, each solved for arbitrary "source terms" coming from all 
those derivative terms which Einstein eqns generate. The general second order expressions are rather 
lengthy: let me just give for illustration the explicit global metric to first order in boundary derivatives 
about y'^ = 

ds^ = 2dvdr- r'^fir) dv^ + r^ dxi dx' -2x^ d^jSf^ dr dx' -2x^d^ pfr'^il - /(r)) dv dx' 

- 4 ^ — dv^ + 2 r^ F(r) a^f dx' dx' + -r dipf^ dv^ + 2r 8,13^ dv dx\ (100) 

where 

F{r) = I dx —, ^— ^^ r = - 

x(x + 1) (2;2 + 1) 4 



r 



, /(l + r)2(l + r2)\ 

In \ — - - 2 arctan(r) + vr 



(101) 



This metric solves Einstein's equations to first order in the neighborhood of x'^ = provided the 
functions 6^°^ and pl satisfy 

9,M°) = ^, d.h^'^ = d,fif\ (102) 

Given g^^^ ^ it is in principle straightforward to use the AdS/CFT correspondence to recover the stress 
tensor. Transforming the metric to Fefferman-Graham form and taking their boundary limit z ~ 1/r — ^ 
0, one can finally read off the desired T^^^ 

T^^ = \-U u^v" + r/'^n -%a^\ (103) 

where the last term contains the first order term which -as expected -has the Navier-Stokes shear 
structure 

o 
with known viscosity-to-entropy ratio. 

The same story is repeated in the second order, with 1-st order hydro eqns needed for the bulk 
solution and the stress tensor on the boundary containing the second order terms 

Tf"" = {-nTY {t]'"'' + 4 u^^u") -2{-nTf a'"' (105) 

+ (vrT)2 ((In 2) T^ + 2T,r + (2 - In 2) [i Tr + T^^ + T^]) (106) 

(107) 

where 



^®In this respect the calculation is similar to a derivation of the effective chiral Lagrangian in its spirit. 
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If, = eafijf,u°'d^u'^. 

The second order coefficients agree with those found previously by Baier et al |177j who worked out 
hnearized approximation. Haack and Yaroni [181j generahzed this analysis for conformal fluids in di- 
mensions 3jDj7. Thus, in a way we know a bit more about the CFT plasma than about such usual 
liquids as water! 

The second paper |178] of this (somewhat enlarged) group of authors provides further important 
developments and clarifications. One of them is the issue of finding the true dynamical horizon of the 
solution - which in dynamical situation is of course not located at zero of the 5^00- (This can be seen 
from the fact that true horizon must be a null surface.) Its position can be found also via e expansion, 
in a quasi- local form. I will not here go into real calculation of it, but following the authors themselves 
explaining the issue with a simpler example. The so called Vaidya metric describes massless matter 
accreting into the usual (3+ld) black hole 

ds^ = -(1 - '^'^)de + 2dtdr + r^dn'' (108) 

r 

The horizon is at some position r = r/j(t) which must define a null surface: this condition means that 
it must satisfy the following eqn 

rh{t) = 2m{t) + 2rh{t)rh{t) (109) 

Without the last term we have static Schwartzschild black hole with the usual Vh = 2m. With time- 
depending accretion, it is a differential eqn, so in general the horizon position is a global entity depending 
on the whole history of accretion process m{t). However this is drastically simplified in the case of slow 
accretion with rh = 0(e), to the extent that one can write down a local ansatz 

rh{t) = 2m{t) + am{t)fn{t) + hm{t)m{tf + cm{tfm{t) (110) 

with a, 6, c being just numbers fixed from the eqn to be 8,64,32, etc. The second 0(e) term is positive: 
thus the event horizon is above the zero of the qoq and calculation of its area is not affected by it. In a 
general case, with a horizon slowly depending on all 4 coordinates, the formulae are more involved but 
the principle is the same. 



A qualitative picture of evolution of the perturbed horizon is given in the Fig. 29 with time ripples 
on the horizon disappear 

The last issue they successfully addressed is that of entropy production. Even in dynamical situation 
one expects that the Bekenstein entropy is still ultimately is given by the horizon area, as it is a UV 
concept and cannot depend on soft variation of the horizon shape. Since the horizon surface was located 
explicitly, one may calculate this area. Furthermore, one may check the local form of the second law 
of thermodynamics, namely that at any spatial point entropy's time derivative is strictly non-negative. 
This is indeed demonstrated to be the case, to the e^ order solution obtained. 

(That can be of course be checked from the second order hydro eqns: this is the usual form of 
entropy production dissipative terms. What is shown in |178] is more than that: the positivity of 
entropy production is true not only in our world, on the pool's surface, but all the way to the "hell" 
with its infinite temperature and even locally at any point!) 

We have pictorially argued in the preceding sections that gravity dual to dissipative relaxation 
processes means falling into the gravitational abyss. Now, for pure gravity falling, |178] found how to 
trace the entropy production locally in time (and z coordinate), relating it to r^ and its derivative terms 
- the declined/time dependent bottom. Similarly one may ask how the entropy is produced if other 
forms of bulk matter - e.g. classical strings - are falling as well. 
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Figure 29: (from |178j ).The event horizon r = rnix^) sketched as a function of the time t and one of 
the spatial coordinates x (the other two spatial coordinates are suppressed). 

6.4 Entropy production in hydro: beyond the derivative expansion 

With derivative expansion being now under theoretical control, one may try to go beyond it and think 
about some kind of re-summations of all dissipative effects. This is not yet done: but let me provide 
two different reasons for doing that, one based on phenomenology and one on theory. 

Phenomenology of hydro, discussed at the beginning of the review, hints that its applicability starts 
very early, at a time of about 1/2 fm/c at RHIC At this time the longitudinal derivatives are not 
formally small enough yet to make the dissipative (Navier-Stokes etc) terms negligible even for the 
smallest rj/ s = l/47r; and yet the ideal hydro seem to do a good job even keeping the entropy. Can it be 
that the higher order terms somehow compensate the first (viscosity), making its total effect smaller? 

The theory issue pointed out by Lublinsky and myself |158] , comes from the behavior of the quasi- 
normal modes. Although there is nothing linearized about the early stages of heavy ion collisions, let 
me still use this example to illustrate the issue. 

Writing down the Navier-Stokes eqn and solving it for sound dispersion produces the following 
dispersion relation 



u; = ^±^c2P-r2fcV4 (111) 

where Tg = ri/{e + p) is the so called sound absorption length. Reliable hydro region is of course at 
small k, so one should disregard the last higher order term under the square root and get the textbook 
answer. Let me for pedagogical reason keep this last term and see what it does: it curves the path of 
ui{k) on the complex plane so that that left and right moving sound waves make half circle and collide 
at a finite k when the square root vanishes: apparently there is no sound beyond that point, the Navier 
and Stokes would say. 

What we just did is of course not justified because higher order gradient terms were not included. In 
the preceding section we discussed the second order gradients recently determined: it turns out those 
are 0{k^) real correction for the sound: they are increasing its speed a bit but keep sound dissipation 
unchanged. For the diffusion (shear) mode there is a correction to the imaginary part 
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.2 , (2 - ln2) , , 



Imuj{k) = {l/4)k' + ^ ^r + ... (112) 



^^We remind that this expression is written in units of ttT for w, k. 
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Figure 30: (Fraction of entropy produced during the hydro phase as a function of initial proper time. 
The initial temperature To = 300 MeV. The left (blue) points correspond to the first order (shear) 
viscosity approximation. The right (red) points are for the all order re-summation. 



which has the same sign as the first term, increasing the dissipation as k (or gradients) grows. 

Thus both examples seem to tell us that dissipation grows with k monotonously, even more so when 
we try to include corrections. Yet it is not clear how Imuj{k) depends on k away from the origin, or 
even if it is growing monotonously. Condensed matter physics is full of such examples, with peaks in 
the dissipation at some characteristic momenta/frequencies. 

In fact the lowest quasinormal mode from AdS/CFT does turn around, see FigJ27|(b) and its imag- 
inary part is reduced after the maximum at k/2'KT ~ 1. This means that the expansion in derivatives 
should have a radius of convergence qjl-nT ~ 1. 

The exercise Lublinsky and myself did was to introduce effective momentum-dependent viscosity 
ry(/c) and assume that it is decreasing with momentum k, reproducing the path of the lowest quasinormal 
modq_J We found that changing constant viscosity to such effective one turns out to be very important 
at very early stages, see Fig, 30 While the first order effects (left) more than double the initial entropy 
of the flow, the modified viscosity (right) produce much better controlled hydrodynamics and produces 
an order of magnitude smaller amount of entropy. 



7 AdS/CFT out of equilibrium 

7.1 Black hole creation: qualitative ideas 

The most challenging frontier in the AdS/CFT-based studies are those attempting to use it out of 
equilibrium, addressing initial equilibration and entropy production. Non-equilibrium state between the 
collision moment and ("reasonably") equilibrated QGP is now called "glasma" - which so far is modeled 
by random classical glue by Venugopalan and collaborators, via classical Yang-Mills equations in weak 
coupling. However the corresponding "saturation scale" Qs at RHIC is only about 1-1.5 GeV - not far 
from parton momenta in sQGP, the perfect liquid as we now know - so one may wander if a strongly 
coupled regime should be tried instead. This is what I propose to call "sglasma" frontier. 

If AdS/CFT is the tool to be used, one has to start with high energy collision inside cold T = AdS^ 
(the vacuum, or a bottomless pool) and then dynamically solve two difficult problems: (i) explain how 



and with what accuracy "the collision debris" act like a "heater" imitating black/hot patch of Fig 25 
(ii) find a consistent solution with "falling bottom", Zh{time), and find its hologram describing hydro 

^°For clarity: we don't understand its physics and do not claim it is the case: we did it to see what possible effect this 
would have, if true. 
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Figure 31: Schematic view of the colhsion setting. Setting of the sGLASMA studies: (a) a single pair of heavy 
quark jets, moving with velocities v and —v and creating falling string. Multiple strings create a 3-d falling 
membrane (2d shown), which is (b) first far from trapped surface and then very closed to it (c). 



explosion/cooling. 

Since the temperature/entropy only appear when a horizon is dynamically created, leading to the 
information loss, one may say that "falling debris" upgrade the extreme black hole into a non-extreme 
one. 

A lot of work has been done on gravitational collapse in asymptotically flat 4d space, for collapse 
into real black holes in our Universe. With modified multidimensional gravity, people are thinking 
about their possible formation in LHC experiments, and this also created extensive literature. In the 
AdS/CFT language we are sure that black hole production is not a rare event: in fact it must happen 
in each and every RHIC heavy ion collision event, with an effective gravity (imitating QCD) in the 
imaginary (unreal) 5-th dimension. 

In heavy ion context. Sin, Zahed and myself |182] first argued that exploding/cooling fireball on the 
brane is dual to a departing small black hole (separate from the AdS center), formed by the collision 
debris and then falling toward the AdS center. A specific solution discussed in that paper was a brane 
departing from a static black hole, which generated a "spherical" solution (no dependence on all 3 spa- 
tial coordinates) with a time-dependent T (which however is more appropriate for cosmology but not 
heavy ion applications). We also discussed several idealized settings, with d-dimensional stretching, cor- 
responding for d=l to a collision of two infinite thin walls and subsequent Bjorken rapidity-independent 
expansion, with 2d and 3d corresponding to cylindrical and spherical relativistic collapsing walls. 

Before we embark on some details, let me explain what (I think) should be done to solve this 



problem, schematically shown in Fig, 31 The first figure (a) shows the first step - a string which 
belongs to expanding dipole, with charges moving away from each other. The string solution and the 
corresponding hologram will be discussed in the next subsection: let us now just note that in this case 
no horizon is produced and thus no temperature, entropy or hydrodynamics are expected. The second 
figure (b) display multiple strings with departing ends: those we think can be a good approximation for 
heavy ion collisions, at least if we think of heavy ions as being made of heavy quarks. If the combined 
mass of all the strings (or other collision debris) is large enough, their gravity will induce a bubble 
of trapped surface (a part of is shown in the lower part of figure (b)). From the viewpoint of distant 
observer, the two membranes - the falling debris and rising horizon ones - will be after some time glued 
and moved together, as shown in figure (c). This is what I call the two membranes paradigm, to which 
we turn in the last subsection. 

Completing this introductory subsection, let me mention a separate direction by Kajantie et al |183j 
addressing the same issues in the 1+1 dimensional world. It is easier to work out math in this case: 
and one can indeed see how a black hole is produced. However, since the shear viscosity is absent in it 
because of 1+1 dimensions and the bulk viscosity is prohibited by conformity, so it is a cute toy case 
without any dissipation. 
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Figure 32: The hologram of a fahing string. The contours show the magnitude of the energy density (a) 
and the Poynting vector T°* in the transverse plane (b) at some time. The direction of the momentum 
flow is indicated by arrows. 



7.2 Falling strings and "jets" at strong coupling 

Rather early in development of QCD, when the notion of confinement and electric flux tubes - known 
also as the QCD strings - were invented in 1970's, B.Andersen and collaborators |184] developed what 
gets to be known as the Lund model of hadronic collisions. Its main idea is that during short time 
of passage of one hadron through another, the strings can get reconnected, and therefore with certain 
probability some strings become connected to color charges in two different hadrons. Those strings 
get stretched longitudinally and then break up into parts, making secondary mesons and (with smaller 
probability) baryons. Many variants of string-based models were developed, and some descendants -e.g. 
PYTHIA - remain widely "event generators" till today. 

AdS/CFT version of the Lund model has been developed in two works by S.Lin and myself |1851I146] . 
Before we describe some of its results, let me mention another possible usage of those results. Some of the 
dreams high energy theorists have about possible discoveries at LHC or beyond is existence of "hidden" 
interactions which we don't see because the so called "mediators" - particles which are charged both 
under Standard Model and hidden theory - are heavy. Perhaps the hidden sector is strongly coupled: 
if so one may wander how production of a pair of new charges would look like. Further discussion of 
this issue (and more references) can be found in a paper by Hofman and Maldacena |186] . who made 
further steps toward understanding the "strongly coupled collider physics" . 

Let us now return to [185^ 1146] and consider two heavy quarks departing from each other with 
velocities ±f and connected by a flux tube (classical string): in AdS setting it is not breaking ^^ but 



rather falling into the 5-th z direction, see FigJSlfa). Its action is the familiar Nambu-Goto action 
in AdS^ background, and if one ignores two transverse coordinates X2,X3 and uses as two internal 
coordinates the t, x (time and longitudinal coordinate) the string is described by by one function of two 
variables z{x,t). The corresponding string action is t 



S 



R' 



2-Ka' 



dt 






;ii3) 



^"'^This is so for classical string minimizing the action. However account for fluctuations of the string allows it to touch 
the flavor brane and break: in fact Peeters, Sonnenschein and Zamaklar [147J have calculated the holographic decays of 
large-spin mesons this way. 
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Before solving the corresponding equation in full, we will first discuss "scaling" solutions in the separable 
form 

^(r,r^) = y^ (114) 

where r, rj are the proper time and space-time rapidity we discussed when we discussed Bjorken solution. 
This form, suggested by conformal properties of the theory, were used in literature |187] , in a different 
- Euclidean - context, for AdS/CFT calculation of the anomalous dimensions of "kinks" on the Wilson 
lines (of which our produced pair of charges is one). The corresponding solution was obtained, but we 
found that it can only exists for sufficiently smal|__] velocities. Moreover, the analysis of the classical 
stability of such scaling solution revealed that it gets unstable for Y > Ym ~ 1/4, where quark rapidity 
is related to their velocity in the usual way v = tanh{Y). For larger velocity of the quarks the scaling 
solution has to be substituted by a non-scaling one, depending on both variables in a nontrivial way, 
which was analyzed numerically. 

The physical picture an observer at the boundary would see is again given by a (gravitational) 
hologram of the falling string calculated in the second paper |146j . One feature of the result should 
not be surprising for the readers who followed the description of the hologram of the static Maldacena 



string above: indeed, no trace of a string is in fact visible. The results - shown in fig, 26 a) are near- 
spherical explosion. This means that there are no jets at strong coupling] Instead of giving rather 
lengthy formulae let me just comment that we found this explosion to be non-thermal and thus non- 
hydrodynamical, in the sense that the stress tensor found (although of course conserved and traceless) 
cannot be parameterized by the energy density and isotropic pressure. 

7.3 Entropy production and the "double membranes paradigm" 

As it was explained in the previous section, "top-down" |173] and "down-up" significantly clarified 
how hydrodynamics can be derived in AdS/CFT. But we also would like in principle to understand 
how initial equilibration happened and the "dynamical horizon" get formed: this means to solving the 
gravitational collapse problem following Einstein equations from some initial "debris" produced by the 
collision all the way to black hole. 

The initial question is what can an appropriate representation be of the colliding nuclei. One 
straightforward approach is by Romatschke and Grumiller |188j who literally used the "shock waves" 
which holographically leads to the delta-funciton like boundary T^,^. Their metric for one shock wave is 

(115) 

note that it grows in into the z direction indefinitely, thus collision is that of "icebergs" which have 
their main weight deep down. Not surprisingly, the solution is developing in time "bottom-up", from 
infrared to the boundary, with boundary T^iy growing with the proper time. The solution is first found 
analytically at small proper time, expanding in its powers. Perhaps, if the authors would be able to 
do the second part of their program - solve numerically Einstein equations for finite proper time, they 
will be able to see approach to equilibrium and eventually transition to hydrodynamical expansion and 
cooling. Even then, the solution is not expected to go to Bjorken fiow because the setting is not rapidity 
independent. 

A significant leap forward had been done recently by Gubser, Pufu and Yarom |189j , who proposed 
to look at heavy ion collision as a process of head-on collision of two point-like black holes, separated 



from the boundary by some depth L - tuned to the nuclear size of Au to be about 4 fm, see Fig, 33 By 



using global AdS coordinates, these authors argued that (apart of obvious axial 0(2) symmetry) this 

^^In particular, for velocities going to zero one finds a strongly coupled version of Ampere's law, the interaction of two 
nonrelativistic currents. 
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*x^ 



Figure 33: From |189j :A projection of the marginally trapped surface that we use onto a fixed time 
slice of the AdS geometry. The size of the trapped surface is controlled by the energy of the massless 
particles that generate the shock waves. These particles are shown as dark blue dots. 



case has higher - namely 0(3)- symmetry with the resulting black hole at the collision moment at its 



center, thus in certain coordinate 



AzL 



(116) 
(Here 



the 3-d trapped surface C at the collision moment should be just a 3-sphere, at constant q = Qc- 
x± are two coordinates transverse to the collision axes.) The picture of it is shown in FigJ29|^b) 

If so, one can find the radius at which it is the trapped null-surface and determine its energy and 
Bekenstein entropy. For large Qc these expressions are 



2^3 



E 






3^2 



s 



~g: 



;il7) 



from which, eliminating qc, the main result of the paper follows, namely that the entropy grows with 
the collision energy as 

S ^ E^/^ (118) 

Note that this power very much depends on the 5-dimensional gravity and is different from the 1950 's 
prediction of Fermi/Landau, that this power should be 1/2. Simplistic comparison with data ignore the 
nonzero baryon number, important especially at lower collision energies, which should be of course be 
removed from the produced entropj]^ 

While the model is far from being realistic, the math of this paper is extremely interesting. First of 
all, the "depth" of the colliding objects should not at all be related to the nuclear size (which can in 
principle be taken to be infinite, for two colliding walls) , but rather to the typical "saturation scale" of 
parton's momenta. Second, this paper does not try to answer the most difficult issue - the dynamical 
formation of a horizon. Indeed, both colliding objects are already black holes, with nonzero temperature 
and entropy: at the collision time those are just joint together into a new shape. And, last but not least, 
it is not clear whether the resulting black hole will retain the same "spherical" shape of the horizon as 
it falls into the AdS center: if not, further entropy is produced. All of that would of course be subject 



of subsequent works: whether the result (118) will survive or not remains to be seen. 

Further work toward a more realistic "gravity dual" to heavy ion fireball is ongoing. A sketch of 
our current thinking has already been presented in Fig|3Tj If many strings are falling together their 
combined gravity is non-negligible - they are partly falling under their own weight. So one should solve 
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This point was clarified after this question was asked by E.Kiritsis at some meeting. 
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nonlinear Einstein eqns, which tell us that (from the viewpoint of distant observer) extra weight may 
actually slow down falling, eventually leading to near-horizon levitation. The trapped surface is moving 
first upward (shown at the bottom of FigjSlJ^b)) toward the falling membrane, till two collide, get close 
and fall together, see FigjSltc). After that distant observer finds a thermal hydrodynamical explosion 
as a hologram. This is the case at mid-rapidity but never in the fragmentation regions. 

So far we cannot solve it in realistic geometry, and used a simplifications instead. Lin and myself 
|191] considered the case when falling shell (or membrane, made of collision debris) is flat {xi,X2,X3- 
independent). In this collapsing shell case one finds a quasiequilihrium solution: the metric above the 
falling membrane is static thermal AdS in spite of the fact that the membrane is falling. We derived 
and solved equation of motion of it - from the so called Israel junction condition- and the metric below 
the falling membrane, which is simple (vacuum) AdS. 

The main question is by which experiments an observer on the boundary can distinguish true 
thermal state from "quasiequilibrium" . What we found is that a "one-point observer" would simply 
see equilibrium pressure and energy density, while more sophisticated "two-point observer" who can 
measure correlation functions will see deviations from the equilibrium ones in their spectral densities. 
Solving for various two-point functions in the background with falling shell/membrane we found such 
deviations: they are oscillating in frequency or peak at certain "echo" times : see more on this interesting 
phenomenon in |191j . 

8 Brief summary and two open questions 

As the reader have seen from this review, the field is actively developing and is very much in flux. At 
such stage it would be dangerous to make any firm conclusions: instead let me explain in general terms 
where we are now. 

The near-Tc temperature region got a lot of new attention: from old and experimentally discredited 
view of being a "mixed phase" - a mixture of dense QGP and much more dilute hadronic phase - it 
got a new name and picture, it is now viewed as a magnetic plasma region. In it monopoles already 
dominate the bulk of plasma and expelled electric fields/objects, but they are not yet Bose condensed, 
as in the confined phase. Its theory is rapidly developing, both based on classical MD and quantum 
phenomena eventually leading to BEG of monopoles and confinement. 

The " quasi- conformal region" T > 2Tc would perhaps be amenable to AdS/CFT approach, on which 
we have many intriguing results described above. In particular we have gained amazing insights into 
the nature of hydrodynamics (deriving it directly from Einstein eqn, order by order in derivatives) and 
dissipation in general, relating those with classical information loss into the black holes. 

Here comes the main open experimental question, originating from the following fact: the RHIC 
experiments had never ventured into the " quasi- conformal region" . Indeed one needs to run heavy ion 
collision at LHC to do that. Thus all of us are waiting for the first heavy ion run and the first ALICE 
data, which will tell us if hydrodynamics (elliptic flow especially) will be as good there as it was at 
RHIC. The positive answer will put AdS/CFT - and our relations with string community - to further 
heavy use. The negative answer would mean that the "perfect liquid" at RHIC is a near-Tc phenomenon, 
perhaps induced by an interplay of electric and magnetic quasiparticles we discussed above. 

Let me end with the central theoretical question of the day: is there a direct relation between two 
explanations of "perfect liquid", the one based on electric/magnetic duality /transition and that based 
on AdS/CFT duality? In essence, the former is based on the gauge coupling g'^/4:7f passing from the 
value < 1 to "strong coupling" or magnetic domain where it is > 1. AdS/CFT duality operates with 
the 't Hooft coupling A = g'^Nc- Since Nc is thought of as infinitely large, A is large always and the 
value of g itself is of no importance. Whether we are in electric or magnetic domain does not really 
matter! What this all means is that a progress in AdS/CFT at finite N,, is badly needed, before we can 
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unite the two main theory directions - based on two famous duahties - into one consistent theory. 

9 Acknowledgements 

The work was partially done during my participation in Galileo Galilei Institute spring 2008 program, 
and I am grateful to the Institute for support and to my fellow co-organizers, especially to Adriano 
Di Giacomo and to Valentin Zakharov, who made this interesting program happen. Many new things 
I learned from its participants are embedded in this review. Many more people helped me, especially 
my collaborators at Stony Brook Ismail Zahed, Derek Teaney, Jinfeng Liao and Shu Li. The work is 
partially supported by the US-DOE grants DE-FG02-88ER40388 and DE-FG03-97ER4014. 

References 



[1] 

[2] 
[3 

[4] 
[5] 

[6] 

[7] 

[8 

[9 

[10 

[11 
[12 

[13 

[14 

[15 

[16 

[17 

[18 

[19 



J. C. Collins and M. J. Perry, Phys. Rev. Lett. 34, 1353 (1975). 

E. V. Shuryak, Sov. Phys. JETP 47, 212 (1978) [Zh. Eksp. Teor. Fiz. 74, 408 (1978)]. 

E. V. Shuryak, Phys. Lett. B 78, 150 (1978) [Sov. J. Nucl. Phys. 28, 408.1978 YAFIA,28,796 (1978 
YAFIA,28,796-808.1978)]. 

J. O. Andersen and M. Strickland, Annals Phys. 317, 281 (2005) |arXiv:hep-ph/0404164 1 . 



B. A. Gelm an, E. V. Shuryak and I. Zahed, Phys. Rev. C 74, 044908 (2006) | arXiv:nucl- 
th/0601029| . 



S. Mrowczynski and M. H. Thoma, Ann. Rev. Nucl. Part. Sci. 57, 61 (2007) [arXivinucT 
[th707 01002|. 

K.M. O'Hara, S.L. Hemmer, M.E. Gehm, S.R. Granade, and J.E. Thomas, Science 298 (2002) 
2179. 



B. A. Gelman, E. V. Shuryak and I. Zahed, arXiv:nucl-th/0410067 

A. Turlapov, J. Kinast, B. Clancy, Le Luo, J. Joseph, J. E. Thomas, Kovtun:2005ev 

E. V. Shuryak, |arXiv:nucl-th/0606046l 



D. T. Son, arXiv: 0804. 3972 [hep-th] 



K. Balasubramanian and J. McGreevy, arXiv:0804.4053 [hep-th]. 

E. V. Shuryak, Phys. Lett. B 34, 509 (1971). 

E. V. Shuryak and O. V. Zhirov, Phys. Lett. B 89, 253 (1979). 

L. D. Landau, Izv. Akad. Nauk Ser. Fiz. 17, 51 (1953). 

E. Fermi, Phys. Rev. 81, 683 (1951). 

I. Y. Pomeranchuk, Dokl. Akad. Nauk Ser. Fiz. 78, 889 (1951). 

J. Bjorken, Phys. Rev. D27(1983)140 



D. Teaney, J. Lauret and E. V. Shuryak, arXiv:nucl-th/0110037 

69 



[20] C. B. Chin, E. C. G. Sudarshan and K. H. Wang, Phys. Rev 012(1975)902 

[21] M. I. Gorenshtein, V. A. Zhdanov and Yu. M. Sinjukov, ZhETF 74(1978)833 

[22] E. V. Shuryak, Yadernaya Fizika 16(1972)395 

[23] C. M. Hung and E. V. Shuryak, Phys. Rev. C 57, 1891 (1998) |arXiv:hep-ph/9709264| . 

[24] T. Hirano, Acta Phys. Polon. B 36, 187 (2005) |arXiv:nucl-th/0410017| . 

[25] C. Nonaka and S. A. Bass, Phys. Rev. C 75, 014902 (2007) fa rXiv:nucl-th/0607018| . 



[26] K. Adcox et al. [PHENIX Collaboration], Nucl. Phys. A 757, 184 (2005) |arXiv:nucl-ex/0410003|. 



[27] S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. Lett. 91, 182301 (2003) |arXiv:nucl- 
|ex/030501^ . 



[28] C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. 87, 182301 (2001) |arXiv:nucl-ex/0107003] . 
[29] T. Hirano and K. Tsuda, Phys. Rev. C 66, 054905 (2002) |arXiv:nucl-th/020 5043] . 
[30] P. F. Kolb and U. W. Heinz, |arXiv:nucl-th /0305084l 



[31] P. Huovinen, P. F. Kolb, U. W. Heinz, P. V. Ruuskanen and S. A. Voloshin, Phys. Lett. B 503, 
58 (2001) [arXiv:hep-ph/0101136] . 

[32] S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. C 69, 034909 (2004) |arXiv:nucl-ex/0307022| . 

[33] P. F. Kolb, J. Sollfrank and U. W. Heinz, Phys. Rev. C 62, 054909 (2000) |arXiv:hep-ph7 0006129] . 



[34] B. B. Back et al. [PHOBOS Collaboration], Phys. Rev. C 72, 051901 (2005) [arXivinucl- 
|ex/0407012] . 



[35] D. Teaney, Phys. Rev. C 68, 034913 (2003) |arXiv:nucl-th/0301099| . 

[36] P. Romatschke and U. Romatschke, Phys. Rev. Lett. 99, 172301 (2007) |arXiv:0706. 15221 [nucl-th]] . 

[37] K. Dusling and D. Teaney, Phys. Rev. C 77, 034905 (2008) [arXiv:0710.5932. [nucl-th]]. 



[38] U. W. Heinz and H. Song, arXiv:0806.0352 [nucl-th]. 



[39] D. Molnar, |^rXiv:0707.12"5r [nucl-th]. 



[40] T. Lappi and R. Venugopalan, Phys. Rev. C 74, 054905 (2006) | arXiv:nucl-th/060902l] . 
[41] F. Gelis and R. Venugopalan, Acta Phys. Polon. B 37, 3253 (2006). 



[42] H. B. Meyer, Nucl. Phys. B 795, 230 (2008) |arXiv:07ii:0738| [hep-lat]] . 



[43] R. Baier, Y. L. Dokshitzer, A. H. Mueller, S. Peigne and D. Schiff, Nucl. Phys. B 484, 265 (1997) 
[arXiv:hep-ph/9608322]. 



[44] E. V. Shuryak and I. Zahed, Phys. Rev. D 67, 054025 (2003) |arXiv:hep-ph/0207163 



[45] D. E. Kharzeev, arXiv:0806.0358 [hep-ph]. 



[46] E. V. Shuryak, Phys. Rev. C 66, 027902 (2002) |arXiv:nucl-th/01 12042 . 

70 



[47 
[48 

[49 
[50 

[51 
[52 
[53 
[54 
[55 

[56 

[57 

[58; 

[59 
[60 
[61 

[62; 

[63 
[64 
[65 

[66 

[67 
[68 
[69 
[70 



G. D. Moore and D. Teaney, Phys. Rev. C 71, 064904 (2005) |arXiv:hep-ph/0412346| . 

H. van Hees, M . Mannarelli, V. Greco and R. Rapp, Phys. Rev. Lett. 100, 192301 (2008) 
arXiv:0709.2884| [hep-ph]]. 



H. Stoecker, Nucl. Phys. A 750, 121 (2005) |arXiv:nucl-th /0406018|. 

J. Gasalderrey-Solana, E. V. Shuryak and D. Teaney, J. Phys. Conf. Ser. 27, 22 (2005) [Nucl. Phys. 
A 774, 577 (2006)] | arXiv:hep-ph/0411315] . 

B. Bauchle, L. Csernai and H. Stocker, |arXiv:0710.1476 [nucl-th]. 



J. Casalderrey-Solana and E. V. Shuryak, arXiv:hep-ph/05 11263 



F. Antinori and E. V. Shuryak, J. Phys. G 31, L19 (2005) |arXiv:nucl-th/0507046| . 



D. Kharzeev and H. Satz, Phys. Lett. B 356, 365 (1995) | arXiv:hep-ph/9504397| . 

A. Andronic, P. Braun-Munzinger, K. Redlich and J. Stachel, Phys. Lett. B 571, 36 (2003) 
|arXiv:nucl-th/03 03036]. 

C. Young and E. Shuryak, a rXiv:0803.2866 [nucl-th]. 

T. Matsui and H. Satz, Phys. Lett. B 178, 416 (1986). 

M. Asakawa and T. Hatsuda, Phys. Rev. Lett. 92, 012001 (2004) |arXiv:hep-lat/ 0308034| . 

S. Datta, F. Karsch, P . Petreczky and L Wetzorke, Nucl. Phys. Proc. Suppl. 119, 487 (2003) 
arXiv:hep-lat/0208012| . 



A. Mocsy, P. Pet reczky and J. C asalderrey-Solana, Nucl. Phys. A 783, 485 (2007) [Nucl. Phys. A 
785, 266 (2007)] |arXiv:hep-ph/0609205|. 



A. Mocsy and P. Petreczky, PoS LAT2007, 216 (2007) [arXiv:0710.5205 [hep-lat]]. 

R. L. Thews, Nucl. Phys. A 783, 301 (2007) |arXiv:hep-ph /0609121|. 

L. Grandchamp, R. Rapp and G. E. Brown, J. Phys. G 30, S1355 (2004) |arXiv:hep-ph/0403204 |. 

F. Karsch, D. Kharzeev and H. Satz, Phys. Lett. B 637, 75 (2006) fa rXiv:hep-ph/0512239l . 

G. Cossu, M. D'Eha, A. Di Giacomo, G. Lacagnina and C. Pica, Phys. Rev. D 77, 074506 (2008) 
|arXiv:0802.1795 [hep-lat]]. 



M. A. Ste phanov, K. Rajagopal and E. V. Shuryak, Phys. Rev. Lett. 81, 4816 (1998) |arXiv:hep 
|ph/9806219]. 



L. McLerran and R. D. Pisarski, Nucl. Phys. A 796, 83 (2007) |arXiv: 0706. 2191 [hep-ph]]. 



L. Y. Glozman, |arXiv:0803.1636| [hep-ph]. 

J. Liao and E. V. Shuryak, Phys. Rev. D 73, 014509 (2006) |arXiv:hep-ph/0510110| . 

L. Castillejo, P. S. J. Jones, A. D. Jackson, J. J. M. Verbaarschot and A. Jackson, Nucl. Phys. A 
501, 801 (1989). 



[71] R. Rapp, E. V. Shuryak and I. Zahed, Phys. Rev. D 63, 034008 (2001) |arXiv:hep-ph/0008207] . 

71 



[72] M. G. Alford, A. Schmitt, K. Rajagopal and T. Schafer, [arXiv:0709.4635l [hep-ph]. 



[73] F. Karsch [RBC Collaboration and HotQCD Collaboration], |arXiv:0804.4148 [hep-lat]. 



[74] M. N. Chernodub, K. Ishiguro, A. Nakamura, T. Sekido, T. Suzuki and V. I. Zakharov, PoS 
LAT2007, 174 (2007) |arXiv :0710.2547| [hep-lat]]. 

[75] A. M. Polyakov, Phys. Lett. B 72, 477 (1978). 

[76] A. Nakamura, T. Saito and S. Sakai, Phys. Rev. D 69, 014506 (2004) |arXiv:hep-lat/0311024| . 

[77] O. Kaczmarek and F. Zantow, PoS LAT2005, 192 (2006) |arXiv:hep-lat/05 10094] . 



78] G. S. Bah, arXiv:hep-ph/9809351 



[79] O. Kaczmarek and F. Zantow, Phys. Rev. D 71, 114510 (2005) |arXiv:hep-la t/0503017| . 



[80] P. Petreczky and K. Petrov, Phys. Rev. D 70, 054503 (2004) |arXiv:hep-lat/0405009] . 
il] D. Antonov, S. Domdey and H. J. Pirner, Nucl. Phys. A 789, 357 (2007) |arXiv:hep-ph/0612256| . 



[82] E. V. Shuryak and I. Zahed, Phys. Rev. C 70, 021901 (2004) jarXiv:hep-ph/ 0307267|. 



[83] E. V. Shuryak and I. Zahed, Phys. Rev. D 70, 054507 (2004) |arXiv:hep-ph/0403127] . 

[84] L.D. Landau and E.M.Lifshits, Quantum Mechanics, chapter 90. 

[85] W. M. Alberico, A. B eraudo, A. De Pace and A. Molinari, Phys. Rev. D 75, 074009 (2007) 
farXi v:hep-ph/06120'62] . 

[86] J. Liao and E. V. Shuryak, Nucl. Phys. A 775, 224 (2006) |arXiv:hep-ph/0508035| . 

57] A. V. Smilga, Annals Phys. 234, 1 (1994). 

P. de Forcrand and O. Jahn, Nucl. Phys. B 651, 125 (2003) |arXiv:hep-lat/0211004| . 



[89] M. Unsal and L. G. Yaffe, arXiv: 0803. 0344 [hep-th]. 



[90] M. Shifman and M. Unsal, arXiv:0802.1232 [hep-th]. 



[91] E. M. Ilgenfritz and E. V. Shuryak, Nucl. Phys. B 319, 511 (1989). 



[92] T. Schafer and E. V. Shuryak, Rev. Mod. Phys. 70, 323 (1998) |arXiv:hep-ph/96 10451] . 

[93] A. M. Polyakov, Nucl. Phys. B 120, 429 (1977). 

[94] M. Unsal, [arXiv:0709.3269 [hep-th]. 

[95] D. J. Gross, R. D. Pisarski and L. G. Yaffe, Rev. Mod. Phys. 53, 43 (1981). 

[96] R. D. Pisarski, Phys. Rev. D 74, 121703 (2006) |arXiv:hep-ph/0608242j. 



[97] A. Bazavov, B. A. Berg and A. Dumitru, |arXiv:Q805.0784| [hep-ph]. 



[98] K. Fukushima, Phys. Lett. B 591, 277 (2004) |arXiv:hep-ph/0310121] 



[99] E. Megias, E. Ruiz Arriola and L. L. Salcedo, Phys. Rev. D 74, 065005 (2006) | arXiv:hep- 
,ph/0412308] . 

72 



[100] C. Ratti, S. Roessner, M. A. Thaler and W. Weise, Eur. Phys. J. C 49, 213 (2007) |a rXiv:hep- 
iph/0609218] . 

C. R. Allton et al, Phys. Rev. D 71, 054508 (2005) |arXiv:hep-la t/0501030]. 

Y. M. Shnir, Berlin, Germany: Springer (2005) 532 p 

S. Mandelstam, Phys. Rept. 23, 245 (1976). G. 't Hooft, Nucl. Phys. B 190, 455 (1981). 

N. Seiberg and E. W itten, Nud. Phys. B 426, 19 (1994) [Erratum-ibid. B 430, 485 (1994)] 
|arXiv:hep-th/94070'87] . 



101 
102 
103 
104 

105 
106 
107 

108 
109 
110 
111 
112 
113 
114 
115 
116 
117 
118 
119 
120 

121 

122 
123 
124 
125 
126 



D. Tong, arXiv:hep-th/0509216 



M. Shifman and A. Yung, Phys. Rev. D 77, 125016 (2008) |arXiv: 0803. 0158 [hep-th]]. 



M. Shifman and A. Yung, Prepared for 7th Workshop on Continuous Advances in QCD, Min- 
neapolis, Minnesota, 11- 14 May 2006 

T. A. DeGrand and D. Toussaint, Phys. Rev. D 22, 2478 (1980). 



A. D'Alessandro and M. D'Eha, arXiv:0711.1266 [hep-lat]. 



H. B. Meyer, arXiv: 0808. 1950 [hep-lat]. 



T. C. Kraan and P. van Baal, Nucl. Phys. B 533, 627 (1998) |arXiv:hep-th/9805168| . 



D. Diakonov and N. Gromov, Phys. Rev. D 72, 025003 (2005) |arXiv:hep-th/0502132]. 



D. Diakonov and V. Petrov, Phys. Rev. D 76, 056001 (2007) |arXiv:0704.3181| [hep-th]]. 

J. Liao and E. Shuryak, Phys. Rev. C 75, 054907 (2007) |arXiv:hep-ph/0611131| . 

M. N. Chernodub and V. I. Zakharov, Phys. Rev. Lett. 98, 082002 (2007) |arXiv:hep-ph/0611228| . 



J. Liao and E. Shuryak, Phys. Rev. C 77, 064905 (2008) |arXiv:0706.4465 [hep-ph]]. 



J. Liao and E. Shuryak, arXiv:0804.4890, [hep-ph]. 



S. C. Huot, S. Jeon and G. D. Moore, Phys. Rev. Lett. 98, 172303 (2007) |arXiv:hep-ph/0608062| . 
P. M. Chesler and A. Vuorinen, JHEP 0611, 037 (2006) |arXiv:hep-ph/0607148| . 



G. Pohcastro, D. T. Son and A. O. Starinets, Phys. Rev. Lett. 87, 081601 (2001) |arXiv:hep- 
pi/0104066 |. 



J. Casalderrey-Solana and D. Teaney, |hep-ph/0605199 



K. Dusling and C. Young, |arXrv:0707 .2068 [nucl-th]. 

R.P. Feynman, Phys. Rev. 90, 1116 (1953) 

R.P. Feynman, Phys. Rev. 91, 1291 (1953) 

D.M.Ceperley, RMP 67 , 279, (1995). 

S. Jang, S. Jang and J.A. Voth, J. Chem. Phys. 115, 7832 (2001) 

73 



127; 

128 
129 
130 

131 
132 

133 

134 

135 

136 

137 
138 

139 
140 

141 
142 
143 

144 
145 
146 
147 
148 



R. Kikuchi et al., Phys. Rev. 119, 1823 (1960) 

R.A. Aziz, et al. J. Chem. Phys. 70, 4330 (1979) 

V. Elser, PhD. Thesis, U.C. Berkeley, (1984) 

M.Cristoforetti and E.Shuryak, Bose condensation in liquids and confinement of monopoles, in 
progress 

A. Gorsky and V. Zakharov, Phys. Rev. D 77, 045017 (2008) |arXiv: 0707. 12841 [hep-th]]. 



L. Del Debbio, M. Faber, J. Greensite and S. Olejnik, Phys. Rev. D 55, 2298 (1997) | arXiv:hep- 
lat/9610005j . 



P. Petreczky, F. Karsch, E. Laermann, S. Stickan and I. Wetzorke, Nucl. Phys. Proc. Suppl. 106, 
513 (2002) | arXiv:hep-lat/0110111] . 

J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998) [Int. J. Theor. Phys. 38, 1113 (1999)] 



arXiv:hep-th/9711200] . 



J. M. Maldacena, Phys. Rev. Lett. 80, 4859 (1998) |ar Xiv:hep-th/9 803002|. S. J. Rey and J. T. Yee, 
Eur. Phys. J. C 22, 379 (2001) ||arXiv:hep- th/9803001]. 

K. S. . Thorne, R. H. . Price and D. A. . Macdonald, NEW HAVEN, USA: YALE UNIV. PR. 
(1986) 367p 



M. Parikh and F. Wilczek, Phys. Rev. D 58, 064011 (1998) | arXiv:gr-qc/9712077| . 

S. S. Gu bser, I. R. Klebanov and A. A. Tseytlin, Nucl. Phys. B 534, 202 (1998) |arXiv:hep- 



th/9805156i . 



E. Shuryak, 'i^iv:0711.0004| [hep-ph]. 



A. Karch, D. T. Son, E. Katz and M. A. Stephanov, Prepared for 7th Workshop on Continuous 
Advances in QCD, Minneapolis, Minnesota, 11- 14 May 2006 



U. Gursoy and E. Kiritsis, JHEP 0802, 032 (2008) | arXiv:0707.1324| [hep-th]] . 

E. Shuryak and I. Zahed, Phys. Rev. D 69, 046005 (2004) |arXiv:hep-th/0308073| . 



I. R. Klebanov, J. M. Maldacena and C. B. . Thorn, JHEP 0604, 024 (2006) |arXiv:hep- 
[th/06022 55|. 



J. M. Maldacena, Phys. Rev. Lett. 80, 4859 (1998) | arXiv:hep-th/9803002| . 

S. J. Rey and J. T. Yee, Eur. Phys. J. C 22, 379 (2001) [ arXiv:hep-th/9803001| . 



S. Lin and E. Shuryak, Phys. Rev. D 76, 085014 (2007) | arXiv:0707.3135, [hep-th]]. 

K. Peeters, J. Sonnenschein and M. Zamaklar, JHEP 0602, 009 (2006) |arXiv: hep-th/05 1 1 044] . 

G. W. Se menoff and K. Zarembo, Nucl. Phys. Proc. Suppl. 108, 106 (2002) |arXiv:hep- 



th/0202156|. 



[149] A. Mikhailov, arXiv:hep-th/0305196 



74 



150 
151 

152 
153 
154 
155 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 
168 
169 

170 

171 
172 
173 
174 
175 



S. J. Sin and I. Zahed, Phys. Lett. B 608, 265 (2005) |arXiv:hep-th/0407215| . 

C. P. Herzo g, A. Karch, P. Kovtun, C. Kozcaz and L. G. Yaffe, JHEP 0607, 013 (2006) |arXiv:hep- 
th/0605158| . 



S. S. Gubser, Phys. Rev. D 74, 126005 (2006) |arXiv:hep-th70605182] . 



A. Buchel, Phys. Rev. D 74, 046006 (2006) |a rXiv:hep-th/0605178| . 

S. J. Sin and I. Zahed, Phys. Lett. B 648, 318 (2007) |arXiv:hep-ph /0606049|. 



J. J. Friess, S. S. Gubser and G. Michalogiorgakis, JHEP 0609, 072 (2006) [arXiv:hep-th/0605292] 



P. M. Chesler and L. G. Yaffe, j|arXiv:0712.0050 [hep-th]. 



S. S. Gubser and A. Yarom, arXiv:0803.008l1 [hep-th]. 



M. Lubhnsky and E. Shuryak, Phys. Rev. C 76, 021901 (2007) | arXiv:0704.1647| [hep-ph]]. 
G. Pohcastro, D. T. Son and A. O. Starinets, JHEP 0209, 043 (2002) |arXiv:hep-th/0205052l . 
P. K. Kovtun and A. O. Starinets, Phys. Rev. D 72, 086009 (2005) |arXiv:hep-th/0506184| . 



D. Teaney, Phys. Rev. D 74, 045025 (2006) | arXiv:hep-ph/0602044| . 

P. Kovtun and A. Starinets, Phys. Rev. Lett. 96, 131601 (2006) |arXiv:hep- th/0602059|. 

P. K. Kovtun and A. O. Starinets, Phys. Rev. D 72, 086009 (2005) |arXiv:hep-th/0506184 |. 



G. Siopsis, arXiv: 0804. 2713 [hep-th]. 



D. Kharzeev and K. Tuchin, [arXiv: 0705. 4280 [hep-ph]. 



F. Karsch, D. Kharzeev and K. Tuchin, Phys. Lett. B 663, 217 (2008) |arXiv:0711.0914| [hep-ph]]. 
A. Buchel, Phys. Rev. D 72, 106002 (2005) |arXiv:hep-th/0509083| . 



A. Parnachev and A. Starinets, JHEP 0510, 027 (2005) |arXiv:hep-th/0506144|. 



P. Benincasa, A. Buchel and A. O. Starinets, Nucl. Phys. B 733, 160 (2006) [arXivihep- 
[th/0507026 |. 



U. Gursoy, E. Kiritsis, L. Mazzanti and F. Nitti, arXiv: 0804. 0899 [hep-th] 



S. S. Gubser, A. Nellore, S. S. Pufu and F. D. Rocha, arXiv:0804.1950 [hep-th]. 



S. S. Gubser, S. S. Pufu and F. D. Rocha, arXiv:0806.0407, [hep-th]. 



R. A. Janik and R. B. Peschanski, Phys. Rev. D 73, 045013 (2006) [ arXiv:hep-th/0512162 

S. Nakamura and S. J. Sin, JHEP 0609, 020 (2006) ) arXiv:hep-th/0607123| . 

P. Benincasa, A. Buchel, M. P. Heller and R. A. Janik, Phys. Rev. D 77, 046006 (2008) 
|arXiv:0712.2025 [hep-th]]. 



[176] M. P. Heller, P. Surowka, R. Loganayagam, M. Spalinski and S. E. Vazquez, arXiv:0805.3774 
[hep-th] . 

75 



[177' 

[178 

[179 
[180 
[181 
[182 
[183 

[184 
[185 
[186 
[187 

[188 
[189 
[190 
[191 
[192 
[193 



R. Baier, P. Romatschke, D. T. Son, A. O. Starinets and M. A. Stephanov, JHEP 0804, 100 

(2008) |a rXiv:0712.2451 [hep-th]]. 

S. Bhattacharyya, V. E. Hubeny, S. Minwalla and M. Rangamani, JHEP 0802, 045 (2008) 
^arXiv:0712.2456 [hep-th]]. 



S. Bhattacharyya et al, arXiv: 0803. 25261 [hep-th]. 



M. Natsuume and T. Okamura, Phys. Rev. D 77, 066014 (2008) |arXiv:0712.2916, [hep-th]]. 



M. Haack and A. Yarom, arXiv:0806.4602 [hep-th]. 



E. Shuryak, S. J. Sin and I. Zahed, J. Korean Phys. Soc. 50, 384 (2007) |arXiv:hep-th/0511199 



K. Kajantie, J. Louko and T. Tahkokalho, Phys. Rev. D 76, 106006 (2007) [arXiv: 0705. 1791 

[hep-th]]. 

B. Andersson, G. Gustafson, G. Ingelman and T. Sjostrand, Phys. Rept. 97, 31 (1983). 



S. Lin and E. Shuryak, Phys. Rev. D 77, 085013 (2008) |arXiv:hep-ph/0610168] 



D. M. Hofman and J. Maldacena, JHEP 0805, 012 (2008) | arXiv:0803.1467| [hep-th]]. 

D. J. Gross, A. Mikhailov and R. Roiban, Annals Phys. 301, 31 (2002) farXiv:hep-th/0205066l. 
Y. Makeenko, P. Olesen and G. W. Semenoff, Nucl. Phys. B 748, 170 (2006) [arXiv:hep-, 
[th/0602100] . 



D. Grumiller and P. Romatschke, |arXiv:0803.3226 [hep-th] 



S. S. Gubser, S. S. Pufu and A. Yarom, arXiv: 0805.1 55 1| [hep-th]. 



C. Ratti, S. Roessner and W. Weise, Phys. Lett. B 649, 57 (2007) |arXiv:hep-ph/0701091] . 



S. Lin and E. Shuryak, WXiv:0808.0910| [hep-th]. 



B. Betz, M. Gyulassy, J. Noronha and G. Torrieri, arXiv:0807.4526 [hep-ph]. 



M. Gyulassy, J. Noronha and G. Torrieri, |arXiv:0807.2235, [hep-ph]. 



76 



